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FOREWORD 


Within  the  past  five  years,  several  major  engineering  excavation 
projects  have  been  accomplished  with  the  use  of  large  burled  explosive 
charges.  The  U.S.  Army  Corps  of  Engineers  Waterways  Experiment  Station 
Explosive  Excavation  Research  Laboratory  (EERL)  has  been  Intimately 
Involved  In  the  development  and  application  of  this  technology.  The  use 
of  large  explosive  cratering  charges  on  dry  land  Is  now  relatively  well 
understood  and  detailed  design  procedures  are  available.  However,  a  number 
of  potential  applications  Involve  excavation  In  submerged  environments. 
Underwater  cratering  dynamics  are  not  well  understood.  Accordingly,  It 
was  decided  to  Initiate  a  program  of  numerical  hydrodynamic  calculations 
to  model  the  effects  of  burled  underwater  explosions  and  achieve  a  better 
understanding  of  the  uynamlc  mechanisms  Involved.  This  program  was  Imple¬ 
mented  under  the  direction  of  EERL,  and  calculations  were  performed  by  the 
Earth  Sciences  (K)  Division  of  the  Lawrence  Livermore  Laboratory  (LLL). 

This  report  sunmarlzes  the  results  of  18  one-dimensional  hydrodynamic 
calculations  plus  several  subsidiary  comparison  calculations  and  acoustic 
wave  calculations  for  high-explosive  detonations  in  saturated  media  and  In 
saturated  media  overlain  by  a  layer  of  water.  The  hydrodynamic  calcu¬ 
lations  provide  an  early-time  description  of  close-range  shock  wave  trans¬ 
mission  within  the  rock  medium,  Including  time  histories  of  pressure, 
stress,  and  particle  velocity,  shock  Interactions  at  the  Interfaces,  spall 
velocities,  shear  stress  and  failure  effects,  etc.  Safety-related  pheno¬ 
mena  may  likewise  be  estimated  from  this  Information.  These  calculatlonal 
results  form  a  basis  for  future  two-dimensional  multilayer  cratering 
calculations  which  may  be  used  to  design  underwater  excavation  events,  both 
nuclear  and  conventional. 

The  Investigation  described  herein  was  sponsored  by  U.S.  Army  Research 
Office  grant  ARO-D  No.  I-479-E  under  R  &  D  Project  No.  20061 102B33G.  In 
accordance  with  ARO-D  program  requirements,  this  report  documents  results 
of  the  one-dimensional  underwater  calculations  and  related  studies  to  date. 
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ONE-DIMENSIONAL  ANALYSES:  SHOCK  WAVE  PROPAGATION 
FROM  UNDERWATER  CRATERING  DETONATIONS 


SUMMARY 

Research  has  recently  been  conducted  on  the  use  of  large  buried 
chemical  explosive  charges  for  engineering  excavation.  Investigations 
in  this  area  were  Initially  concerned  with  crater  formation  in  relatively 
homogeneous  rock  media,  a  process  which  is  now  well  understood. 

However,  it  has  become  apparent  that  the  greatest  cost  and  operational 
advantages  may  be  realized  by  applying  the  technique  to  rock  excavation 
In  an  underwater  environment.  Cratering  dynamics  In  an  underwater  (two- 
layer)  configuration  have  not  been  well  understood.  Small-scale 
modeling  tests  and  large-scale  excavation  projects  have  revealed  two 
new  factors  which  significantly  Influence  underwater  cratering  processes: 
(1)  Early-time  dynamic  effects  caused  by  the  presence  of  the  rock-water 
interface  and  water  layer;  and  (2)  very  late-time  water  washback  and 
slope  failure  effects  in  the  crater  vicinity.  This  report  addresses  the 
first  of  these  two  effects  by  means  of  hydrodynamic  computer  calcula¬ 
tions.  The  one-dimensional  spherically-symmetric  "SOC"  computer  program 
has  been  used  to  calculate  the  dynamic  effects  of  a  spherical  explosive 
charge  emplaced  in  rock  media  and  in  rock  with  overlying  water  layers  of 
various  depths.  Calculations  have  been  conducted  for  three  different 
rock  types,  coral,  weak  basalt,  and  granite.  Significant  correlations 
of  dynamic  behavior  with  medium  properties  and  layer  depths  have  been 
established.  Some  comparisons  with  experimental  data  have  also  been 
accomplished.  The  results  indicate  that  computer  design  calculations 
for  future  underwater  events  are  quite  feasible. 
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ONE  DIMENSIONAL  ANALYSES: 

SHOCK  WAVE  PROPAGATION  FROM 
UNDERWATER  CRATERING  DETONATIONS 

Charles  M.  Snell 
I.  INTRODUCTION 

Widespread  Interest  In  underwater  excavation  with  large  explosive 
charges  was  first  manifested  during  the  early  stages  of  the  Plowshare 
program  for  peaceful  uses  of  nuclear  explosives.  ^  This  Interest  was 
primarily  a  result  of  three  factors:  (1)  The  enormous  commercial  value 
of  deep-draft  harbors  and  navigation  facilities;  (2)  The  difficulty  and 
expense  of  underwater  excavation  In  rock  with  conventional  techniques; 
and  (3)  The  apparent  potential  of  nuclear  explosives  for  rapid,  effi¬ 
cient  large-scale  excavation.  A  variety  of  possible  sites  and  applica¬ 
tions  were  considered.  Including  some  projects  which  could  not  be 
economically  accomplished  with  conventional  construction  techniques.  An 
experimental  harbor  excavation  model.  Project  Chariot,  was  planned  for 
a  remote  section  of  the  Alaskan  coast,  but  was  never  Implemented. 

Other  conceptual  designs  Included  the  excavation  of  navigation  channels, 
widening  and  deepening  extant  channels,  removal  of  submerged  obstacles, 
creation  or  Improvement  of  harbor  entranceways ,  etc. 

During  the  same  period,  extensive  Investigations  of  dry-land 
cratering  with  nuclear  explosives  were  carried  out.  Both  theoretical 
and  experimental  studies  were  performed,  culminating  In  the  series  of 
Plowshare  cratering  events  at  the  Nevada  Test  Site.  A  physical  model 
was  developed  to  describe  the  phenomena  of  crater  formation.  Hydro- 
dynamic  computer  programs  were  written  to  predict  effects  of  burled 
detonations.  It  soon  became  apparent  that  computer  techniques  offered  a 


-2- 


means  of  analyzing  the  results  of  past  experiments  and  of  predicting  both 
crater  dimensions  and  other  physical  effects  of  future  events.  Thus,  the 
design  uncertainties  and  the  need  for  extensive  field  experiments  could  be 
greatly  reduced.  Successful  calculations  were  performed  for  the  Plowshare 
events  and  for  high-explosive  cratering  tests  In  various  media.  Certain  of 
the  calculations  made  during  the  Trans-Isthmian  Canal  Studies  Involved  wet  or 
fully  saturated  rock,  similar  to  the  media  which  might  be  encountered  In 
underwater  excavation  projects.  However,  no  studies  were  performed  for  sub¬ 
merged  configurations. 

After  conclusion  of  the  Plowshare  test  series.  Interest  In  nuclear 
cratering  applications  waned.  Nonetheless,  some  of  the  modeling  tests 
associated  with  the  Plowshare  program  demonstrated  that  high-explosive  crater¬ 
ing  charges  present  certain  benefits  for  smaller  excavation  tasks.  Within 
the  past  five  years,  several  major  engineering  projects  have  been  performed 
with  large  conventional  explosive  emplacements.^5,6^  Two  of  these  projects 
Involved  excavation  In  an  underwater  environment.^-7,8,9^  Due  to  the  high 
cost  of  excavation  In  submerged  rock,  this  area  Is  perhaps  the  most  promising 
for  future  use  of  large-charge  blasting.  In  addition  to  the  cost  advantage, 
explosive  excavation  may  also  provide  environmental  benefits  over  other 
methods,  greater  speed  and  efficiency  of  construction,  applicability  to 
emergency  situations,  and  the  ability  to  carry  out  projects  which  would  be 
extremely  difficult  or  dangerous  to  accomplish  with  standard  techniques.  The 
major  problem  associated  with  underwater  designs  Is  the  lack  of  experience 
and  empirical  data  on  underwater  cratering.  The  earliest  submerged  high- 
explosive  cratering  tests  were  Intended  to  model  weapons  effects;  the  charges 
were  located  at  the  water-medium  Interface^10^,  or  were  burled  at  very  shallow 
depths  In  the  medium  (too  shallow  to  be  of  Interest  for  excavation  design). 
More  recent  excavation  shots  were  emplaced  In  saturated  sand^  or  In  weak  low- 
density  coral L  ,  with  no  large-scale  tests  In  submerged  rock  or  other  media. 
The  ability  to  accurately  predict  dry-land  crater  dimensions  In  a  variety 
of  media  Is  dependent  on  the  combination  of  a  broad  data  base  and  a  fundamental 


understanding  of  the  crater  formation  process.  Indeed,  no  major  surprises 
have  been  encountered  during  any  recent  experiments.  The  opposite  situation 
prevails  for  underwater  cratering:  Dimensions,  crater  shape,  and  collateral 
effects  cannot  be  reliably  predicted  before  the  shot.  In  the  past.  It  has 
been  necessary  to  conduct  expensive  single-charge  calibration  events  at  the 
project  site,  or  to  utilize  an  untested  design  which  may  not  precisely  achieve 
project  requirements. 

Most  underwater  cratering  emplacements  consist  essentially  of  a  two- 
layer  configuration,  dense  saturated  solid  material  overlain  by  a  much  less 
dense  layer  of  water.  The  considerable  differences  between  dry  land  and 
submerged  craters  are  thus  attributable  to  two  factors:  (1)  Water  washback 
or  liquefied  flow  of  the  fallback  material;  washback  effects  occur  at  very 
late  times,  after  the  ejecta  mound  has  collapsed.  And, (2)  the  Influence  of 
the  water  layer  on  shock  propagation  and  material  dynamics  at  early  times, 
as  the  shock  wave  travels  outward  and  Is  reflected  at  the  rock  Interface 
and  the  water  surface.  The  subsidiary  effects  of  saturation  may  also  reduce 
the  strength  of  the  rock  and  Increase  the  efficiency  of  shock  propagation 
(most  submerged  events  are  expected  to  occur  In  fully- saturated  material). 

Water  washback  effects  are  known  to  be  most  critical  In  fine-grained 
coheslonless  media  such  as  sand  or  crushed  coral.  For  small-scale  (s  one 
pound)  cratering  tests  In  saturated  sand,  It  has  been  observed  that  crater 
dimensions  and  shapes  change  dramatically  at  relatively  shallow  water  layer 
depths  between  0.05  and  0.2  times  the  charge  burial  depth  In  the  solid 
^Dv/Dr  *  °*05’0*2)*  Craters  In  this  depth  range  are  modified  from  typical 
"dry-land"  contours  Into  broad,  shallow,  saucer-shaped  depressions  with  little 
or  no  elevated  lip  material.  Appearances  Indicate  that  the  normal  lip  material 
Is  overtopped  by  water  washback;  slope  failure  occurs  and  ejected  material  flows 
Into  the  crater,  producing  the  flattened  shape.  The  crater  radius  actually 
Increases  sharply  at  0W/Dr  s  0.1,  while  the  crater  depth  and  lip  height  decrease 
abruptly.  For  deeper  water  layers,  Dw/Df  ~  0. 3-0.4,  the  crater  radius  begins 
to  decrease  gradually,  returning  to  approximately  the  "dry-land"  value  at 
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V°r  z  w^®  tht  crater  depth  also  continues  to  decrease  slowly.  Small 
lips  tend  to  reappear  for  deeper  water  layers,  but  the  shallow  crater  depths 
and  flattened  contours  are  maintained.^]  Current  evidence  Indicates  that 
shallow  flat-bottomed  craters  are  also  formed  by  large-scale  («  one-ton)  events 
In  a  submerged  sand  environment  (0^  <  0.6). ^  Finally,  somewhat  similar 

craters  are  produced  by  large  charges  (1-10  tons)  In  weak  submerged  coral, 
with  the  added  modification  that  collapse  and  compaction  of  the  porous  coral 
matrix  cause  very  large  crater  radii  and  very  large  volumes.  The  flattened 
shape  and  shallow,  crater  depth  are  again  apparent.  The  cratering  process 
In  coral  Is  Interpreted  as  one  of  shock  crushing  and  liquefaction  of  the  weak 
matrix  to  great  ranges,  beyond  the  distance  at  which  the  lip-forming  material 
would  normally  come  to  rest.  The  lip  fallback  and  Its  liquefied  nonsupporting 
base  thus  fall  radially  and  flow  Into  the  crater  with  the  water  washback. 

The  net  volume  of  the  coral  material  Is  reduced  by  fracturing,  compaction, 
and  later  settlement  of  the  finely-crushed  fallback,  accounting  for  the  rela¬ 
tively  large  volumes  and  radii  of  the  coral  craters.  It  Is  expected  that  flat- 
bottomed  llpless  craters  will  occur  In  coral  except  possibly  when  the  ratio  of 
water  depth  to  burial  depth  Is  extremely  small J7,5^ 

Very  weak  or  fine-grained  materials  are  the  exception  rather  than  the 
rule  In  most  underwater  blasting  situations.  Indeed,  the  greatest  advantages 
of  explosive  excavation  methods  may  be  attained  for  the  more  difficult  tasks 
of  excavating  Intermediate  or  high-strength  rock.  Unfortunately,  little 
Information  Is  available  on  underwater  cratering  In  media  of  this  type.  Small- 
scale  tests  (s  one  pound)  have  been  conducted  In  rather  weak  saturated  concrete 
with  no  overlying  water  and  with  water  depths  Dw/Df  s  0.3-2./.  This  material 
fractures  Into  large  chunks,  similar  to  rock  ejecta.  The  limited  data  Indicate 
that  crater  dimensions  are  not  sensitive  to  the  presence  of  water  overburden, 
although  the  crater  depths  may  be  slightly  reduced  at  all  water  layer  depths.^-13*7] 
It  appears  probable  that  dense  media  which  break  Into  large  granular  sizes  when 
cratered  are  not  subject  to  water  washback  and  slope  failure  effects,  at  least 
for  shallow  to  Intermediate  depth  water  layers.  The  conditions  under  which  slope 
failure  begins  will  depend  on  particle  size  (or  size  gradation),  particle 
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density*  water  depth*  and  the  exact  configuration  of  the  experiment.  The 
relative  Importance  of  these  factors  has  not  been  Investigated.  However* 
for  typical  high-explosive  cratering  configurations,  washback  effects 
are  likely  to  be  predominant  In  coheslonless  or  finely-crushed  materials 
(particle  sizes  *  a  f ew  Inches)*  and  much  less  Important  or  negligible  In 

most  rock  media  (particle  sizes  £  several  Inches  to  several  feet,  densities 
>  2  g/cm3). 

Since  engineering  experience  Is  not  currently  adequate  to  enable 
confident  design  of  submerged  events*  It  Is  appropriate  to  ask  If  the 
computational  methods  evolved  during  the  earlier  Plowshare  studies  may  be 
applied  to  this  area.  In  order  to  answer  this  question*  It  Is  necessary 
to  examine  briefly  the  characteristics  of  the  available  techniques.  Several 
versions  of  Lagranglan  hydrodynamic  computer  codes  -re  used  to  calculate  the 
effects  of  burled  detonations.  These  codes  utilize  finite  difference  schemes 
to  numerically  model  adiabatic  shock  wave  propagation  through  media  for  which 
constitutive  relations  are  known.  Elastic,  plastic.  Inelastic,  fluid,  and 
gaseous  behavior  of  a  material  may  be  taken  Into  account.  The  available 
codes  Include  SOC^**15^  (one-dimensional,  spherically  symmetric  shock  propa¬ 
gation)  and  TENSOR^®*^}  (two-dimensional,  cyllndrlcally-symmetrlc  shock 
propagation).  Both  codes  carry  up  to  26  variables  which  describe  the  state 
of  the  material  at  each  point  In  the  Lagranglan  grid.  The  basic  variables 
Include  displacements,  velocities.  Internal  energies,  pressures,  stress  tensor 
normal  components,  local  compression,  etc.  The  one-dimensional  code  may  be 
used  to  calculate  free-fleld  shock  propagation  at  a  constant  depth  In  a 
nearly- Infinite  medium  ("horizontal  problem");  alternatively,  vertical  shock 
propagation  In  a  constant  gravitational  field  may  be  modeled,  and  a  free 
surface  or  a  rigid  wall  may  be  Introduced  at  any  desired  range  to  reflect 
the  shock  wave  back  toward  the  center  of  the  problem  ("vertical  problem"). 

The  two-dimensional  code  Is  normally  used  to  determine  the  effects  of  a  single 
explosive  source  In  a  cyllndrlcally-symmetrlc  configuration.  For  example, 
a  spherical  source  emplaced  below  a  horizontal  free-surface  may 
be  used  to  model  the  formation  of  a  single-charge  crater.  Both  codes  have 
provision  for  regions  of  various  materials  with  different  properties. 


Certain  limitations  Inherent  In  the  Lagranglan  calculations  must 
be  considered.  The  codes  use  material  zones  of  finite  size  which  limit 
calculatlonal  resolution  of  "discontinuous"  phenomena  such  as  a  steep-fronted 
shock.  Resolution  may  be  Improved  by  finer  zoning,  at  the  cost  of  Increased 
calculation  time.  Likewise,  a  numerical  technique  known  as  linear  artificial 
viscosity,  used  to  damp  nonphysical  numerical  noise  In  the  problem,  may  some¬ 
times  cause  other  nonphysical  effects  in  the  solution.  The  influence  of  zone 
size  and  viscosity  parameters  may  sometimes  be  analyzed  and  taken  Into  account 
by  simple  one-dimensional  parameter  studies.  Another  limitation  of  the  hydro- 
dynamic  calculations  Is  related  to  the  time-scale  of  the  problem.  The  codes 
are  most  efficient  and  accurate  when  used  to  calculate  strong  shock  Inter¬ 
actions,  over  time  scales  on  the  order  of  one  to  several  shock  transit  times 
throughout  the  problem  region.  They  are  not  suited  to  the  analysis  of  very 
late-time  "non-hydrodynamlc"  phenomena,  after  shock  Interactions  are  essentially 
completed.  In  some  cases,  the  results  of  Lagranglan  calculations  may  be 
combined  with  simpler  methods  to  simulate  material  dynamics  at  late  times. 

The  successful  application  of  hydrodynamic  codes  hinges  on  a  detailed 
knowledge  of  the  constitutive  relations  of  all  materials  In  the  problem, 
throughout  the  expected  range  of  pressures  and  densities.  The  density,  det¬ 
onation  velocity,  energy  release,  and  pressure-volume  adiabat  of  the  explosive 
source  region  must  be  determined.  Similarly,  for  solid  media,  it  Is  necessary 
to  define  hydrostatic  compressibility  curves,  the  high-pressure  Hugonlot 
relation,  the  elastic  moduli  of  the  bulk  material,  and  shear  strength  at 
failure  as  a  function  of  confining  Dressure.  The  codes  are  therefore  most 
conveniently  applied  to  media  with  already-established  constitutive  relations. 
Extensive  materials  testing  may  be  required  to  derive  this  Information  for  a 
new  explosive  or  geologic  medium.  However,  the  cost  of  testing  may  be  more 
than  offset  by  the  reduced  need  for  expensive  site-calibration  experiments, 
and  by  the  decreased  likelihood  of  design  failure  or  unexpected  results. 

In  summary,  hydrodynamic  codes  are  well-suited  to  the  simulation  of 
shock  propagation  and  early-time  dynamics  in  multiple-layer  material  problems. 
The  effects  of  relevant  medium  parameters,  such  as  degree  of  saturation  and 
failure  strength  in  shear,  may  be  taken  Into  account.  Late-time  interactions 
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whlch  occur  after  the  phase  of  ballistic  material  ejection  cannot  be 
examined  with  this  approach.  Lagranglan  codes  are  thus  directly  applicable 
to  the  calculation  of  shock  transmission  and  material  ejection  from  under¬ 
water  cratering  events.  In  addition,  two-dimensional  codes  may  be  used  to 
determine  final  crater  dimensions  for  most  events  In  rock  media  (ballistic 
ejecta  deposition  not  greatly  modified  by  late-time  effects).  Crater  forma¬ 
tion  In  fine-grained  coheslonless  materials  Is  Influenced  by  washback,  slope 
failure,  and  settlement  effects  which  occur  during  or  after  deposition  of 
the  fallback.  For  this  reason,  the  final  crater  contours  cannot  be  accurately 
predicted,  although  an  Idealized  Initial  or  "ballistic"  fallback  distribution 
can  be  calculated.  It  Is  suggested  that  the  results  of  Lagranglan  and  ballistic 
ejecta  calculations  might  be  combined  with  a  late-time  slope  stability  and 
washback  analysis  In  order  to  determine  final  crater  dimensions.  This  approach 
has  not  yet  been  attempted,  but  might  be  justified  for  the  design  of  a  large- 
scale  excavation  project  In  a  weak  medium. 

The  hydrodynamic  codes  obviously  offer  considerable  benefits  within 
the  field  of  explosive  excavation  design.  This  Is  particularly  true  of  under¬ 
water  applications,  for  which  empirical  data  are  not  available  and  the  crater 
formation  process  Is  poorly  understood.  Perhaps  the  greatest  long-range 
advantage  to  be  derived  from  numerical  modeling  Is  the  Improved  understanding 
of  cratering  mechanisms  which  can  be  gained  by  studying  the  stress  and  velocity 
fields  In  the  excavated  material.  Detailed  Information  about  these  fields 
cannot  be  obtained  from  experimental  measurement  programs  at  any 
cost.[18] 

In  view  of  the  continuing  Interest  In  underwater  explosive  excavation 
and  the  need  for  more  quantitative  design  procedures,  It  was  decided  to 
conduct  a  calculatlonal  study  of  submerged  blasting.  The  initial  phase  of  this 
study  Involved  a  suite  of  one-dimensional  (spherically  symmetric)  calculations 
using  a  modified  version  of  the  SOC  code  which  contains  a  high-explosive  "burn" 
option.  Both  horizontal  and  vertical  calculations  were  performed  for  a  10-ton 
high-explosive  energy  source  emplaced  In  three  different  fully-saturated  media. 
The  vertical  calculations  Included  rock  "free-surface"  configurations  (no 
overlying  water  layer),  and  configurations  at  the  same  burial  depth  In  rock 
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wlth  overlying  water  layers  of  varying  thickness.  The  three  media  were 
selected  on  the  basis  of  the  following  criteria:  (1)  Availability  of  material 
properties  data  and  known  or  calculable  compressibility  curve;  (2)  Applica¬ 
bility  to  typical  high-explosive  excavation  projects;  (3)  Existence  of 
past  experimental  data;  and  (4)  Moduli  representative  of  a  range  of  geologic 
media  likely  to  be  encountered  In  field  work.  Detailed  discussion  of  the 
material  properties  Is  given  In  Chapter  V. 

The  primary  goals  of  the  calculation  series  were  fourfold:  (1)  To 
demonstrate  the  feasibility  of  hydrodynamic  code  application  to  multiple 
layer  cratering  configurations;  (2)  To  Investigate  the  properties  of  typical 
saturated  media  for  underwater  blasting  calculations;  (3)  To  achieve  a  basic 
understanding  of  the  submerged  cratering  process  (within  the  limitations  of 
the  one-dimensional  code);  and  (4)  To  develop  some  of  the  knowledge  required 
to  Implement  two-dimensional  design  calculations  for  future  engineering  projects. 
This  report  presents  a  summary  of  the  one-dimensional  study  and  related 
underwater  shock  Investigations  to  date.  Peak  pressures,  stresses,  and 
velocities  for  all  configurations  are  given.  Proper  Interpretation  of 
multiple  layer  dynamic  Interactions  requires  a  knowledge  of  the  waveforms  as 
well  as  peak  values.  For  this  reason,  a  data  abstract  of  the  calculated 
waveforms  Is  Included.  This  Information  should  also  be  very  useful  In 
developing  and  verifying  future  two-dimensional  calculations.  Since  the 
calculatlonal  suite  was  limited  to  three  typical  media,  the  analysis  Is  con¬ 
cerned  primarily  with  shock  Interactions,  early-time  dynamics,  and  water 
layer  effects  for  these  materials.  The  fundamental  dynamic  parameters  may  vary 
dramatically  with  medium  properties  and  with  layer  depths.  Thus,  caution 
should  be  exercised  In  extrapolating  from  or  generalizing  these  results; 
actual  crater  design  will  require  two-dimensional  calculations  for  a  specific 
medium  and  depth  configuration.  - 

II.  CRATER  FORMATION  MECHANISMS 

The  dynamics  of  explosive  crater  formation  In  "single-layer"  media  have 
previously  been  examined  In  some  detail Detonation  of  a  chemical 
explosive  charge  Is  completed  In  a  time  on  the  order  of  a  millisecond,  the  exact 
value  depending  on  charge  size  and  shape,  type  of  explosive,  and  other 
factors.  When  the  detonation  front  reaches  the  explosive  Interface,  a 
relatively  steep- fronted  shock  wave  Is  transmitted  Into  the  surrounding  medium. 


Inelastic  dissipative  effects  gradually  diffuse  the  sharp  peak  of  the 
outward-moving  wave.  Stresses  and  velocities  decay  or  attenuate  because 
of  dissipative  energy  loss  and  geometric  divergence  of  the  wave.  At  ranges 
very  close  to  a  strong  energy  source,  environmental  material  may  be  vaporized 
or  liquefied  by  the  shock.  These  effects  are  normally  absent  or  unimportant 
for  high-explosive  detonations  In  rock,  since  peak  pressures  near  the  explosive 
are  on  the  order  of  100  kbar  (energy  content  :0.1  Mbar-cm3/cm3) .  At  lower 
pressure  levels,  energy  deposition  continues  by  compresslonal  crushing, 
cracking,  or  deformation  of  the  material.  Finally,  as  the  peak  stresses 
near  the  shock  front  drop  to  very  low  levels,  major  dissipative  processes 
cease  to  occur.  The  transmitting  material  begins  to  respond  elastically, 
and  the  low-amplitude  shock  wave  propagates  as  a  sonic  pulse  In  the  local 
medium. 

This  Ideal  picture  of  shock  propagation  In  an  Infinite  homogeneous 
medium  Is  greatly  modified  by  the  Introduction  of  a  free  surface.  Boundary 
conditions  require  that  pressure  and  velocity  across  an  Interface  be  continuous. 
In  the  case  of  a  free  surface,  the  requirement  that  pressure1^"*  0  applies 
at  the  Interface.  Ignoring  density  changes  due  to  Irreversible  compaction, 
the  density  also  remains  unchanged  at  the  boundary.  Finally,  the  displace¬ 
ments  experienced  In  typical  shock  waves  are  sufficiently  small  so  that  changes 
In  gravitational  potential  are  negligible  compared  with  shock  compression 
energy.  Using  these  assumptions,  a  reflected  wave  must  occur  to  restore  all 
Initial  conditions  except  for  the  velocity  acquired  by  the  material  near  the 
Interface.  In  the  limiting  case  of  normal  Incidence,  the  boundary  condition 
P*0  Is  satisfied  by  a  reflected  wave  of  rarefaction  (negative  pressure) 
travelling  back  from  the  Interface.  The  rarefactlve  or  tensile  wave  Is  a 
mirror  Image  of  the  upward-travelling  compression  wave,  and  exactly  cancels 
the  wave  at  the  boundary.  The  rarefaction  front  propagates  Into  material 
which  has  a  local  particle  velocity  MVp"  (Induced  by  the  compresslonal  shock) 
toward  the  boundary.  The  rarefaction  must  change  this  velocity  by  an  amount 
"  -Vp"  relative  to  Its  direction  of  advance  (away  from  the  Interface),  leaving 
the  material  with  a  local  velocity  "2Vp"  toward  the  boundary.  Thus,  relative 
to  the  Initial  or  undisturbed  coordinate  system,  the  rarefaction  front  has 
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a  velocity  Vs-2vp  away  from  the  boundary  (downward),  where  "Vs"  Is  the 
upward  velocity  of  the  Incident  shock  wave;  the  boundary  Is  translated  up¬ 
ward  relative  to  the  Initial  coordinate  system  with  velocity  2Vp.  Velocity 
doubling  therefore  occurs  when  a  normally- Incident  wave  of  compression  Is 
reflected  at  a  free  surface.  A  similar  velocity  reinforcement  effect  may 
be  demonstrated  for  plane  waves  at  non-normal  Incidence.  The  velocity  vector 
of  the  Incident  wave  Is  normal  to  the  wavefront.  In  the  upward  direction 
of  propagation.  The  vector  of  the  reflected  tensile  wave  Is  normal  to  the 
wavefront,  opposite  the  downward  direction  of  propagation.  The  vertical 
components  of  the  Incident  and  reflected  waves  at  the  Interface  reinforce, 
doubling  the  vertical  velocity.  The  horizontal  components  cancel.  Thus, 
the  vertical  velocity  at  the  Interface  Is  doubled  for  oblique  Incidence, 
and  the  material  velocity  Is  strongly  redirected  toward  the  vertical.  This 
analysis  Is  oversimplified,  since  It  Ignores  Irreversible  (dissipative) 
effects  at  the  shock  front,  the  tendency  of  rarefaction  waves  to  "shock 
down"  (spread)  due  to  shock  velocity  variations  with  pressure,  and  failure 
phenomena  In  solid  media.  Likewise,  the  acoustic  paradox  for  grazing  Incidence 
at  a  boundary  cannot  be  taken  Into  account  (a  grazing  wave  produces  only 
tangential  velocity  changes,  and  yet  must  restore  the  material  to  Its  Initial 
pressure).  In  spite  of  these  problems,  the  simple  rarefactlve  scheme  provides 
an  accurate  description  of  Interactions  at  or  near  a  free  surface,  except  In 
the  case  of  very  high  pressures  or  oblique  angles  of  Incidence. 

Tensile  reflection  at  a  free  surface  gives  rise  to  one  of  the  two 
primary  cratering  mechanisms,  material  spallation.  Most  rock  media  cannot 
support  large  tensile  stresses  (l.e.,  the  Intact  rock  mass  will  break  apart 
easily  under  tension).  The  tensile  strength  of  In-sltu  rock  Is  further  reduced 
by  the  presence  of  joints  or  other  discontinuities.  Cratering  detonations 
are  burled  sufficiently  shallow  that  the  peak  shear  and  tensile 
stresses  near  the  surface  far  exceed  the  rock  strength.  Material  close  to 
the  free  surface  Is  thus  broken  by  the  reflected  tensile  wave.  If  not  already 
shattered  by  the  congressional  shock.  The  momentum  content  of  the  upward-moving 
compression  wave  and  the  downward-travelling  tensile  wave  Is  trapped  In  the 
cohesionless  fractured  material,  which  spalls  upward  at  approximately  twice  the 
vertical  peak  particle  velocity  Induced  by  the  incident  congressional  shock. 
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Thls  material  enters  a  state  of  freefall  and  describes  a  ballistic 
trajectory  (as  modified  by  air  resistance  or  collisions  between  neighboring 
particles).  As  the  head  of  the  tensile  wave  propagates  downward,  deeper 
layers  of  material  successively  undergo  spallation.  However,  the  peak  particle 
velocity  Is  not  doubled  for  this  material,  since  the  tensile  wave  arrives 
on  the  decaying  tall  of  the  compresslonal  shock  and  Is  attenuated  by  divergence 
and  Inelastic  propagation  effects.  The  deeper  spall  layers  enter  a  freefall 
phase  which  may  however  be  quickly  terminated  by  Impact  against  overlying 
layers  and  recompaction  by  underlying  material.  Spallation  may  be  thus  thought  of 
as  a  depth-and  range-dependent  phenomenon  whereby  layers  of  mound  material 
above  the  explosion  receive  a  vertical  velocity  reinforcement  from  the 
reflected  tensile  pulse.  If  a  detonation  Is  very  deeply  burled,  the  rock 
near  the  free  surface  may  remain  unbroken.  In  this  case,  the  Initially- 
doubled  velocity  peak  at  the  surface  may  be  damped  because  the  material  does 
not  break  apart  and  enter  freefall.  Spallation  parting  may  occur  at  greater 
depths,  below  the  surface,  and  thick  layers  of  material  may  be  launched  with 
somewhat  reduced  freefall  velocity. 

The  second  primary  cratering  mechanism,  gas  acceleration,  begins 
when  the  reflected  tensile  wave  reaches  the  explosive  cavity.  Prior  to 
this  time,  cavity  growth  Is  essentially  spherical  (slight  deviations  from 
spherical  free-fleld  growth  may  occur  due  to  depth-dependent  variations  In 
hydrostatic  pressure  around  the  cavity).  As  the  rarefaction  Impinges  on  the 
roof  of  the  cavity,  local  stresses  In  the  rock  are  relieved,  and  accelerated 
asymmetrical  growth  toward  the  free  surface  begins.  The  signal  from  the 
cavity  Is  transmitted  upward  at  less  than  the  preshot  acoustic  velocity  because 
of  the  rarefied  mound  material.  Indeed,  the  gas-acceleration  signal  may  never 
reach  the  surface  layers,  which  have  spalled  away  from  the  underlying  material. 
Nonetheless,  gas  acceleration  Is  often  an  Important  mechanism  for  transferring 
cavity  energy  to  the  lower  parts  of  the  mound.  The  effect  of  gas  acceleration 
is  to  progressively  recompact  high-velocity  material  adjacent  to  the  cavity 
with  overlying  material;  thus,  the  closest  adjacent  material  Is  decelerated 
and  mound  velocities  are  compacted  over  a  relatively  narrow  range.  Energy 


-12- 


1s  transmitted  from  the  Immediate  cavity  region  to  the  mound  as  a  whole. 

In  contrast  with  the  spall  phenomenon,  which  redirects  velocity  vectors 
toward  the  vertical,  gas  acceleration  forces  above  the  cavity  are  more 
nearly  radial  In  direction.  Gas  acceleration  Is  generally  most 
effective  for  rather  deeply-burled  events,  since  the  cavity  signal  must 
have  time  to  be  transmitted  throughout  a  significant  part  of  the  mound.  The 
Intense  spallation  wave  of  a  shallow  event  rapidly  launches  material  Into 
ballistic  trajectories,  and  mound  breakup  occurs  before  the  gas-acceleration 
"push"  can  become  effective. 

The  hydrodynamic  phases  of  the  explosive  cratering  process  are 
completed  when  the  coherent  material  mound  begins  to  disintegrate  and  the 
explosive  cavity  gases  vent  to  the  atmosphere.  By  or  before  this  time, 
significant  energy  transfer  and  velocity  changes  within  the  mound  have 
ceased  to  occur.  Most  of  the  material  ejected  from  the  crater  region  Is 
essentially  In  ballistic  freefall,  and  the  trajectories  may  be  calculated 
kinematically.  This  technique  Is  used  In  two-dimensional  hydrodynamic  calcu¬ 
lations,  which  are  carried  until  the  mound  dynamic  parameters  (such  as  kinetic 
energy  and  mound  mass)  approach  an  asymptotic  limit. The  hydrodynamic 
code  Is  then  linked  to  a  ballistic  throwout  code,  which  determines  the  final 
distribution  of  ejecta  and  the  crater  size.  This  procedure  may  create  slight 
errors  In  the  calculated  crater,  since  the  deep  material  below  the  ejecta 
mound  Is  still  somewhat  compressed  at  the  time  of  linking  and  will  later 
rebound  to  some  fraction  of  Its  original  volume.  Note,  however,  that  the 
all-important  dynamic  calculations  for  the  material  mound  above  the  shot 
are  correct  and  accurate.  The  most  crucial  aspects  of  the  cratering  process 
are  thus  adequately  described  by  the  hydrodynamic  code  model. 

The  above-discussed  crater  formation  process  undergoes  further 
modifications  If  gradual  material  changes  or  discontinuities  are  present 
In  the  cratering  medium.  Of  particular  Interest  Is  the  case  of  horizontally- 
stratified,  Impedance  mismatched  layers  above  the  shot  point.  Even  for 
shallowly  burled  explosions,  layers  of  successively  lower  density  and  sonic 

velocity  are  often  encountered  between  shot  point  and  the  surface.  The 
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surface  layer  Itself  may  consist  of  porous  weathered  material  of  low  sonic  velocity. 
Submerged  cratering  shots  present  a  similar  configuration,  since  the  dense 
rock  layer  Is  overlain  by  a  low  velocity  material  of  lower  Impedance.  The 
chief  difference  lies  In  the  fact  that  the  upper  layer  Is  a  fluid  with  no 
shear  or  tensile  strength.  In  either  case,  partial  tensile  reflection  of 
the  Initial  shock  will  occur  at  an  interface  between  higher-  and  lower- Impedance 
layers.  Vertical  particle  velocity  at  the  Interface  will  be  significantly 
Increased,  but  not  doubled.  Pressure  at  the  Interface  will  be  decreased  by 
the  tensile  wave  cutoff,  but  will  not  decline  to  zero.  A  weak  congressional 
shock  wave  will  be  transmitted  Into  the  low  Impedance  medium.  The  condition 
that  pressure  and  particle  velocity  must  be  continuous  across  the  Interface  will  be 
fulfilled;  thus,  the  approximate  peak  pressure  and  velocity  transmission  at 
the  Interface  can  be  calculated  from  the  acoustic  wave  equations.  Spallation 
may  develop  If  the  material  below  the  Interface  Is  sufficiently  fractured. 

However,  velocities  In  the  material  above  the  Interface  will  decline  as 
the  shock  Is  transmitted  away  from  the  Interface  (unless  the  overlying  layer 
Is  very  shallow,  and  It  too  spalls  almost  Immediately).  Spall  material  below 
the  Interface  may  then  Impact  upon  the  overlying  material,  ending  the  freefall 
phase.  Layer  Interactions  do  not  cease  at  this  point;  the  compress lonal  wave 
transmitted  Into  the  overlying  material  Is  subsequently  reflected  as  a  tension 
wave  at  the  next  Interface,  perhaps  creating  spallation  In  the  upper  layer. 

This  rarefaction  will  then  return  to  the  higher  Impedance  Interface,  where 
It  will  be  partially  reflected  again  (as  a  rarefaction).  Upper  layer  tensile 
pressures  will  be  Increased  by  this  re- reflect Ion,  but  the  upward  velocity 
component  will  be  decreased  (the  downward  velocity  component  of  the  upward- 
moving  tensile  re-reflection  partly  cancels  the  upward  component  of  the  downward- 
moving  tensile  wave).  Pressure  and  velocity  at  the  Interface  are  continuous, 
so  a  tensile  wave  with  an  upward  velocity  component  will  be  transmitted  Into 
the  lower  layer,  producing  an  Increase  In  the  vertical  velocities.  The  details 
of  this  second  shock  arrival  at  the  Interface  can  only  be  determined  by 
hydrodynamic  calculations,  since  the  pressures,  densities,  and  velocities  In 
both  layers  have  been  modified  by  earlier  dynamic  Interactions.  If  the  upper 
layer  has  not  yet  been  dispersed  by  spallation,  the  re-reflected  tensile  wave 
may  be  transmitted  upward  once  more,  suffering  yet  another  reflection  (as  a 
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downward- travel ling  congressional  wave).  The  compression  pulse  In  turn  may 
be  partially  reflected  (as  a  compression)  at  the  lower  Interface*  and  a 
compress lonal  wave  with  a  downward  velocity  component  will  be  transmitted 
Into  the  lower  layer.  Similar  reflections  will  be  generated  by  any  addi¬ 
tional  layers  In  the  configuration.  Even  for  simple  acoustic  waves*  the 
potential  complexity  of  multi-layer  Interactions  Is  quite  evident.  For 
finite-amplitude  shocks*  spall  dispersion  within  a  layer  may  decrease  or 
block  transmission  of  reflections  after  the  first  tensile  wave  (unless  the 
layer  Is  recompacted  when  a  reflection  arrives).  The  situation  Is  further 
clouded  by  the  time-dependence  of  material  velocities  and  of  compression  or 
dispersion  within  each  layer,  the  effects  of  material  strength, 
etc.  Finally*  the  acoustic  transmission  and  reflection  analysis  breaks  down 
completely  after  the  first  tensile  reflection  returns  to  the  cavity  and  the 
gas  acceleration  signal  propagates  Into  the  medium.  The  gas  acceleration 
pulse  Itself  will  be  subject  to  Interactions  at  various  Interfaces  as  It 
recompacts  the  successive  layers*  and  subsidiary  gas  acceleration  pulses  may 
be  generated  as  later  tensile  reflections  reach  the  cavity.  Hydrodynamic 
calculations  are  required  to  treat  the  general  problem  of  finite-amplitude 
multi pie- layer  shock  propagation  In  nonlinear  media  (material  moduli  change 
with  pressure). 


III.  CALCULATION/*!.  SUITE 

Due  to  the  complex  Interrelationships  between  multilayer  shock  Inter¬ 
actions  and  material  dynamics,  no  generalized  model  can  be  expected  to  simulate 
the  cratering  process  for  all  media  and  layer  depths.  Indeed*  all 
dynamic  Interactions  after  the  first  Interface  reflection  will  be  affected 
by  medium  properties  and  shock  strength.  For  this  reason*  the  sample  materials 
and  layer  depths  were  selected  to  conform  closely  with  realistic  high-explosive 
cratering  designs.  Eighteen  one-dimensional  (spherlcally-symmetrlc)  SOC  calcu¬ 
lations  were  performed  for  high-explosive  detonations  In  three  saturated  media: 
coral*  basalt,  and  granite.  Further  calculations*  suggested  by  the  results  of 
this  first  series,  were  later  added  for  comparison. 


mmrrw'1'  * 


m 


-15- 


The  energy  source  In  ell  calculations  was  a  sphere  of  aluminized 
anmonlum  nitrate  slurry  explosive  (AANS)  weighing  approximately  10  tons. 

The  spherical  source  region  had  a  radius  of  1.13  m  (3.71  ft)  and  was 
divided  Into  30  calculatlonal  zones  of  equal  thickness.  This  source  was 
"detonated"  In  both  horizontal  and  vertical  configurations  for  each  of  the 
three  media.  Horizontal  calculations  were  first  performed  to  simulate  free- 
fleld  shock  propagation  In  a  semi-infinite  rock  medium.  The  problem  size 
(100  m  «  328  ft)  was  much  larger  than  the  region  of  Interest  for  close- 
range  effects,  and  the  calculations  were  run  until  late  time  (75  msec).  The 
horizontal  calculations  assumed  a  constant  overburden  pressure  of  2.02  bars 
at  all  ranges  (l.e.,  constant-depth  free-fleld  environment).  This  pressure 
Is  equivalent  to  about  36  ft  of  saturated  coral  overburden,  or  a  somewhat 
lesser  depth  of  dense  rock.  This  small  overburden  pressure  had  little  effect 
on  close-range  shock  calculations.  The  horizontal  free-fleld  calculations 
are  useful  for  defining  shock  propagation  In  the  absence  of  a  free  surface, 
but  resolution  Is  somewhat  degraded  and  the  front  Is  smeared  by  the  coarse 
zone  size  required  In  a  large  problem.  Vertical  calculations  were  next 
conducted,  with  the  explosive  source  emplaced  below  a  rock  surface  with  or  with¬ 
out  water  overburden  (the  upper  boundary  In  each  case  being  a  free  surface). 
The  depth  of  burial  beneath  the  rock  Interface  "Dr"  was  36  ft  or  30  ft 
(Dr  *  range  from  the  center  of  the  explosive  to  the  rock  Interface).  The 
36-ft  depth  is  near  optimum  depth  for  crater  size  from  10-ton  charges,  and 
also  simulates  the  Project  Tugboat  underwater  cratering  events  In 
coral.  Water  layer  depths  "Dw"  above  the  rock-water  Interface  were  0  ft  (no 
water,  rock  Interface  Is  a  free  surface),  5  ft,  10  ft,  and  36  ft.  Due  to 
time  limitations,  every  possible  combination  of  rock  and  water  depth  could 
not  be  calculated  for  each  of  the  three  media.  Therefore,  the  D  =  36  ft 
configurations  were  emphasized.  Since  the  size  of  the  vertical  problems  was 
quite  small,  very  fine  constant  zone  size  (about  0.04  m  or  0.13  ft)  was  used 
In  the  rock  layer.  This  zoning  achieved  excellent  resolution  of  the  shock 
front.  Vertical  problems  were  run  for  about  20-28  msec  to  allow  ample  time 
for  the  Initial  compressive  wave  transmission  and  return  of  the  tensile 
reflection.  Strong  dynamic  interactions  In  the  vertical  direction  are 
essentially  completed  by  this  time.  Figure  1  schematically  depicts  the 
horizontal  and  vertical  problem  configurations  (not  to  scale). 
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Horlzontal  calculations  (run  to  75  msec  ;  overburden  =2.02  bars) 


Figure  1.  Horizontal  and  vertical  problem  configurations  for 
the  one -dimensional  spherical  (SOC)  calculations. 
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Abbreviations  have  been  devised  for  convenience  In  discussing  the 
various  calculations.  Medium  names  are  denoted  as  follows:  COR  (coral), 

BAS  (basalt),  and  GRN  (granite).  The  medium  name  Is  followed  by  the  configu¬ 
ration,  either  HRZ  (horizontal),  or  by  the  vertical  calculation  rock  depth/ 
water  depth.  In  feet  (l.e.,  30/10  represents  Dr  =  30  ft,  Dw  ■  10  ft).  For 
example,  the  coral  horizontal  calculation  Is  designated  CORHRZ.  The  basalt 
calculation  at  a  depth  of  36  ft  with  a  36  ft  overlying  water  layer  Is 
BAS36/36.  One  additional  horizontal  calculation  In  coral,  using  a  modified 
constitutive  relation,  Is  denoted  "PS10HRZ'1.  The  nomenclature  scheme  Is 
summarized  In  Table  1. 

Previous  studies  with  SOC  and  other  hydrodynamic  codes  have  demonstrated 
that  generation  of  physically  meaningless  numerical  noise  often  occurs  at 
Impedance-mismatched  Interfaces  within  a  problem.  Such  noise  usually  takes 
the  form  of  a  "ringing"  or  echo  effect  originating  at  the  location  of  a  sudden 
change  In  zone  size  (for  single-material  problems)  or  at  the  Interface  between 
two  different  materials.  Noise  generation  Is  most  severe  for  sharply- rising 
Incident  shock  fronts,  and  for  large  Impedance  mismatches.  Difficulties  of 
this  sort  may  be  reduced  by  avoiding  sudden  changes  In  zone  size  (within  a 
material),  and  by  Impedance  matching  of  adjacent  zones  across  an  Interface 
(zone  size  proportional  to  the  acoustic  Impedance  of  the  medium;  coarser  zones 
In  the  "soft"  low  Impedance  medium). 

The  latter  technique  has  been  used  at  the  rock-water  discontinuity 
In  these  problems.  Matching  across  a  very-high  to  very-low  Impedance  Interface 
may  degrade  resolution,  slrce  much  larger  zones  are  required  for  the  overlying 
material.  The  shock  front  may  undergo  dramatic  spreading  as  It  propagates 
Into  the  low- Impedance  layer.  For  similar  reasons,  reflection  phenomena  at 
the  Interface  may  be  Inadequately  resolved;  the  reflected  shock  front  rapidly 
transits  several  zones  of  the  dense  material  before  the  overlying  low-impedance 
(low- velocity)  layers  become  fully  active.  Significant  errors  occur  when  the 
distance  travelled  by  the  Incident  shock  front  during  Its  rise  time  and  the 
depth-scale  of  reflected  shock  Interactions  near  the  Interface  are  smaller  than 


Emplacement  Configurations  for  Which  Calculations  were  Performed;  Energy  Source  is 
a  Ten-ton  Spherical  Charge  of  Aluminized  Ammonium  Nitrate  Slurry  (15.9  tons  equiv. 
energy  yield). 
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the  zone  size  In  the  overlying  material.  Resolution  at  the  Interface  and 
In  the  low  Impedance  layer  may  be  Improved  by  using  very  fine  zones  In 
both  materials,  at  the  cost  of  Incr-e^d  calculation  time. 

During  this  study,  It  was  conclusively  shown  that  the  water- transmitted 
shock  front  and  the  Interactions  at  the  Interface  could  be  well  resolved 
for  the  smoothly-rising  shock  waves  characteristic  of  coral  or  basalt.  In 
spite  of  the  coarse  water  zone  sizes  required  for  Impedance  matching.  How¬ 
ever,  the  very  short  rise  times  and  the  sharp  fronts  In  granite  gave  rise  to 
some  Inconsistencies.  The  shock  wave  underwent  abrupt  spreading,  with  a  decrease 
In  particle  velocity  upon  entering  the  water  layer.  Early-time  reflection 
phenomena  In  granite  zones  adjacent  to  the  Interface  were  likewise  Improperly 
calculated.  These  difficulties  are  accentuated  by  the  fact  that  the  Impedance 
discontinuity  Is  greatest  for  the  dense  granite  medium,  requiring  coarse 
water  zones.  Accurate  Impedance  matching  Is  also  most  necessary  for  sharply- 
rising  shock  fronts.  Finally,  small  granite  zone  sizes  cannot  be  used  due 
to  the  long  calculation  time  required  for  a  "stiff"  high-velocity  material. 

In  order  to  define  the  errors  Incurred,  additional  granite  calculations  were 
performed  with  finer  overlying  water  zones  (not  Impedance  matched).  These  calcu¬ 
lations  included  a  set  with  water  zones  half  the  size  of  thfe  "impedance  matched" 
zoning.  It  was  found  that  the  resolution  factor  grossly  affected  peak  shock 
transmission  In  the  water  layer,  but  had  little  Influence  on  long-term  dynamics 
within  the  granite  layer. 

Tables  2a  and  2b  summarize  the  Interface  locations,  rock  zone  sizes, 
and  water  zone  sizes  featured  In  the  set  of  vertical  calculations.  Interface 
locations  (measured  from  the  center  of  the  problem)  and  the  zone  sizes  are 
given  In  both  meters  and  feet.  Note  that  the  Interface  positions  and  rock 
zoning  are  Identical  for  corresponding  problems  in  all  three  media.  However, 
the  water  zone  sizes  (Table  2b)  vary  because  of  the  impedance  matching  require¬ 
ment.  In  the  case  of  granite,  two  different  water  zone  sizes  are  listed.  The 
first  refers  to  the  impedance-matched  calculations  ("coarse  water  zones"), 
while  the  second  corresponds  to  half-size  water  zoning  ("fine  water  zones"). 


TABLE  2(a).  Interface  Locations  and  Approximate  Rock  Zone  Sizes  Used  In  the 
Vertical  Calculations. 


Explosive-Rock 

Interface 

Calculation*  (m)  (ft) 


00 
05 
30/10 
30/36 


Ranae  from  Center 


Rock-Water  Water 

Interface  Surface 

(m)  (ft)  (m)  (ft) 


1.13  3.71 


None 

None 

10.67 

35 

'9.146 

30 

12.195 

40 

20.12 

66 

None 

None 

►10.976 

36 

12.50 

41 

14.02 

46 

21.95 

72 

Rock  Zone 
Size* 

(m)  (ft) 


0.04  0.131 


0.04  0.131 


Identical  Interface  locations  and  rock  zone  sizes  are  used  for  the  corresponding 
calculations  In  all  three  media 


TABLE  2(b).  Approximate  Water  Zone  Sizes  Used  In 
Vertical  Calculations 


Water  Zone  Size 


Coral  0.10 

Basalt  0.254 

Granite  (coarse  water  0.392 

zones) 

Granite  (fine  water 
z 
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Addltlonal  test  calculations  with  somewhat  modified  Interface  locations, 
depths,  or  zone  sizes  were  subsequently  added  to  the  study.  These  calcula¬ 
tions  will  be  discussed  and  compared  with  the  "standard"  configurations 
(Table  2)  In  a  later  chapter. 

IV.  ONE-DIMENSIONAL  VS.  TWO-DIMENSIONAL  CALCULATIONS 

The  spherical ly-dlvergent  geometry  Imposed  In  SOC  calculations  presents 
certain  difficulties  of  Interpretation  when  applied  to  cratering  dynamics. 

Crater  formation  Is  an  Inherently  two-dimensional  process;  It  Involves,  In 
the  simplest  case,  propagation  of  a  spherical ly-dlvergent  shock  front, 
reflection  at  a  planar  Interface,  followed  by  hydrodynamic  acceleration  and 
ballistic  ejection  of  the  crushed  material.  Oblique  reflection  of  the  shock 
wave  plays  an  Important  role  In  fracturing  material  near  the  edge  of  the 
crater  and  In  lateral  ejection.  Gas  acceleration  effects  In  the  lateral 
direction  may  also  contribute  to  crater  formation.  Dynamics  of  this  laterally 
displaced  material  are  not  directly  related  to  the  dynamics  vertically  above 
the  charge;  Instead,  they  are  dependent  on  shock  attenuation,  medium  strength, 
emplacement  geometry,  and  other  factors.  For  these  reasons,  the  volume  of 
material  ejected  and  crater  dimensions  cannot  be  accurately  predicted  from 
one-dimensional  results.  Likewise,  the  "optimum"  depth  of  burial  for  a 
cratering  charge  depends  upon  medium  properties,  and  Is  not  simply  related 
to  the  vertical  spall  velocity,  residual  velocity,  kinetic  energy,  or  other 
one-dimensional  dynamic  parameters.  With  these  cautionary  notes  In  mind.  It 
Is  permissible  to  observe  that  there  Is,  nonetheless,  a  considerable  degree 
of  correlation  between  crater  formation  and  close-range  vertical  dynamics. 

This  correlation  Is  expressed  by  the  relationship  between  the  vertical  spall 
velocity  In  the  near-surface  layers  and  the  required  depths  of  burial  for 
"optimum"  crater  formation  or  for  "mounding"  (no  apparent  crater).  Before  the 
advent  of  two-dimensional  codes,  the  vertical  velocity  was  successfully  used 
as  a  rough  criterion  for  burial  depth  In  a  variety  of  media.  Similarly,  the 
cube-root  scaling  laws  which  describe  shock  front  parameters  (velocity,  stress, 
spall  velocity)  also  apply  approximately  to  burial  depth  and  crater  dimensions 
In  competent  media  at  "small"  length  scales  (l.e.,  typical  high-explosive  events) 
and  over  small  ranges  of  yield.  These  simple  laws  are  of  course  not  valid 


over  great  yield  ranges  and  at  large  yields  where  gravitational  scaling  effects, 
air  viscosity,  and  other  factors  become  Important. ^ *22,233  H0wevert  gjven 

the  conditions  of  small  yield  range,  small  variation  In  the  size  scale  or 
depth  of  the  experiment,  and  similar  media,  one  expects  that  the  other  "non- 
hydrodynamlc"  scaling  effects  will  be  unimportant  and  a  good  correlation  will 
exist  between  vertical  dynamics  and  efficiency  of  material  ejection. 

No  detailed  study  has  been  conducted  to  relate  one-dimensional  (vertical) 
and  two-dimensional  dynamic  parameters  for  burled  cratering  events.  Unfortunately, 
two-dimensional  calculations  are  not  yet  available  for  media  Involved  In  this 
Investigation.  The  nature  of  the  difference  between  one-  and  two-dimensional 
calculations  Is  Illustrated  In  Figure  2a.  This  figure  compares  the  kinetic  energy 
content  In  coral  (for  one-dimensional  calculations)  with  the  kinetic  energy  content 
In  Bearpaw  clay  shale  (for  two-dimensional  calculations). [19'24]  The  coral  results 
(solid  curves)  are  shown  for  a  horizontal  (free-fleld)  configuration  and  a 
vertical  configuration  at  a  burial  depth  of  36  ft  (coral  36/00,  no  water  layer). 

The  two-dimensional  clay  shale  kinetic  energy  curves  (dashed)  are  plotted  for  the 
total_  kinetic  energy  In  the  clay  shale  and  the  kinetic  energy  content  In  the  ejecta 
mound  (material  thrown  above  the  Initial  ground  surface  level).  Clay  shale  Is  a 
weak  fully-saturated  material  with  a  compressibility  curve  quite  similar  to  coral. 
However,  certain  differences  should  be  mentioned:  (1)  Clay  shale  has  a  density 
and  Impedance  about  17%  greater  than  that  of  coral,  while  sonic  velocities  for 
the  two  media  are  about  the  same.  (2)  The  clay  shale  calculations  use  a  different 
explosive,  nltromethane;  the  total  energy  release  for  the  clay  shale  calculations' 

Is  1.23  times  as  great  as  the  energy  release  of  the  AANS  charge  used  In  the  coral 
calculations.  Note  that  the  clay  shale  total  kinetic  energy  curve  lies  above  the 
coral  curve  by  a  factor  of  approximately  1.23  at  early  time.  (3)  The  clay  shale 
calculations  were  at  a  depth  of  burial  of  41  ft,  or  1.14  times  as  deep  as  CORAL 
36/00.  Therefore,  the  shock  arrives  at  the  free  surface  slightly  later.  (4)  Clay 
shale  Is  a  very  weak,  semi-fluid  material  with  a  strength  on  the  order  of  the  over¬ 
burden  pressure;  therefore,  cavity  expansion  would  continue  almost  unhindered  at  late 
time,  even  for  a  spherical  free-fleld  calculation.  Kinetic  energy  continues  to  enter 
the  medium  by  a  process  similar  to  fluid-flow  expansion.  The  ratio  of  kinetic  energy 
to  Internal  energy  In  the  shale  medium  does  not  change  dramatically  with  time, 
ihese  contentions  are  borne  out  by  free-fleld  one-dimensional  calculations 
u*ot  shown),  which  reveal  a  total  kinetic  energy  content  similar  to  Figure  2a, 
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Figure  2a.  Cor.  arison  of  kinetic  energy  in  coral  (one- 
aimensional  vertical  calculations)  and  in 
clay  shale  (two-di::iensional  calculations). 


decreasing  very  slightly  at  late  times.  The  stronger  coral  material.  In 
contrast,  halts  cavity  expansion  much  sooner  (see  coral-horizontal),  the  kinetic 
energy  content  begins  to  decrease,  and  the  Internal  energy  of  the  medium  Increases 
correspondingly.  For  this  reason,  the  total  kinetic  energy  curve  for  clay  shale 
may  not  be  legitimately  compared  with  coral  after  2  or  3  msec.  A  two-dimensional 
coral  calculation  would  show  a  large  relative  decrease  In  total  kinetic  energy 
at  late  times  due  to  severe  damping  of  expansion  beneath  the  lower  half  of  the 
cavity.  However,  the  mound  kinetic  energy  for  clay  shale  may  be  compared  with 
coral  36/00  at  late  times  (after  the  coral  is  almost  completely  fractured  and 
has  entered  freefall  spallation).  The  clay  shale  curve  has  been  multiplied 
by  a  factor  of  1.63  to  facilitate  comparison  (upper  dashed  curve).*  Shock 
arrival  at  the  free  surface  (s  4.5  to  5  msec)  is  Indicated  In  the  figure;  spalla¬ 
tion  begins  at  this  time.  Note  that  the  two  curves  show  a  similar  percentage 
Increase  at  late  times,  but  the  two-dimensional  (clay  shale)  kinetic  energy 
rises  gradually  to  its  maximum  value  between  5  and  20  msec.  The  coral  curve 
ascends  abruptly,  and  spallation  Is  essentially  complete  throughout  the  vertical 
problem  by  10  msec.  It  Is  evident  that  spall  launch  takes  longer  to  propagate 
to  long  slant  ranges,  throughout  the  two-dimensional  mound.  However,  the  ulti¬ 
mate  Increase  In  kinetic  energy  during  the  spallation  phase  Is  encouragingly 
similar  to  the  one-dimensional  result.  Comparing  the  true  clay  shale  mound 
energy  (lower  dashed  curve)  with  coral  36/00,  It  Is  seen  that  the  coral  kinetic 
energy  lies  above  the  clay  shale  by  a  factor  of  1.56  to  1.68  after  spallation 
and  most  dynamic  Interactions  are  complete  (l.e.,  at  times  later  than  about 
10  msec  for  the  vertical  coral  36/00,  or  20  msec  for  the  two-dimensional 
clay  shale  configuration).  Correcting  for  the  input  energy  differences 
between  the  two  problems,  this  corresponds  to  a  factor  of  1.92  to  2.07  between 
the  two-dimensional  and  one-dimensional  (vertical)  results,  roughly  as 
expected  for  energy  division  between  upper  and  lower  hemispheres.  The 
very  accurate  "factor  of  two"  difference  observed  here  occurs  fortuitously 
because  of  the  very  low  strength  properties  and  high  shock  propagation  effi¬ 
ciency  In  clay  shale,  and  also  because  of  other  slight  differences  between  the 
problems.  A  two-dimensional  coral  calculation  would  ultimately  place  some¬ 
what  less  than  half  of  the  corresponding  one-dimensional  vertical  kinetic 

*  This  multiplicative  correction  factor  is  used  because  only  about  half  of  the 
total  energy  In  the  two-dimensional  calculation  propagates  into  the  upper 
hemisphere  (region  of  mound  launch).  However,  the  total  energy  is  1.23  times 
greater.  Thus,  a  rough  compensation  factor  relating  the  two-dimensional 
mound  launch  energy  to  the  one-dimensional  vertical  calculation  KE  is  given 
by:  2x(l/1.23)  «=  1.63. 
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energy  Into  mound  kinetic  energy.  Nonetheless,  behavioral  similarities 
shown  here  Indicate  a  strong  analogy  between  the  dynamics  of  the  one-dimensional 
and  two-dimensional  calculations.  A  more  detailed  Investigation  Is  needed 
to  fully  verify  the  relationship. 

Other  more  subtle  differences  also  complicate  the  Interpretation 
of  one-dimensional  calculations.  Chief  among  these  Is  the  less-than- 
satlsfactory  handling  of  reflected  and  refracted  wave  propagation.  One- 
dimensional  spherical ly-dlvergent  calculations  produce  spherical  convergence 
of  reflected  waves  returning  toward  the  center  of  the  problem.  This  effect 
Is  not  serious  over  time  periods  and  distance  Intervals  which  are  small  compared 
with  the  propagation  time  and  range  from  the  problem  center.  The  shock  Inter¬ 
actions  near  an  Interface  and  far  removed  from  the  center  are  almost  planar. 
Thus,  peak  surface  velocities  and  early  time  spall  dynamics  close  to  an 
Interface  may  be  accurately  calculated.  The  spherical  convergence  effect 
Is  In  one  sense  self-cancelling,  even  over  larger  regions  of  the  problem. 

The  spherical ly-dlvergent  tensile  reflection  for  a  two-dimensional  cratering 
event  spreads  throughout  a  broad  region  of  the  mound  (rather  than  converging 
on  the  center),  but  still  contributes  to  spallation  and  kinetic  energy  of 
the  mound  as  a  whole.  Spall  dynamics  arc  thus  fundamentally  similar  for  one- 
dimensional  and  two-dimensional  cases.  However,  significant  differences  In 
the  material  velocities  and  dynamic  Interactions  at  a  given  radius  are  likely 
to  occur  as  the  one-dimensional  wave  propagates  back  to  the  cavity  and  Is 
again  reflected  outward.  The  worst  "errors"  (relative  to  a  corresponding 
two-dimensional  problem)  are  expected  to  originate  within  the  extreme  Inner 
regions.  For  example,  the  late-time  cavity  gas  acceleration  effects  In  the 
vertical  direction  might  be  overestimated  for  a  one-dimensional  problem,  since 
the  returning  tensile  wave  which  reaches  the  cavity  Interface  Is  convergent 
rather  than  spherically  divergent.  The  effects  will  be  most  severe  and  will 
become  apparent  soonest  at  ranges  close  to  the  cavity.  However,  the  Inner 
regions  of  the  problem  account  for  only  a  small  portion  of  the  total  mass. 
Changes  In  the  overall  material  dynamics  may  be  very  small  and  may  occur  at 
late  times.  Further,  If  the  overlying  material  Is  rarefied  by  spallation, 
reloading  signals  from  the  cavity  propagate  outward  quite  slowly  and  may  not 
reach  the  bulk  of  material  nearer  the  free  surface  within  a  time  Interval 
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relevant  to  crater  formation.  The  problems  considered  here  are  relatively 
shallow  cratering  events  In  fully-saturated  media.  Previous  Investlga- 
t1onsL  *  J  and  results  of  this  study  demonstrate  that  shock  propagation 
In  saturated  media  Is  a  very  efficient  process,  and  that  spallation  Is 
the  predominant  mechanism  of  crater  formation.  The  upper  regions  of 
material  remain  In  freefall  long  after  spallation  begins.  Recompaction  and 
cavity  echo  signals  reach  the  outer  regions  only  at  very  late  times  (t  »  two 
acoustic  transit  times)  and  are  extremely  weak.  Gas  acceleration  and  multiple 
reflections  from  the  central  region  exert  little  Influence  on  the  overall 
dynamics  for  problems  of  this  class. 

The  analogy  between  kinetic  energy  In  one-  and  two-dimensional  problems 
has  already  been  examined  (Figure  2a).  The  one-dimensional  coral  36/00  calcu¬ 
lation  shown  In  Figure  2a  reveals  a  sudden  rise  In  kinetic  energy  from  5  to  10 
msec  (spallation  phase),  followed  by  a  very  gradual  Increase  at  later  times. 
Cavity  gas  acceleration  and  late-time  recompaction  effects  are  relatively 
unimportant  within  the  scope  of  overall  dynamics  to  20  msec.  The  mound  kinetic 
energy  for  the  two-dimensional  clay  shale  calculation  presents  similar  behavior, 
although  the  spall  wave  takes  longer  to  spread  throughout  the  mound  region. 

Gas  acceleration  effects  have  likewise  been  shown  to  be  unimportant  In  the 
two-dimensional  configuration.^  Kinetic  energy  Is  a  velocity-sensitive 
parameter  which  relates  closely  with  material  spallation  and  cratering  dynamics. 
The  correlation  between  one-  and  two-dimensional  results  thus  demonstrates 
that  spherical  convergence  does  not  adversely  affect  problem  dynamics  over 
times  of  Interest  (*  20  msec).  However,  the  one-dimensional  kinetic  energy 
Is  weighted  by  "mass  times  velocity  squared,"  de-emphasizing  the  low-mass 
Inner  zones.  Small  velocity  changes  near  the  center  of  the  problem  are  not 
readily  detected.  In  an  effort  to  analyze  detailed  dynamic  effects,  the 
spall  velocities  and  waveforms  above  shot  point  were  examined.  Due  to  data 
availability,  a  different  set  of  problems  was  used  for  this  examination.  A 
vertical  one-dimensional  coral  problem  at  a  depth  of  41.31  ft  Is  compared 
with  a  two-dimensional  clay  shale  event  at  a  depth  of  41  ft.  The  clay  shale 
problem  In  this  case  consists  of  a  hypothetical  nuclear  source  and  has  a  total 
kinetic  energy  of  only  0.6  times  the  coral  problem  energy  at  late  times.  Thus, 
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the  shock  stresses  and  velocities  far  from  the  source  are  lower  than  for 
coral.  However,  the  nuclear  shock  wave  arrival  times  are  slightly  earlier 
due  to  the  high  Intensity  of  the  close-range  nuclear  shock  when  It  is  decoupled 
from  the  source.  Figure  2b  shows  the  velocity  waveforms  at  a  range  of 
9  m  (29.5  ft)  above  the  problem  center.  Note  the  earlier  arrival  and 
slightly  lower  amplitude  of  the  clay  shale  wave.  Shock  Interactions  at 
spallation  are  similar  for  the  one-  and  two-dimensional  problems.  In  both 
cases,  the  spall  wave  arrives  about  3.2  msec  after  the  Initial  peak  and  causes 
a  similar  Increase  In  velocity.  The  freefall  phase  persists  until  very  late 
time  (>  20  msec),  and  no  gas  acceleration  arrival  occurs.  Approximately  62X 
of  the  total  mass  Is  located  at  or  beyond  9  m  for  the  one-dimensional 
problem.  The  velocity  waveforms  at  other  ranges  have  also  been  compared; 

Figure  2c  graphs  the  vertical  spall  velocities  as  a  function  of  range  from 
the  problem  center.  Coral  spall  velocities  are  slightly  higher  near  the  free 
surface  due  to  the  greater  amplitude  of  the  coral  shock  wave  at  long  ranges. 
However,  at  locations  very  close  to  the  source  (4  to  6  m),  the  spall 
velocities  are  equal.  This  effect  Is  caused  by  the  slower  decline  in  particle 
velocity  at  late  times  for  clay  shale  (fluid  flow  effect  mentioned  earlier). 

The  tall  of  the  upward- travel ling  compression  wave  declines  more  slowly  for 
clay  shale  than  for  coral,  thus  maintaining  a  slightly  higher  velocity  when 
the  spallation  wave  arrives.  It  Is  quite  evident  that  the  expected  Increase 
In  coral  spall  velocity  near  the  problem  center  (due  to  spherical  convergence) 

Is  less  Important  than  this  minor  medium-dependent  effect.  This  result  Implies 
that  convergence  has  relatively  little  influence  on  spallation  throughout  most 
of  the  one-dimensional  problem  region.  Late-time  gas  acceleration  effects 
close  to  the  source  were  examined;  It  was  found  that  the  coral  velocity  after 
gas  acceleration  was  approximately  58  m/sec  at  a  range  of  6  m;  the 
corresponding  velocity  for  clay  shale  was  5t>  m/sec  at  the  same  range  (the 
spall  velocity  before  gas  acceleration  Is  50  m/sec  In  both  cases.  Fig.  2c).  This 
slight  difference  may  Indicate  a  convergence  effect  near  the  center  of  the 
one-dimensional  coral  calculation,  but  the  velocity  change  Involved  Is  negligible 
and  could  be  attributed  to  differences  In  the  media  and  energy  sources.  It 
appears  that  spherical  convergence  has  little  effect  on  one-dimensional  spall - 
dominated  cratering  calculations.  This  conclusion  might  not  apply  to  very 


Vertical  spall  velocities  above  shot  point:  two-dimensional  clay  shale  calculation 
(depth  =  4l  ft)  and  one -dimensional  (vertical)  coral  calculation  (depth  =  4l.31  ft). 
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deeply  burled  events,  for  which  the  gas  acceleration  pulse  recompacts  a 
large  portion  of  the  mound  and  contributes  significantly  to  late-time 
material  acceleration. 

The  spherical  geometry  used  In  one-dimensional  configurations  also 
falls  to  take  account  of  refraction  at  medium  Interfaces.  The  problem 
of  reflection  and  transmission  of  a  low  amplitude  (acoustic)  wave  at  a 
discontinuous  boundary  and  subsequent  propagation  of  the  reflected  and  trans¬ 
mitted  waves  Is  analogous  to  the  problem  of  light  reflection- refraction  In 
geometric  optics.  General  equations  may  be  written  for  simple  hydrodynamic 
compression  waves  (In  a  fluld)^^  or  for  compression  plus  shear  waves. ^27^ 

The  acoustic  transmission  and  reflection  conditions  and  geometric  propagation 
laws  apply  with  a  high  degree  of  accuracy  to  weak  shock  waves  In  fluids,  and 
are  known  to  be  valid  even  for  very  strong  shocks  as  long  as  the  angle  of 
Incidence  does  not  deviate  too  far  from  the  normal. Failure  and  hysteretlc 
effects  complicate  strong  shock  transmission  In  solid  media,  but  the  acoustic 
and  geometric  laws  are  usually  applicable  to  peak  shock  parameters  across  an 
Interface.  As  discussed  earlier,  acoustic  theory  for  a  high-  to  low-impedance 
Interface  requires  that  an  Incident  compresslonal  shock  be  partially  reflected 
as  a  tensile  wave.  The  pressure  at  the  boundary  Is  reduced,  and  the  vertical 
velocity  component  Is  Increased.  The  transmission  and  reflection  coefficients 
depend  only  on  the  characteristic  Impedances  of  the  two  media  and  the  angle 
of  Incidence,  approaching  limiting  values  at  normal  Incidence.  The  Impedances 
may  be  sonic  values  "pCson1c"  for  an  acoustic  wave  (Csonic  =  sonic  velocity), 
or  finite-amplitude  values  "pCshocj("  for  a  discontinuous  shock  front 
^shock  *  s^oc,(  fr°nt  velocity).  The  "acoustic"  equations  apply  In  either 
case,  assuming  that  the  reversible  Ranklne-Hugonlot  jump  conditions  are  valid 
at  the  front.  The  finite  difference  equations  used  In  hydrodynamic  codes  are 
not  completely  compatible  with  the  discontinuous  "jump"  conditions  at  a 
boundary.  Nonetheless,  most  codes  will  accurately  treat  transmission  at  an 
Interface  if  the  problem  is  carefully  formulated.  The  jump  conditions  are 
sometimes  used  to  check  the  accuracy  of  a  solution. 
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Shock  propagation  within  the  second  medium  presents  another  unrelated 
problem  for  the  calculation  scheme.  Refraction  of  the  wavefront  takes  place 
at  an  Interface.  The  angle  of  refraction  Is  dependent  on  the  sonic  velocities 
(or  front  velocities)  In  the  two  media  and  the  angle  of  incidence,  according 
to  the  familiar  law  of  geometric  optics.  Refraction,  In  turn,  modifies  the 
radius  of  curvature  of  the  wavefront  In  two  respects.  The  first  of  these  Is 
spherical  aberration,  caused  by  the  shift  In  the  wavefront's  direction  of 
propagation  relative  to  the  axis  of  symmetry  (vertical  axis  perpendicular  to 
the  Inte-face).  Upon  transmission  at  a  high-velocity  to  low-velocity  Inter¬ 
face,  the  ray  Is  refracted  toward  the  normal  direction.  Increasing  the  apparent 
distance  to  the  axis  of  symmetry  and  the  radius  of  curvature.  As  the  radius 
of  curvature  Increases,  the  shock  attenuation  rate  In  the  second  layer  decreases 
(relative  to  an  unrefracted  spherlcally-dlvergent  wavefront).  The  second 
effect  Is  due  to  the  change  In  the  radius  of  curvature  of  the  wavefront  as 
measured  in  the  radial  plane.  This  effect  likewise  Increases  the  radius  of 
curvature  at  a  high-  to  low-velocity  Interface.  Again,  the  attenuation  rate 
In  the  second  layer  will  be  decreased.  The  change  in  divergence  of  the  wave- 

front  approaches  a  limiting  value  for  normal  Incidence,  defined  by  the  velo¬ 
cities  In  the  two  media. 

Multilayer  wave  propagation  thus  Involves  two  distinct  phenomena: 
transmission  at  an  Interface,  and  attenuation  of  the  transmitted  wave  beyond  the 
Interface.  The  problem  of  Interactions  at  the  Interface  for  a  specific  ray 
may  properly  be  analyzed  with  a  one-dimensional  code.  However,  the  geometrically- 
conditioned  problem  of  shock  propagation  Into  the  upper  layer  Is  Inherently 
two-dimensional.  The  spherical  divergence  condition  Imposed  In  one  dimension 
excessively  damps  the  shock  energy  and  Increases  attenuation  (If  the  wavefront 
divergence  should  decrease,  as  at  a  high-  to  low-velocity  interface), or  coventrates 
the  shock  energy  and  decreases  attenuation  (If  divergence  should  Increase,  as 
at  a  low-  to  high-velocity  Interface).  In  the  case  of  rock  calculations  with 
"low- velocity"  water  overburden,  divergence  of  the  shock  front  decreases 
dramatically  at  the  interface  (radius  of  curvature  Increases).  The  SOC  code 
should  therefore  underpredict  shock  pressures  and  velocities  in  the  upper  layer. 


-32- 


Unllke  spherical  convergence,  which  develops  only  for  reflected  waves, 
this  effect  applies  to  the  Initial  outgoing  shock  wave.  Serious  errors 
may  occur  In  the  first-order  dynamic  parameters  and  In  safety  or  effects 
predictions  for  the  upper  layer.  Initial  shock  front  peak  velocities  and 
pressures  moving  Into  the  water  layer  will  be  progressively  reduced  (due 
to  excessive  divergence  of  the  outgoing  shock  wave  beyond  the  Inter¬ 
face),  and  the  peak  spall  velocity  at  the  water  surface  will,  of  course,  be 
underpredicted.  One  consolation  lies  In  the  fact  that  tensile  wave  spall 
velocity  calculations  below  the  water  surface  will  be  Improved.  This  Improve¬ 
ment  occurs  for  two  reasons:  (1)  Spherical  convergence  of  the  water  surface 
reflection  In  one-dimensional  calculations  causes  spall  velocities  to  be 
overestimated;  however,  divergence  of  the  (two-dimensional)  compression  and 
reflected  tension  waves  In  the  water  layer  Is  reduced  by  the  refraction 
effect.  The  error  caused  by  spherical  convergence  Is  thereby  decreased; 

(2)  The  one-dimensional  error  caused  by  Ignoring  refraction  decreases  velo¬ 
cities  In  the  outgoing  compresslonal  shock  (both  the  peak  and  decaying 
velocities  In  the  "tall"  of  the  pulse).  The  error  created  by  spherical 
convergence  Increases  the  particle  velocities  for  the  reflected  tension  wave 
below  the  water  surface.  The  combined  effect  of  these  errors  on  the  reinforced 
"tensile  wave  plus  remaining  compresslonal  wav*1  spall  velocities  Is  partially 
self-cancelling  near  the  middle  of  the  water  layer  (but  not  at  the  Interface 
or  the  surface). 

A  detailed  examination  has  been  conducted  to  determine  the  effect  of 
refraction  errors  on  the  peak  shock  parameters.  Water  was  treated  as  a 
linearly  elastic  medium  with  acoustic  propagation,  permissible  In  this  case 
because  peak  pressures  In  the  water  layer  are  on  the  order  of  one  kl lobar. 

Refraction  at  the  Interface  and  transmission  In  the  water  layer  were  calculated 
using  the  accoustlc  formulae  for  peak  pressure  of  a  simple  compresslonal  wave.^®^ 
Shear  wave  effects  and  shear  reflection  at  the  Interface  were  not  taken  Into  account. 
It  was  found  that  the  worst  errors  occur  for  large  velocity  mismatch  at  the 
Interface  and  for  deep  water  layers  (long  ranges  from  the  Interface). 

Sample  results  for  rock  and  water  layers  of  equal  thickness  (coral  36/36, 
basalt  36/36)  are  presented  In  Table  3.  The  SOC-calculated  peak  spall  velo¬ 
cities  are  compared  with  spherical ly-dlvergent  acoustic  calculations  (uncorrected 
for  refraction  effects  on  wavefront  attenuation  away  from  the  Interface)  and 
corrected  acoustic  calculations  (refraction  effects  Included).  The  spherical 


;  of  Refraction  Error  on  One-Dimensional 
Calculated  Spall  Velocities. 
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(uncorrected)  acoustic  velocities  are  similar  to  the  SOC  results.  However, 
the  corrected  acoustic  velocities  at  the  water  surface  are  higher  by  about 
21%  (for  a  slight  velocity  mismatch,  coral-water),  and  56%  (large  velocity 
mismatch,  basalt-water).  The  peak  velocities  and  pressures  at  the  shock 
front  are  In  error  by  a  similar  amount  at  ranges  well  beyond  the  rock-water 
Interface.  For  an  extreme  velocity  mismatch,  the  pressures,  particle  velo¬ 
cities,  and  spall  velocities  far  from  the  Interface  may  be  Incorrect  *>y  a 
factor  of  two  or  more.  Results  of  the  Individual  calculations  are  presented 
In  a  later  section.  The  refraction  error  explains  some  discrepancies  which 
have  appeared  In  earlier  applications  of  one-dl mens  Iona 1  codes  to  layered 
configurations.  This  effect  Is  particularly  troublesome,  since  shock  para¬ 
meters  are  underpredicted  for  the  commonly-encountered  situation  of  decreasing 
sonic  velocity  In  the  shallow  (near-surface)  layer. 

The  combined  Influence  of  refraction  and  spherical  convergence  errors 
on  multilayer  calculations  has  not  yet  been  fully  Investigated.  The  situation 
Is  rather  complex  because  the  dynamic  effects  of  the  compresslonal  shock  and 
of  the  rock-water  Interface  reflection  (below  the  rock  surface)  have  already 
occurred  before  the  water  surface  reflection  arrives.  Also,  recompaction 
effects  near  the  rock-water  Interface  and  the  cavity  gas  acceleration  signal 
(which  Is  accentuated  by  the  water  layer  recompaction  of  underlying  material) 
may  precede  the  water  surface  reflection  In  some  regions.  Acoustic  analysis 
Is  Inappropriate  due  to  these  late-time  effects.  In  general,  the  same 
arguments  applied  to  spherical  convergence  errors  are  appropriate  here.  The 
geometrically-caused  errors  may  safely  be  Ignored  for  small  changes  In 
distance  (l.e.,  the  shallow  water-layer  calculations,  such  as  the  36/05  or 
36/10  configurations).  However,  the  shock  front  parameters  and  the  dynamics 
of  an  Individual  zone  may  be  greatly  modified  If  a  large  change  In  problem 
geometry  occurs  (36/36  calculations).  These  changes  may  be  rather  more  apparent 
for  the  refraction  error,  since  the  dynamics  of  the  outer  layers  (rather 
than  the  less  massive  inner  regions)  are  first  affected.  The  overall 
problem  dynamics  are  less  sensitive  to  these  errors  than  are  the  velocities 
at  a  given  location.  Therefore,  mound  kinetic  energies  and  residual  velo¬ 
cities  throughout  the  mound  after  the  shock  fronts  have  passed  provide 
appropriate  criteria  for  comparing  the  Individual  calculations. 
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Another  fortunate  situation  also  reduces  the  Influence  of  geometric 
errors  on  rock  dynamics.  Refraction  difficulties  are  most  significant  for 
a  large  velocity  mismatch  at  the  Interface  (1.e.,much  stlffer  underlying 
medium).  However,  the  reflection  at  the  Interface  Is  also  strongest  for 
these  cases.  In  the  extreme  situation  of  a  very  low-impedance  overlying  material, 
the  shock  Is  almost  perfectly  reflected  as  a  tensile  wave.  Thus,  when  a  strong 
impedance  mismatch  exists,  tensile  wave  spallation  due  to  the  Interface 
reflection  will  be  predominant  In  the  rock  layer.  Likewise,  the  cavity  gas 
acceleration  signal  and  recompaction  effects  propagate  more  rapidly  In  high- 
velocity  materials.  These  effects  may  have  more  Influence  on  problem  dynamics 
near  or  below  the  Interface  than  the  weak  water  surface  reflection.  In  the 
36/36  calculations.  It  was  found  that  rock  layer  spallation,  water  layer 
recompaction  effects,  and  the  later  gas  acceleration  pulse  control  rock  dynamics 
near  the  Interface.  For  coral  36/36  and  basalt  36/36,  the  water  surface  reflection 
arrives  at  the  Interface  Immediately  before  the  enormously  larger  gas  accelera¬ 
tion  pulse.  Its  only  effect  Is  to  terminate  the  very  slight  deceleration  due 
to  water  layer  recompaction  which  Is  still  occurring  at  this  time.  For  granite 

36/36,  the  gas  acceleration  pulse  arrives  at  the  Interface  before  the  water  surface 
reflection  can  return,  and  overwhelms  water  tensile  reflection  effects.  Thus, 

dynamics  In  the  rock  layer  (and  In  the  lower  portion  of  the  water  layer)  are 
dependent  upon  the  Initial  shock  wave.  Interface  spall  reflection,  recompaction 
prior  to  the  water  surface  reflection,  and  the  cavity  gas  acceleration  pulse. 
Long-range  propagation  and  surface  reflection  of  the  water  pulse  Influence 
dynamics  only  In  the  upper  portion  of  the  water  layer,  and  are  unimportant  In 
the  lower  regions  of  the  problem.  The  water  surface  effects  are  more  signifi¬ 
cant  for  very  shallow  water  layers,  but  geometric  errors  In  such  problems  become 
negligible.  This  conclusion  Implies  that  precise  modelling  of  the  "upper  layer" 
shock  propagation  Is  not  crucial  to  the  "lower  layer"  dynamics  In  the  Impedance- 
mismatched  problems  examined  here.  One-dimensional  codes  may  be  legitimately 
used  for  this  particular  class  of  problems.  Equally  Important,  submerged 
cratering  dynamics  at  early  times  are  controlled  by  the  Interface  spall  phenomenon, 
water  layer  recompaction,  and  the  consequent  enhancement  of  the  gas  acceleration 
pulse.  For  this  reason,  a  single-layer  analysis  with  one  Interface  cannot 
adequately  describe  the  process;  multilayer  calculations  of  the  type  employed 
here  must  be  used. 
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V.  MATERIAL  PROPERTIES 

Hydrodynamic  calculations  require  a  knowledge  of  the  constitutive 
relations  for  all  materials  In  the  problem.  Ideally,  these  relations  should 
be  established  over  the  entire  range  of  expected  physical  effects.  Five 
materials  are  encompassed  In  this  study:  Aluminized  anmonlum  nitrate  slurry 
explosive,  water,  and  the  saturated  rock  media  coral,  basalt,  and  granite. 

The  properties  of  these  materials  will  be  examined  separately. 

Explosive.  Detonation  characteristics  and  constitutive  relations 
of  a  wide  variety  of  high  explosives  have  been  extensively  Investigated.  A 
convenient  formulation,  known  as  the  Jones-WIlklns-Lee  equation  of  state. 

Is  often  used  In  hydrodynamic  calculations.^  The  best-studied  materials 
fall  within  the  class  of  "military"  or  high  explosives,  which  propagate  a 
nearly-dlscontlnuous  detonation  front  at  constant  velocity.  Energy  Is  released 
Inmedlately  upon  arrival  of  the  front.  The  properties  of  commercial  explosives 
available  for  civil  excavation  projects  are  somewhat  less  well-defined.  Only 
one  commercial  explosive,  an  aluminized  ammonium  nitrate  slurry,  has  been  studied 
In  detail.*-  •*  A  modified  version  of  the  Jones-WIlklns-Lee  equation  of  state 
has  been  developed  to  describe  this  material.  The  equation  of  state  assumes 
a  total  energy  release  of  0.11  Mbar-cm3/cm3  (1752  g-cal/g).  About  23*  of  the 
explosive  energy  Is  released  Instantaneously  upon  arrival  of  the  detonation 
front  at  a  given  location.  The  remaining  part  of  the  energy  Is  released  at  a 
constant  rate  over  a  period  of  0.025  msec  after  the  detonation  front 
arrival  time.  An  explosive  sphere  of  l.T3-m  radius  was  used  In  these 
calculations.  The  properties  of  the  energy  source  are  as  follows: 

P0  (Initial  density)  -1.50  g/cm3 
Det.  Velocity  ■  5300  m/sec 
Radius  »  1.13  m 

Total  energy  Input  ■  6.65  x  105  Mbar-cm3 

-  1.59  x  1010  g-cal 
*  15.9  tons  of  explosive  energy 
Weight  •  9.066  x  106  g  =  10  tons. 

The  explosive  sphere  was  divided  Into  30  zones  of  equal  radius;  calculatlonal 
time  step  was  restricted  so  that  no  more  than  one  percent  of  the  explosive  energy 
could  be  generated  per  cycle.  With  this  model,  the  peak  of  the  explosive  shock 
wave  entered  the  surrounding  medium  about  0.24  msec  after  the  beginning  of  the 
calculation. 
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Water.  Compressibility  and  tbiarmodynamlc  properties  of  water  and  Ice 
are  well-established  over  a  very  wide  range  of  pressures  and  temperatures. f-3l'38j 
Available  data  have  been  used  to  produce  the  two  water  equation  of  state  curves 
(pressure  vs. specific  volume)  shown  In  Figures  3,  4,  .and  5.  The  upper  curve 
In  Figures  3  and  4  (larger  specific  volume  at  low  pressures)  Is  based  on  a 
measured  Hugonlot  centered  at  20°  C  and  one  atmosphere. ^1*32] 
equation  of  state  was  adopted  for  the  water  layer  calculations.  The  second 
curve  was  used  by  Butkovlch  In  a  model  for  generating  equations  of  state  of  a 
rock-water-air  matrlx^3^;  It  was  applied  In  determining  the  basalt  and  granite 
equations  of  state  for  this  study  (see  below).  The  second  curve  has  no  rela¬ 
tionship  to  the  water  layer  calculations.  The  water  equations  of  state  do  not 
Include  hydrostatic  water-ice  phase  changes  near  11  and  24  kbar;  It  Is  not 
known  whether  these  transitions  occur  under  shock  loading. ^  Peak  pressures 
In  the  water  layer  are  on  the  order  of  1  kbar  for  these  calculations,  so 
phase  changes  would  not  take  place  In  the  water.  Shock  vaporization  of  water 
was  likewise  Ignored,  since  vaporization  occurs  only  at  pressures  above  50  to 

It  Is  noted  that  the  water  layer  compressibility  curve  shows  a  slight 
discontinuity  at  low  pressure  (Figure  3).C31»32^  The  abrupt  Increase  and 
subsequent  decrease  In  compressibility  cause  corresponding  sudden  changes  In 
the  bulk  modulus  and  shock  velocity  for  the  water.  This  behavior.  In  turn, 
gives  rise  to  a  discontinuity  In  the  hydrodynamic  shock  waveforms  at  low  pres¬ 
sures.  The  equation  of  state  has  been  modified  to  achieve  constant  bulk  modulus 
at  low  pressures  (dotted  line.  Figure  3).  The  Initial  bulk  modulus  and  sonic 
velocity  for  the  modified  curve  correspond  with  the  expected  values  for  water 
at  about  20°C.  Calculations  were  attempted  using  both  equations  of  state; 
the  constant  bulk  modulus  eliminated  low-pressure  discontinuities  In  the  water 
waveforms,  but  no  significant  effect  on  the  rock  layer  was  apparent.  Results 
presented  here  are  based  on  the  constant  bulk  modulus  calculations. 

Rock  Media.  The  three  rock  media  were  selected  to  typify  potential 
underwater  excavation  sites  and  to  represent  a  wide  range  of  material  properties. 
The  sites  which  served  as  precise  or  approximate  analogs  for  these  media  are 
listed  below. 


Figure  3«  Equation-of- ufatc  of  v/auu-,  0  to  1  kb  a  • 


Pressure  (kbar) 
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Figure  5.  Equation-of-state  of  water,  0  to  500  kbar 
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1.  Coral:  Tugboat  coral,  Kawalhae  Harbor,  Hawaii;  site  of 
Project  Tugboat,*-9**  large-scale  underwater  explosive  harbor  excavation 
series.  Material  Is  a  porous,  low  strength,  fully-saturated  aragonltlc 
coral.  Engineering  properties  data,  pressure-volume  measurements,  and  failure 
strength  tests  are  available.^9*40]  Tests  were  made  on  Intact  sanples; 

the  overall  density  of  the  reef  Is  slightly  lower  and  water  content  Is 
higher  due  to  branching  of  the  coral.  On-site  measurements  Including 
seismic  velocities  and  experimental  data  Including  alrblast  measurements 
and  spall  velocities  from  explosive  events  have  also  been  obtained.  [Ml  ,42] 

2.  Basalt:  Kaalualu  Harbor  basalt,  Island  of  Hawaii;  potential 
site  of  a  small  boat  harbor.  The  material  Is  a  porous,  low-  to  Intermediate- 
strength  saturated  basalt  representative  of  possible  excavation  sites  In 
competent  rock.  Samples  and  engineering  properties  data  available.  No 
compressibility  or  strength  tests. 

3.  Granite:  Orinoco  River  high-strength  granite,  Orinoco  River, 
Venezuela;  site  of  proposed  underwater  excavation  project.  The  medium  Is 

a  very  dense,  low  porosity  high-strength  granite.  Samples  and  engineering 
properties  data  available.  No  compressibility  or  strength  tests. 

Required  Input  data  for  each  of  these  media  Include  the  pressure-volume 
compressibility  curve  or  equation  of  state  (EOS),  which  specifies  the  bulk 
modulus  as  a  function  of  pressure.  One  additional  modulus  Is  required  to 
specify  material  behavior.  In  these  calculations,  a  constant  Poisson's  ratio 
has  been  used  for  each  medium.  Finally,  failure  envelopes  are  required  to 
determine  the  brittle  and  ductile  failure  behavior  of  the  material  under 
compression,  shear,  and  tension.  Derivation  of  these  Input  parameters  Is 
discussed  below. 

The  coral  medium,  as  previously  noted.  Is  a  very  low  strength  material 
subject  to  extensive  crushing,  compaction,  and  washback  effects.  Final  crater 
sizes  cannot  be  predicted  with  currently  available  codes.  Nonetheless,  coral 
Is  Included  In  the  calculation  suite  because  of  the  availability  of  laboratory 
test  data  and  limited  dynamic  measurements  from  Project  Tugboat.  The  Tugboat 
excavation  events  utilized  cylindrical  charges  of  an  amnonlum  nitrate  slurry 
explosive.  The  emplacement  configuration  Is  compared  with  a  schematic  "SOC" 
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vertical  problem  In  Figure  6.  Note  that  the  Tugboat  explosive  source 
Is  less  dense,  and  therefore  greater  In  volume,  than  the  AANS  explosive 
used  In  the  calculations.  The  L.ieclfled  energy  (z  1700  g-cal/g)  and 
detonation  velocity  (z  5000  m/sec)  of  the  Tugboat  charges  are,  however, 
closely  comparable  to  the  calculations.  The  Tugboat  Phase  I  (single-charge) 
events  will  be  compared  with  SOC  results  In  later  sections  of  the  paper. 

Engineering  properties  data  for  numerous  coral  samples  were  derived 
during  the  course  of  the  experiments,^  while  hydrostatic  pressure- volume 
curves  and  compresslonal  failure  data  were  taken  by  Stephens. Coral 
material  parameters  and  elastic  moduli  derived  from  these  data  are  listed  In 
Table  4.  Stephens'  hydrostatic  loading  and  unloading  data  (for  loading  to 
30  (kbar)  are  presented  In  Table  5,  and  plotted  In  Figures  7  and  8.  The 
loading  and  unloading  curves  coincide  at  low  pressure,  as  would  be  expected 
for  a  saturated  material  (no  air  void  porosity  to  be  removed  by  compression). 
However,  the  effects  of  the  water  phase  transitions  are  apparent  from  10  to 
12  kbar  and  18  to  25  kbar.  These  transitions,  which  may  not  occur  under  shock 
loading,  generate  slight  discontinuities  In  the  hydrodynamic  calculation  wave¬ 
forms.  Therefore,  a  "smoothed"  equation  of  state  curve  was  drawn  through  loading 
and  unloading  data  (Table  5  and  Figure  8).  Horizontal  SOC  calculations  were 
performed  with  the  measured  and  smoothed  data  (PS10HRZ  and  CORHRZ,  respectively). 
Differences  between  the  two  calculations  were  small  and  consisted  primarily  of 
noise  In  the  waveforms  close  to  the  cavity.  The  smoothed  curve  was  used  In 
all  subsequent  calculations. 

Stephens  did  not  acquire  loading  data  above  30  kbar,  and  no  very  high 
pressure  Hugonlot  measurements  are  available  for  coral.  Therefore,  the  pressure- 
volume  curve  has  been  extrapolated  to  higher  pressures  (Table  5  and  Figure  9). 

The  extrapolation  Is  partially  based  on  measured  Hugonlot  data  for  Sabana  shale 
(Figure  9),  a  very  weak  saturated  material  with  a  somewhat  higher  density  and 
lower  porosity  than  coral,  and  on  data  for  other  saturated  shale  media. C43-45»3^ 
Hydrostatic  and  Hugonlot  data  are  also  available  for  dense,  nonporous  calcium 
carbonate  and  magnesium  carbonate  rocks.  These  data  derive  In  part  from  shock 
experiments  using  marble,  limestone,  and  other  carbonate  media  with  a  compo¬ 
sition  similar  to  the  coral  matrix,  but  with  little  or  no  air  void  porosity 
or  water  content.  The  compressibility  of  a  coral -water  matrix  may  be 
simulated  by  using  a  calcium  carbonate  equation  of  state  to  describe  the 
rock  matrix  and  a  normal  water  compressibility  curve  for  the  water  component 
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Figure  6.  Emplc cement  configurations  for  10-ton  (20,000  lb)  charges 
TUGBOAT  detonations  and  spherical  S0C  calculations. 


TABLE  4-  Measured  and  Assumed  Material  Parameters  for  Saturated  Submerged  Cratering  Media  (Saturation >99.91  in  an  v..— 
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TABLE  5.  Measured,  Smoothed  and  Extrapolated  Equation  of  State  Data 
for  Saturated  Tugboat  Coral. 


Pressure 

P(Mbar) 

Measured 

Loading  Curve 

Spec.  Vol . 

(cm3/q) 

Smoothed 
Loading  Curve 
Spec.  Vol. 

(cm3/q) 

Measured 
Unloading  Curve 
Spec.  Vol. 

(cm3/a) 

0 

.531915 

.531915 

.531915 

.001 

.521277 

.521277 

.521277 

.002 

.512341 

.512341 

.512341 

.003 

.504894 

.504894 

.504894 

.005 

.492766 

.492766 

.492766 

.007 

.483245 

.483245 

.483245 

.010 

.471436 

.471436 

.470213 

.011 

.467926 

./!S7926 

.461064 

.012 

.458351 

.4640 

.458351 

.015 

.451755 

.4550 

.451755 

.018 

.446649 

.446649 

.446649 

.020 

.443617 

.4402 

«•  •  •» 

.022 

.440851 

.4345 

.426436 

.025 

.423138 

.4264 

.423138 

.027 

.420745 

.4218 

— 

.030 

.417553 

.417553 

— 

.040 

.406873* 

.060 

.387038* 

— 

.080 

.373814* 

— 

.100 

.361000* 

— 

.150 

.334144** 

— 

.200 

.315631** 

— 

.250 

.302305** 

— 

.300 

.291931** 

— 

.400 

.275148** 

— 

* 

Extrapolation,  30  -  100  kbar 

Based  on  Hugonlot  data  for  saturated  Sabana  shale,  normalized  to  the 
measured  specific  volume  of  coral  at  30  kbar. 
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Figure  8.  Hydrostatic  loading  and  unloading  curves 
and  smoothed  EOS  curve  for  TUGBOAT  coral , 
0  to  30  kbar. 
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(slmple  mixture  volumetric  compressibility  model).  This  analysis  predicts 
stlffer  bulk  moduli  and  lower  compressibility  than  Figure  9  would  Indicate, 
both  In  the  low  pressure  regime  (measured  data)  and  along  the  high  pressure 

Hugonlot  (extrapolation).  Volume  errors  of  several  percent  occur  In  the  region 
of  the  hydrostatic  measurements.  However,  the  bulk  modulus  of  coral 

Increases  significantly  at  the  upper  end  of  the  measured  pressure  range 

(25-30  kbar).  Thus,  coral  m^y  be  less  compressible  than  the  compromise  fit  in 
Figure  9  at  high  pressures.  If  the  stiff  behavior  Indicated  by  the  carbonate 
Hugonlot  Is  correct,  the  extrapolation  will  be  Increasingly  Inaccurate  at 
higher  pressures  $  100  kbar.  Since  maximum  pressures  encountered  furlng  high- 
explosive  calculations  are  on  the  order  of  100  kbar,  the  Implied  Inaccuracies 
are  not  Important  here. 

Strength  data  for  the  Tugboat  coral  material  were  obtained  to  a 
maximum  confining  pressure  of  3.5  kbar  by  trlaxlal  compression  and  Brazil 
tests  on  Intact  core  samples. Results  are  displayed  In  Figure  10. 

The  abscissa  and  ordinate  of  the  plot  conform  with  the  parameters  used  In 
the  code:  failure  shear  stress  k  [■  (a1 -a3)/2]  and  adjusted  mean  pressure 
F  [-(o^+o3)/2].  Ho  data  were  obtained  for  failed  (crushed)  material, 
but  an  estimated  failure  envelope  has  been  drawn  below  the  "Intact" 
data  (and  assuming  no  tensile  strength  for  failed  material).  The  data  shown 
In  Figure  10  Indicate  an  abrupt  Increase  In  strength  to  about  60  bars  at 
very  low  pressures,  followed  by  an  apparent  slight  decrease  at  higher  F  values. 
This  unusual  behavior  has  been  explained  as  a  work-hardening  effect  Involving 
collapse  of  the  coral  matrix. ^  The  code  does  not  accept  a  decreasing  strength 
function,  so  a  constant  strength  Is  assumed  between  0.2  and  2.86  kbar  (dashed) 
line  In  Figure  10).  The  failure  surfaces  have  been  extrapolated  to  a  maximum 
shear  strength  of  101  bars  at  and  beyond  5  kbar.  Tests  Indicate  that  the  behavior 
of  coral  Is  ductile  except  at  very  low  confining  stresses;  a  transition  from 
brittle  failure  to  ductile  behavior  at  0.2  kbar  Is  assumed. 

Engineering  data  for  Kaalualu  basalt  and  Orinoco  River  granite  are 
presented  In  Table  4.  The  measurements  Include  bulk  density,  rock  grain  density, 
and  porosity. The  elastic  moduli  and  other  properties  not  available  from 

fl 7  47-51 i 

measurements  were  estimated  using  data  for  similar  media.1  *  J  Decause  the 
elastic  constants  of  an  Isotropic  linearly  elastic  material  are  Interrelated, 
the  density  and  any  two  constants  are  sufficient  to  define  the  material  (see 
"Relationships  Among  Elastic  Constants",  table  at  the  beginning  of  this  report). 
Such  parameters  as  the  bulk  modulus,  congressional  velocity,  shear  modulus. 
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and  Poisson's  ratio  have  been  determined  for  various  basalts  and  granites. 

These  measurements.  In  conjunction  with  reasonable  expectations  for  the  known 
sample  densities,  have  been  used  to  produce  the  values  In  Table  4.  In  parti¬ 
cular,  the  bulk  moduli  and  sonic  velocities  are  known  to  be  reasonable  for 
comparable  geologic  materials.  The  shear  modulus  (or  Poisson  ratio)  is  less 
accurately  fixed.*  The  Poisson's  ratio  for  porous  basalts  does  not  vary 
dramatically,  and  the  assumed  value  Is  probably  accurate.  However,  a  very 
low  Poisson  ratio  (high  shear  modulus)  Is  accepted  for  the  granite.  This 
Is  In  accord  with  the  relatively  high  bulk  modulus  and  sonic  velocity.  Hardhat 
granite  has  been  used  as  a  model  for  the  elastic  parameters  In  this  case.^^ 

The  Hardhat  material  Is  a  very  dense  high-strength  rock  from  a  deep  layer 
(z  1000  ft).  It  has  low  porosity,  with  a  high  sonic  velocity  and  shear  modulus. 
The  granite  Is  thus  being  selected  to  represent  a  dense,  highly  competent 
medium  at  one  extreme  of  the  material  properties  spectrum.  This  choice  may  not 
be  Inappropriate,  In  view  of  the  low  porosity  and  probable  high  strength  of 
the  Orinoco  rock. 

Compressibility  curves  and  failure  surfaces  for  the  basalt  and  granite 
are  likewise  unknown.  Butkovlch  has  recently  developed  a  method  for  estimating 
the  equation  of  state  of  silicate  rocks. E39»52»53]  jhl s  technique,  which  Is 
based  on  simple  mixture  theory  for  a  rock-water-air  matrix,  generates 
accurate  pressure-volume  relations  for  a  variety  of  silicate  rocks,  over  a 
broad  range  of  porosity  and  water  content.  The  equation  of  state  for  the 
grain  density  rock  Includes  allowance  for  high-pressure  (Hugonlot)  phase  transi¬ 
tions  known  to  occur  In  silicate  minerals.  The  water  equation  of  state 
discussed  previously  is  used  for  the  water  component.  A  slightly  modified 
version  of  the  Butkovlch  technique,  suitable  for  fully-saturated  materials, 
has  been  used  to  generate  basalt  and  granite  equations  of  state.  Following 
empirical  expectations  and  experimental  evidence  for  saturated  media, ^  the 
bulk  modulus  Is  assumed  to  be  slowly  and  monotonlcally  Increasing  at  low 

*  The  shear  modulus,  for  example,  Is  highly  sensitive  to  the  assumed 
sonic  velocity.  Selsmlcally  determined  sonic  velocities  may  differ  from 
values  measured  for  small  intact  rock  samples,  thus  causing  significant 
uncertainty  if  shear  modulus  is  calculated  from  this  parameter. 
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pressures  (%  0  to  2  kbar).  The  low-pressure  curve  Is  then  smoothly  connected  with 
a  rock-water  matrix  equation  of  state  at  higher  pressures.  Loading  and 
unloading  curves  are  assumed  to  be  Identical  (no  nonrecoverable  air  void 
porosity).  The  resultant  equations  of  state  for  basalt  and  granite  are 
displayed  In  Figures  11-12  and  13-14,  respectively.  In  this  case,  pressure 
has  been  plotted  against  the  dimensionless  volumetric  compression  V 
(v  *  f”  *1“  rr  -1;  where  p  ■  density,  or  v  =  specific  volume).  Data 
for  similar  basalts  and  granites  are  presented  for  comparison.  The  basalt  curve 
appears  consistent  with  measurements  of  comparable  porous  and  consolidated 
basalts  (pQ  %  2.6)  at  low  pressures.  It  falls  somewhat  below  (less  stiff  than) 
dense  consolidated  basalt,  and  the  canal  zone  basalts. ^3  fa  very  high 
pressures  (JJj  300  kbar),  the  large  water  content  causes  the  generated  curve  to 

become  less  stiff  than  the  measured  basalts.  This  region  Is  of  no  Interest 
for  high-explosive  calculations.  The  granite  curve  compares  closely  with 
measured  data. 

Failure  surfaces  have  also  been  estimated  on  the  basis  of  available 
<?ata.  The  basalt  medium,  as  Indicated  by  extensive  measurements  for  porous 
and  saturated  basalts, £17*31^  Is  allowed  a  tensile  strength  of  0.1  kbar  and  an 
ultimate  strength  of  0.4  kbar.  The  granite  C31»54*493  is  arbitrarily  assumed 
to  have  a  high  ultimate  strength  of  5  kbar,  with  a  tensile  strength  of  0.38  kbar. 
There  Is  some  evidence  that  strength  may  decrease  considerably  with 
saturation,  even  In  the  case  of  a  very  low  porosity  granite. t48»143 

However,  the  high  strength  Is  retained,  In  keeping  with  the  desire  for  one 
"extreme  case"  calculation.  Failure  envelopes  of  the  assumed  ultimate  strength 
have  been  smoothly  fitted  using  an  empirical  hyperbolic  equation. ^39^  The 
resultant  failure  surfaces  for  Intact  rock  are  shown  In  Figures  15-16  (basalt) 
and  17-18  (granite).  Low-pressure  failure  surfaces  for  crushed  material  have 
been  drawn  at  approximately  one-half  of  the. corresponding  failure  strengths  for 
Intact  material.  The  condition  of  no  tensile  strength  at  F=0  has  been  Imposed 
at  the  lower  end  of  these  curves,  as  Is  reasonable  for  fully-crushed  rock.  The 
crushed  failure  surfaces  are  then  smoothly  connected  with  the  Intact  failure 
envelopes  at  very  high  pressures. 
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Figure  13.  Pressure  vs  compressibility  for  dense  granites 
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The  parameters  which  control  multilayer  acoustic  wave  transmission 
are  the  Impedances  and  sonic  velocities  of  the  layers.  Acoustic  parameters 
for  the  five  media  discussed  above  are  summarized  In  Table  6.  The  ratio  of 
the  medium  Impedance  to  the  Impedance  of  water  "P0C/pCwa^er"*  and  the  pressure 
transmission  factor  normal  to  the  medium-water  Interface  "Tpgo,"  are  also  9lven. 
Acoustic  transmission  parameters  are  not  valid  at  a  strong  shock  front. 

However,  typical  pressures  near  the  rock-water  Interface  are  on  the  order  of 
one  kllobar  In  these  calculations;  the  acoustic  transmission  analysis  may 
be  used  with  a  high  degree  of  accuracy. 

VI.  RESULTS:  ARRIVAL  TIMES,  VELOCITIES,  STRESSES,  AND  WAVEFORMS 

This  section  presents  detailed  results  of  the  "SOC"  horizontal 
and  vertical  calculations.  Since  the  study  Involved  a  broad-scale  survey 
of  several  media  and  configurations,  any  one  of  which  might  prove  applicable 
to  a  future  design.  It  Is  considered  worthwhile  to  compile  a  fairly  extensive 
library  of  Individual  events.  Accordingly,  the  peak  velocities,  pressures,  and 
stresses  have  been  plotted  as  a  function  of  range  for  all  calculations.  Selected 
velocity  and  pressure  time-dependent  waveforms  at  several  ranges  from  the  center 
of  each  problem  are  also  presented. 

Before  reviewing  the  dynamic  results.  It  Is  appropriate  to  consider 
the  order  of  effects  of  Interest.  We  are  examining  typical  cratering  events  of 
15.9  tons  energy  yield,  at  burial  depths  of  about  36  ft  (143  ft/kt1/3).  Shock 
parameters  at  ranges  In  this  vicinity  are  relevant  to  crater  formation.  For 
the  vertical  calculations,  material  velocities  near  tiie  free  surface  (or  rock- 
water  Interface)  are  a  primary  factor  In  material  ejection.  Thus,  velocity 
waveform  changes  near  the  Interface  attributable  to  the  water  layer  are 
significant.  Peak  particle  velocities  near  the  Interface  are  typically  on 
the  order  of  70  m/sec  (l.e.,  near-optimum  burial  depth  cratering  events); 
peak  pressures  near  the  Interface  and  In  the  water  layer  are 
%  1  kbar.The  Interface  region  Is  In  the  range  of  assumed  brittle  failure 
behavior  for  basalt  and  granite,  but  Is  ductile  for  coral.  The  ultimate  shear 
strength  of  5  xbar  for  granite  Is  never  reached  at  the  Interface,  or  at  ranges 
$  4  m  from  the  center  of  the  problem.  The  ultimate  strength  of  basalt  may 
or  may  not  be  exceeded  In  the  vicinity  of  the  Interface.  Finally,  the 
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TABLE  6.  Acoustic  Parameters  for  the  Media  and  for  Aluminized  Ammonium  Nitrate 
Slurry  Explosive  (AANS). 


Medium 

K 

Initial 

Bulk 

Modulus 

(kbar) 

po 

‘  Impedance  * 

poC 

OoC 

;•  ★ 

TP90 

(medium- 

water) 

(g/cm3) 

/psl-sec, 

(  ft  > 

pC(water) 

Coral 

49 

6,500 

1 

1.880 

164.53 

2.4816 

0.5744 

Basalt 

210 

11,810 

2.605698 

414.32 

6.249 

0.2759 

Granite 

380 

17,800 

2.718234 

651.43 

9.825 

0.1847 

Water  (fresh. 

Temp.  Rr  20-30°C) 

22.568 

4,933.4 

0.9982 

66.30 

1.0 

(1.0) 

AANS  explosive** 

Cdet 

=17,389 

1.500 

^det 

“351.18 

(5.2968) 

(0.3176) 

*  Tp90  =  Pressure  transmission  factor  for  an  acoustic  wave  propagating  normal  to  the 
medium  -  water  interface. (Reflection  factor  RpgQ  =  1  -  Tpgg;  since  water  is 
of  lower  Impedance  than  the  other  media,  the  transmitted  wave  is  compress  Iona  1 
and  the  reflected  wave  Is  tensile.) 

**  The  AANS  material  used  In  these  calculations  Is  a  high-energy  non-ideal  aluminized 
ammonium  nitrate  slurry  explosive.  The  velocity  of  the  propagating  detonation  front 
"Cdet"  is  17,389  ft/sec  (5300  m/sec);  the  total  energy  release  is  assumed  to  be 
0.11  Mbar-cm3/cm3.  The  explosive  Is  described  in  greater  detail  in  Chapter  V. 
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ultlmate  strength  of  0.101  kbar  for  coral  is  exceeded  In  virtually  all 
zones,  throughout  the  vertical  problems.  It  Is  of  Interest  to  note  that 
peak  pressures  at  or  near  the  Interface  lie  In  the  relatively  low  pressure 
range,  where  a  slowly-changing  bulk  modulus  Is  assumed  for  saturated  materials 
and  water.  Thus,  nearly  linear-elastic  behavior  for  peak  transmission  at  and 
beyond  the  Interface  may  be  expected.  The  peak  shock  parameters  will  be 
analyzed  acoustically  and  compared  with  SOC  runs  In  Chapter  VIII. 

Time  of  Arrival .  Shock  time  of  arrival  measurements  provide  a  useful 
means  of  checking  code  calculations.  Arrival  times  have  been  obtained  for 
only  one  underwater  cratering  event,  Tugboat  Phase  II  UKL  (a  four-charge 
array  of  10-ton  charges;  Dr  *  36  ft,  Dw  *  5  ft).  Several  dynamic 
gages  were  emplaced  In  the  water  layer,  directly  above  each  of  the 
four  charges.  Little  usable  Information  was  recovered  from  these  gages, 
but  the  Initial  arrival  was  detected  In  most  cases.  Unfortunately,  the 
recording  system  used  did  not  give  accurate  "real  time"  measurements  (detonation 
to  shock  arrival),  and  the  locations  of  the  gages  were  not  very  precisely 
determined. 

Vertical  propagation  from  each  charge  of  the  Tugboat  IJKL  event 
may  be  compared  with  the  almost-ldentlcal  Coral  36/05  vertical  calculation. 

Initial  wavefront  arrivals  and  peak  arrival  times  for  Coral  36/05  are  shown 
In  Figure  19.  The  peak  arrival  In  the  extreme  upper  portion  of  the  water  layer 
cannot  be  determined  because  the  spall  reflection  truncates  the  peak.  Note 
that  the  wavefront  velocity  in  coral  approaches  2000  m/sec  (slope  of  straight 
line),  in  agreement  with  the  assumed  sonic  velocity  of  1980  m/sec.  Tugboat  IJKL 
measured  arrivals  are  also  plotted  In  Figure  19.  Considerable  scatter  Is  evident, 
but  most  of  the  measurements  lie  near  the  calculated  curve.  Positions  for  three 
of  the  gages  listed  In  the  Tugboat  report^  lie  above  the  water  surface  level, 
evidently  an  error  (all  gages  were  underwater). 

Arrival  times  for  the  deep  water  layer  calculation  coral  36/36  are 
depicted  in  Figure  20.  Initial  and  peak  arrivals  in  the  coral  are  identical 
to  Figure  19.  The  velocity  in  the  water  layer  approaches  the  sonic  value 
of  1503  m/sec,  indicating  nearly-acoustic  behavior  of  the  wave. 
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The  SOC  Initial  shock  arrival  times  at  the  rock-water  Interface  and 
at  the  water  surface  for  all  vertical  calculations  are  compiled  In  Table  7a. 

These  values  may  be  used  to  aid  In  analysis  of  the  waveforms  presented  later. 

The  acoustic  velocity  pulse  transit  times  from  problem  center  to  problem 
boundary  (rock  surface  or  water  surface)  are  also  given  In  Table  7b.  The 
acoustic  transit  times  are  In  all  cases  a  fraction  of  a  millisecond  longer 
than  the  SOC-calculated  arrival  time  at  the.correspondlng  boundary. 

Peak  Parameters.  Appendix  A  presents  graphs  of  peak  particle 
velocities,  pressures,  radial  stresses  and  tangential  stresses  as  functions  of 
Initial  range,  for  both  horizontal  and  vertical  calculations  (Figures  A1-A47). 

These  are  the  maximum  values  of  each  variable  which  occur  In  a  given  zone  during 
the  course  of  the  calculation.  Thus,  a  peak  velocity  for  a  vertical  calculation 
may  represent  either  the  peak  at  the  shock  front,  or  the  spall  velocity  (If  It 
exceeds  the  Initial  peak).  Figures  A1-A15  show  the  peak  curves  for  the  horizontal 
(free-fleld)  calculations.  Note  that  the  two  coral  calculations  CORHRZ  and 
PS10HRZ,  for  the  "smoothed"  and  measured  equations  of  state,  are  almost  Identical 
(A1-A4).  Figures  A5,  A10,  and  A15  show  the  relationships  between  peak  pressure, 
radial  stress,  and  tangential  stress  for  each  of  the  rock  media.  Note  that  the 
three  curves  for  each  material  overlap  at  very  close  ranges,  where  the  solid  shows 
a  ductile  rapidly- falling  behavior,  and  the  development  of  shear  stress  Is 
severely  limited  by  the  ultimate  strength  of  the  material.  The  radial  and 
tangential  stresses  are  large  relative  to  the  permitted  shear  stress  and  are 
nearly  equal.  At  longer  ranges,  the  shear  stress  becomes  appreciable  (relative  to 
the  magnitudes  of  the  radial  and  tangential  components).  Thus,  the  tangential 
stress  drops  below  the  radial  value;  the  three  curves  of  radial  stress,  pressure, 
and  tangential  stress  approach  the  parallel  behavior  expected  for  an  elastic 
material  of  constant  Poisson's  ratio. 

Figures  A16-A19  plot  the  peak  velocities,  pressures,  and  stresses  for  all 
three  media  on  the  same  graph.  The  particle  velocities  (A16)  are  highest  for 
coral,  followed  by  basalt,  then  granite  (coral  Is  the  lowest  Impedance  medium, 
and  granite  Is  the  highest).  Attenuation  rates  of  velocity  and  stress  are  dis¬ 
tinctly  higher  for  granite  than  for  the  other  two  media,  particularly  at  ranges 
>  4  to  10  m.  This  Is  the  expected  effect,  since  shear  strength  Is  a  dominant 
dissipation  mechanism  In  these  saturated  (nonporous)  media.  The  three 
tangential  stress  curves  plotted  In  Figure  A19  all  show  precipitous  declines,  but 
at  somewhat  different  ranges:  granite  declines  between  3.5  and  4.8  m,  basalt 
between  12  m  and  16  m,  and  coral  near  25  m.  At  ranges  within  each  of  these  respective 
Intervals,  the  peak  shear  stresses  decrease  below  the  ultimate  strength  of  the 


TABLE  7(a).  Approximate  Initial  Arrival  Times  of  the  Pressure  Pulse  at  the  Rock 
Surface  (or  Rock-Water  Interface)  and  at  the  Water  Surface. 


Configuration 


Coral 

Coral 

Coral 

Coral 

Coral 

Coral 

Coral 

Coral 


30/00 

30/05 

30/10 

30/36 

36/00 

36/05 

36/10 

36/36 


Initial  Arrival 


Water 

(msec) 


Configuration 


Basalt  36/00 
Basalt  36/10 
Basalt  36/36 
Granite  36/00 
Granite  36/10 
Granite  36/36 


Initial  Arrival 


Rock 

(msec) 


Configuration 


Coral 

Coral 

Coral 

Coral 

Coral 

Coral 

Coral 

Coral 


30/00 

30/05 

30/10 

30/36 

36/00 

36/05 

36/10 

36/36 


Transit 

ms  pel 


4.05 

5.1 

6.1 
11.3 

4.97 

6.0 


Configuration 


Basalt  36/00 
Basalt  36/10 
Basalt  .’6/36 
Granite  36/00 
Granite  36/10 
Granite  36/36 


Transit 
msec 


Material  zones  adjacent  to  the  explosive  radius  first  become  significantly 
active  about  0.22  msec,  and  the  peak  shock  enters  the  medium  at  0.24  msec. 
Shock  front  velocities  somewhat  exceed  the  acoustic  velocity  within  the  first 
few  meters  near  the  explosive. 


material.  Thus,  the  shear  stress  Is  no  longer  severely  limited  by  material  strength, 
and  radial  stresses  larger  than  the  tensile  components  begin  to  build  up  (l.e.,  the 
material  responds  In  a  more  nearly  elastic  fashion).  This  change  occurs  at  close 
range  for  granite  because  of  Its  high  strength. 

Figures  A20  through  A27  are  comparative  plots  of  the  36/00  and  36/10 
vertical  problems  for  all  three  media.  In  the  case  of  velocities,  note  the 

thin  surface  spall  1«yer  near  the  Interface  (A20);  note  also  the  tensile  wave 
pressure  "cutoff"  effect  below  the  Interface  (A21).  For  the  36/10  water 
overburden  calculations  (A24),  a  spall  layer  occurs  immediately  below  the 
rock-water  Interface;  the  shock  wave  then  enters  the  water  layer,  attenuates 
further,  Is  reflected  at  the  water  surface,  and  creates  another  spall  layer 
near  the  water  Interface.  (The  medium  above  the  Interface  has  been  denoted 
as  "air"  In  the  diagrams;  the  upper  boundary  In  the  calculations  Is  actually 
a  free  surface,  with  no  constraint.)  Comparing  the  media,  we  find  that 
spall  velocities  are  highest  for  coral,  and  lowest  by  a  large  margin  for 
granite.  Comparing  the  36/00  and  36/10  configurations,  we  see  that  spall 
velocities  at  and  near  the  rock  Interface  are  decreased  by  the  presence  of 
the  water  layer  (A20,  A24).  Spall  velocities  In  the  water  layer  appear  to 
be  high  (relative  to  the  rock)  because  water  Is  a  low-impedance  medium.  The 
water  layer  pressures  and  velocities  also  decrease  In  the  order:  coral,  basalt, 
granite.  Finally,  an  abrupt  decline  In  the  water  velocities  very  near  the 
Interface  Is  apparent  for  granite  36/10  (Figure  A24).  The  velocity  continues 
to  attenuate  rapidly  In  the  water  layer  and  water  surface  spall  velocity  Is 
low.  Likewise,  the  pressure  for  granite  does  not  show  a  "rebound"  In  the 
water  layer  beyond  the  reflection  region  at  the  rock  Interface  (A25).  This 
damping  of  the  water  layer  pulse  for  granite  will  later  be  seen  to  be 
partly  due  to  a  calculatlonal  difficulty. 

Figures  A28  through  A31  present  velocities  and  pressures  for  the  coral 
vertical  configurations  at  Df  =  36  ft.  (The  coral  horizontal  calculations, 
shown  for  comparison,  appear  to  lie  above  the  "vertical"  curves  at  very  close 
ranges  (A28,  A30).  This  apparent  discrepancy  was  caused  by  dezoning  or  combina¬ 
tion  of  close-range  zones  at  late  times  during  the  calculation,  and  by  a  computer 
plotting  routine  which  draws  a  rough  "averaged"  fit  through  the  remaining  zones.) 

An  Interesting  effect  may  he  noted  In  Figures  A28-A29:  the  rock  Interface 
spall  velocities  (at  10.97  m)  are  Identical  for  coral  36/05,  36/10,  and  36/36, 
but  are  higher  for  the  free-surface  event  36/00.  It  Is  evident  that  water  layers 
of  all  depths  exert  a  "constant  Impedance"  effect  on  the  Interface  spall  reflec¬ 
tion:  Water  layer  depth  causes  no  apparent  changes  In  the  Initial  spall  velocity. 

Note  that  velocity  behavior  at  the  Interface  and  attenuation  In  the  water  layer 
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spall  "cutoff"  layer  for  coral  36/05  extends  almost  back  to  the  rock 
Interface  (A29).  Thus,  a  water  surface  reflection  might  increase  the  rock 
interface  peak  velocity  with  a  water  layer  shallower  than  5  ft.  A  rock 
spall  velocity  due  to  water  surface  reflection  which  exceeds  the  Initial 
spall  velocity  (due  to  the  rock-water  Interface  reflection)  was  not  encountered 
for  any  of  these  calculations,  and  Is  expected  to  occur  only  for  very  shallow 
water  layers. 

Additional  calculations  were  Implemented  to  test  the  effect  of  dezoning 
on  a  vertical  problem.  The  coral  36/05  event  was  run  to  20  msec  with  and 
without  automatic  dezoning  after  the  shock  peak  has  passed.*  Numerical  results 
In  the  two  cases  did  not  differ  significantly  (Figures  A29,  A31),  but  the 
dezoned  calculation  required  far  less  computer  time. 

Figures  A32  through  A35  depict  peak  velocities  and  pressures  for  the 

coral  Df  *  30  ft  calculations.  Once  more,  the  spall  velocities  at  the  rock 

Interface  (R  -  9.146  m)  are  equal  for  30/05,  30/10,  and  30/36,  but  are  considerably 

higher  for  30/00  (no  water).  Due  to  the  shallower  depth  of  burial,  all  of 

the  rock  and  water  spall  velocities  are  much  higher  than  the  corresponding 

values  for  the  *  36  ft  series, 
r 

Figures  A36  and  A37  show  velocity  and  pressure  for  the  basalt  events 
at  Df  *  36  ft.  Again,  the  basalt  "horizontal"  curve  deviates  because  of 
averaging  by  the  computer  program.  The  basalt  rock  Interface  spall  velocities 
are  high  for  36/00  (no  water),  lower  for  36/10  and  36/36.  However,  the 
apparent  difference  Is  less  than  for  coral  (A28).  This  trend  Is  expected, 
since  the  Impedance  mismatch  Is  greater  for  basalt  and  tensile  reflection 
at  the  rock-water  Interface  Is  more  efficient.  The  peak  values  in  the  water 
layer  are  smooth  and  well-behaved,  with  an  attenuation  rate  similar  to  the 
coral  calculations.  Similar  velocity  and  pressure  plots  for  granite  (Dr  =  36  ft) 
are  given  In  Figures  A38  and  A39.  In  this  case,  the  Impedance  mismatch  Is 
still  stronger,  and  the  36/00  calculation  deviates  only  slightly  from  the 
underwater  calculations  at  the  Interface.  Severe  attenuation  effects  In  the 
water  layer  are  evident. 

•Dezoning  refers  to  a  code  routine  which  gradually  combines  small  adjacent  zones 
In  the  Inner  portion  of  the  problem  far  behind  the  shock  front.  This 
technique  Increases  the  time  step  and  enables  the  calculation  to  run  more 
efficiently  without  seriously  degrading  resolution. 
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A  parameter  study  was  carried  out  to  test  the  effects  of  artificial 
viscosity  In  these  media.  Horizontal  calculations  In  coral  (smoothed 
equation  of  state)  and  basalt  were  run  with  artificial  viscosity  v  ■  0.1  and 
0.2.  The  effects  on  peak  velocity,  pressure,  and  stress  are  shown  In 
Figures  A40  through  A45.  Increasing  viscosity  reduces  all  of  the  peak  sho„k 
parameters;  the  reduction  Is  greatest  at  long  ranges.  Moderate  changes  are 
apparent  In  the  coral  calculation  beyond  5  m.  Changes  In  the  basalt 
calculation  are  extremely  small  Inside  11  m  (l.e.,  within  the  range  of 
Interest  for  vertical  calculations).  Greater  differences  appear  at  long 
ranges.  The  coral  calculation  was  observed  to  be  noise-free  at  v  »  0.1, 
while  a  larger  damping  factor  was  required  for  best  results  In  basalt. 
Accordingly,  a  viscosity  of  0.1  was  used  for  all  coral  calculations  presented 
In  this  study,  while  v  «  0.2  was  used  for  basalt  and  granite. 

In  conjunction  with  the  artificial  viscosity  test,  trials  were 
conducted  to  determine  the  Influence  of  zone  size  on  vertical  calculations. 
Results  for  typical  configurations  demonstrated  that  zone  size  has  virtually 
no  Influence  on  peak  shock  parameters,  waveforms,  or  late  time  Interactions. 
This  result  holds  true  because  of  the  extremely  fine  constant  rock  zone  size 
used  for  the  suite  of  vertical  calculations  (much  finer  than  zoning  normally 
used  for  horizontal  free-field  configurations^).  Shock  resolution  is 
essentially  not  limited  by  rock  zone  size  In  the  vertical  calculations. 

The  effect  of  water  zone  size  on  granite-to-water  shock  transmission 
was  mentioned  previously.  In  an  effort  to  reduce  'the  errors  encountered, 
granite  36/10  and  36/36  calculations  were  rerun  with  water  zone  sizes  reduced 
by  50%  (  fine  zone"  calculation;  not  Impedance  matched  at  granite-water 
Interface).  Results  are  depicted  in  Figures  A46  and  A47.  The  abrupt  drop 
In  water  pressures  and  velocities  Immediately  beyond  the  interface  is  reduced, 
and  the  water  spall  velocities  are  Increased.  However,  rapid  attenuation 
within  the  water  layer  persists  In  both  calculations.  The  anomalous  sudden 
damping  of  the  shock  near  the  Interface  Is  thus  caused  by  zone-size-dependent 
resolution  loss  In  the  coarsely  zoned  water  layer.  Continuing  rapid 
attenuation  In  the  water  Is  a  more  complex  phenomenon,  not  entirely  explained 
by  zone-size  effects.  This  topic  will  be  further  examined  In  Chapters  VIII 
and  X. 
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Interface  spallation.  Spall  launch  at  an  interface  is  caused  by 
Interaction  of  the  outgoing  compression  wave  with  the  reflected  tension 
wave.  At  a  free  surface  (used  to  represent  the  medium-air  Interface  In  the 
vertical  calculations),  perfect  reflection  occurs  and  the  peaks  of  the 
two  waves  coincide.  The  particle  velocity  Is  doubled.  As  the  tension  wave 
moves  downward.  It  Interacts  with  the  declining  tail  of  the  compression  wave, 
producing  lesser  reinforcement.  At  some  depth  within  the  medium,  dependent 
on  attenuation  and  the  rapidity  with  which  the  compressional  wave  .decays, 
spall  launch  velocity  will  be  lower  than  the  Initial  peak  particle  velocity 
at  the  same  location.  Spallation  Is  Illustrated  by  Figure  21,  which  shows 
water  layer  velocity  waveforms  In  the  zones  located  0.9  m  and  1.6  m  below 
the  coral  30/10  water  surface.  The  spall  velocity  after  the  tensile  arrival 
is  greater  than  the  Initial  peak  at  0.9  m.  However,  the  later  spall  arrival 
at  1.6  m  produces  a  lower  velocity  and  falls  to  match  the  Initial  peak.  This 
effect  Is  responsible  for  the  appearance  of  the  peak  plots  (l.e..  Figure  A20): 
The  spherical ly-dlvergent  compression  wave  determines  the  "peak"  velocity 
except  within  a  narrow  region  near  the  interface.  Since  the  peak  velocity 
at  the  wavefront  Is  a  transient  phenomenon,  It  does  not  control  dynamics 
within  the  mound.  Instead,  one  must  examine  the  spall  velocities  and  late¬ 
time  residual  velocities  long  after  the  peak  has  passed. 

Spall  Interactions  In  solid  media  are  frequently  somewhat  more  complex 
due  to  the  less  regular  waveforms  and  the  failure  and  tensile  strength  effects 
involved.  Figures  22-23  present  typical  pressure  waveforms  with  their  corres¬ 
ponding  velocity  waveforms  near  the  interface  for  coral  36/00  and  basalt  36/00. 
Freefal 1  begins  when  the  reflected  tensile  wave  returns  the  pressure  to  zero. 
Due  to  its  longer  rise  time,  flattened  peak,  and  slower  decay,  the  basalt 
pulse  produces  significant  spall  reinforcement  over  a  greater  depth  interval 
than  coral.  The  basalt  spall  velocity  exceeds  the  initial  peak  velocity  at 
R  =  9.01  m  (or  1.9  m  below  the  Interface),  while  the  coral  velocity  reinforce¬ 
ment  has  fallen  below  the  Initial  peak  at  R  =  9.49  m  (1.5  m).  By  the  same 
token,  spall  reinforcement  within  the  broad  basalt  spall  layer  is  somewhat 
less  Intense  and  well-defined  than  for  coral  (Figure  A20). 
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Figure  22a.  Pressure  waveforms  and  material  spallation 
near  interface  (coral  36/0° )• 


Pressure  waveforms  and  material  spallation 
near  interface  (basalt  36/OO). 


Figure  23b.  Pressure  waveforms  and  material  spallation 
near  interface  (basalt  36/OO). 


Waveforms.  Appendices  B  through  F  contain  velocity  and  pressure 
waveforms  at  selected  radii  for  all  of  the  SOC  calculations.  A  key  to 
the  nomenclature  precedes  Appendix  B  (p.  325).  Results  of  the  horizontal 
free-fleld  calculations  In  all  media  are  presented  as  Appendix  B,  followed 
by  the  waveforms  of  the  vertical  coral  calculations  (Dp  «  30  and  36  ft)  In 
Appendix  C,  the  vertical  basalt  diagrams  (Dr  -  36  ft)  In  Appendix  D,  and 
the  vertical  granite  calculations  (Dr  *  36  ft)  In  Appendix  E.  Finally,  examples 
of  the  waveforms  transmitted  Into  the  overlying  water  layer  for  vertical 
calculations  are  given  In  Appendix  F.  It  should  be  pointed  out  that  most 
of  the  plots  are  logarithmic;  therefore,  the  apparently  high  noise  level I 
appearing  In  some  pressure  waveforms  Are  actually  a  factor  of  50  to  100  lower 
than  the  peaks.  In  order  to  facilitate  Interpretation  of  the  vertical  data 
(Appendices  C-E),  Important  dynamic  Interactions  have  been  noted  on  the 
velocity  waveforms.  The  following  abbreviations  are  used:  R  >  rock  surface 
or  rock-water  Interface  spall  reflection  arrival;  W  >  water  surface  reflec¬ 
tion  arrival;  C  >  cavity  gas  acceleration  pulse  arrival.  In  addition, 
certain  other  significant  events  such  as  freefall  phases  or  deceleration 
during  recompaction  of  overlying  material  have  been  Indicated. 

Horizontal  waveforms  at  each  range  for  the  three  media  (Appendix  B) 
have  been  plotted  together.  These  calculations  Include  a  2.02  bar  overburden 
pressure  (Figure  B2),  which,  however,  has  little  Influence  on  the  close-range 
results.  Horizontal  calculations  for  the  two  different  coral  equations  of 
state  (CORHRZ  and  PS10HRZ)  have  been  compared  In  Figures  B28-30.  Trivial 
differences  are  noted  at  ranges  close  to  the  cavity.  The  artificial  viscosity 
calculations  for  v  *  0.1  and  v  >  0.2  In  coral  and  basalt  are  shown  In 
Figures  B31-B34.  Small  but  noticeable  differences  occur  at  ranges  far  from 
the  cavity;  for  lower  viscosity,  the  Initial  arrival  occurs  later,  pulse 
rise  time  Is  shorter,  the  peak  Is  higher,  and  post-peak  "noise"  Interactions 
are  shifted  to  slightly  later  times. 

In  some  Instances,  It  is  desirable  to  use  calculational  results  for  a 
comparison  of  cratering  phenomenology  In  different  media.  Velocity  waveforms 
for  the  coral,  basa’t,  and  granite  calculations  36/00,  36/10,  and  36/36  have 
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been  replotted  on  the  sane  scale  (Figures  24-29).  Locations  adjacent  to 
the  rock  interface  were  selected,  since  these  regions  are  crucial  to 
material  (dynamics  and  water  layer  effects  show  up  most  readily  near  the 
Interface.  Late-time  spall  velocities  for  the  36/00  configuration  (R  *  9.49 
and  10.6  m)  clearly  show  that  the  coral  and  basalt  events  are  comparable, 
while  the  granite  velocity  Is  low  by  more  than  a  factor  of  two  (Figures  24-28). 
The  coral  velocity  Is  slightly  higher  than  basalt  at  ranges  very  close  to 
the  interface  (10.6  m)  because  of  the  sharper  shock  front  In  coral.  Under¬ 
water  configurations  36/10  bear  out  this  trend  (Figures  26-27):  Time- 
dependent  dynamic  Interactions  for  the  media  are  distinctly  different, 
but  coral  and  basalt  ultimately  attain  similar  velocities.  Granite  remains 
much  lower.  Calculations  "36/36"  present  a  similar  picture,  except  that 
the  late-time  gas  acceleration  pulse  reaches  the  rock  surface  somewhat  sooner 
for  dense  high-velocity  granite  (Figures  28-29).  It  is  Interesting  to  note 
that  the  late-time  velocities  of  coral  and  basalt  near  the  Interface  for  the 
36/36  calculations  (deep  water  layer)  are  actually  greater  than  the  corres¬ 
ponding  granite  36/00  spall  velocities  (no  water  layer).  Late-time  velocities 
less  than  or  on  the  order  of  30  m/sec  clearly  Indicate  that  the  granite 
events  are  buried  well  below  optimum  cratering  depth.  It  Is  evident  that 
the  zeal  for  choosing  a  very  stiff,  high-strength  rock  medium 
was  carried  a  little  too  far.  The  parameters  and  calculatlonal  results  may 
be  unrepresentative  of  a  "real-world"  medium  near  the  ground  surface,  where 
flaws,  joints,  and  weathering  reduce  the  effective  strength  of  the  material. 

Calculated  waveforms  may  also  be  used  to  define  the  detailed  dynamic 
effects  of  water  overburden  on  events  In  a  given  medium.  Velocity  waveforms 
and  arrival  times  for  all  calculations  in  a  single  series  are  intercompared. 

In  some  cases,  the  free-fleld  waveforms  are  very  complex  and  must  be  lacluded 
in  the  comparison  as  well  (see,  for  example.  Figures  B7-B11,  granite).  A 
single  example  will  serve  to  Illustrate  the  technique.  Let  us  examine  water 
layer  effects  throughout  the  mound  for  coral  calculations  at  Dr  *  36  ft. 
Appropriate  data  are  compiled  in  Appendix  C  (Figures  C1-C54).  In  order  to 
facilitate  comparisons  of  this  sort,  all  of  the  velocity  waveforms  have  been 
replotted  on  an  expanded  linear  scale,  and  corresponding  waveforms  for  a 
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Lgure  2b.  Velocity  waverorms  lor  3t>/00  calculations  Ln  coral, 
basalt,  and  granite;  Range  R  =  ^.Ur‘  rn  (1.U8  m  l  rom 
medium-water  interlace). 
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glven  medium  have  been  placed  on  the  same  graph  where  possible  (Figures  30- 
55).  Thus,  the  velocity  effects  for  coral  36/00  and  coral  underwater 
calculations  (Df  8  36  ft)  may  be  Intercompared  In  Figures  30-39.  Wave¬ 
forms  are  shown  for  ranges  between  2.01  and  10.89  r.i.  (Coral  Df  *  30  ft 
waveforms  have  also  been  Included  In  Figures  30-3'>;  the  Interface  Is  located 
at  30  ft  *  9.14  m,  and  no  waveforms  are  given  beyond  this  range. )  A 
summary  of  the  comparative  water  layer  effects  for  Dr  8  36  ft  Is  presented 
below,  beginning  near  the  Interface  (R  8  10.97  m),  where  wate  overburden 
effects  are  most  evident,  and  moving  Inward. 

Range  8  10.89  m,  near  the  rock-water  Interface  or  rock  surface; 

Coral  36/00  (Figure  38a)  shows  the  classic  spall  waveform  very  close  to  a 
free  surface;  the  Initial  peak  Is  cut  off  by  the  tensile  reflection,  and 
the  velocity  at  the  free  surface  Is  approximately  doubled.  Material  enters 
the  f reefal 1  phase  (signified  by  "R")  almost  Inmedlately  (horizontal  line 
of  constant  velocity;  deceleration  by  gravity  Is  very  slight  over  the  brief 
time  Interval  shown).  Peak  velocity  due  to  reflection  at  the  rock-water 
Interface  for  coral  36/05,  36/10  and  36/36  Is  somewhat  lower  (Figure  38b-c), 
about  70  m/sec  Instead  of  100  m/sec.  This  decrease  occurs  because  the  rock- 
water  Interface  Is  not  a  free  surface  (tensile  reflection  Is  not  perfect). 
Note  that  all  three  underwater  calculations  have  the  same  Initial  peak  spall 
velocity  Independent  of  water  depth,  as  expected  for  a  constant  degree  of 
impedance  mismatch.  However,  the  overlying  water  layer  damps  motion  at  the 
Interface  after  the  initial  spall  reflection,  and  the  coral  material  doe  not 
Immediately  enter  freefall  spallation.  Subsequent  velocity  history  depends 
on  the  water  layer  depth;  freefall  does  not  begin  until  the  water  surface 
tensile  reflection  ("W")  returns  to  the  rock  layer  and  ends  the  damping 
phase.  The  reflection  transit  time  and  the  degree  of  damping  depend  on  the 
depth  of  the  water  layer.  Figure  38b  shows  that  the  spall  velocity  near  the 
Interface  Is  less  than  50  m/sec  for  coral  36/05,  and  lower  than  30  m/sec  for 
coral  36/10.  The  water  surface  reflection  for  coral  36/36  returns  to  the 
Interface  at  very  late  time  (Figure  38c).  Damping  Is  almost  complete,  and 
the  only  apparent  effect  of  the  spall  relief  Is  to  end  the  slight  late-time 
noise  and  Initiate  true  freefall.  The  damped  late  time  velocity  (or  the 


550 


(oas/m) 


Time  (msec ) 

Figure  30a.  Coral  R  =  2.01  m  (sample  calculations  at  very  close  range 
for  Coral  36/00,  3&/10,  30/00,  and  30/10  only). 
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Figure  36a.  Coral,  R  =  9*^9®  (compare  with  7.98  m,  coral  30  ft). 
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igure  37a.  Coral,  R  =  10.6  m  (con-are  with  8.96  n,  coral  30  ft). 
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Figure  37b.  Coral,  R  =  10.6  m  (compare  with  8.98  m,  coral  30  ft) 
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Figure  38a.  Coral,  R  =  10. 89  m  ( carfare  with  8.98  m,  coral  30  ft). 
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Figure  40a.  Basalt 
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Figure  42a.  Basalt,  R  =  6.01  m 
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Figure  42b.  Basalt 
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Figure  46a.  Basalt,  R  =  10.6 
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Figure  48b.  Granite,  R  *  3.01  m. 
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Figure  49®.  Granite,  R  =  4.01  m. 


Figure  49b.  Granite,  R  =  4.01  m 
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Figure  52a.  Granite,  R  =  9. 01  m. 
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Mb,  ir anitc,  fi  = 


or**  Ur 


Time  (msec) 


Figure  55b.  Granite,  R  =  10. 89  si 
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Figure  55c.  Granite,  R  =  10. 89  m 
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spall  velocity)  for  coral  36/36  Is  about  21  m/sec.  Shortly  after  the  36/36 
surface  reflection,  a  powerful  gas  acceleration  pulse  "C"  (moving  outward 
from  the  cavity)  arrives  and  Increases  the  velocity  to  about  30  m/sec.  Gas 
acceleration  pulse  transmission  near  the  Interface  Is  accentuated  for  deep 
water  layers,  due  to  the  Increased  velocity  damping  and  material  compaction 
near  the  Interface.  It  Is  evident  that  the  presence  of  water  (regardless 
of  depth)  Influences  the  Initial  reflection  at  the  rock  Interface;  the  depth 
of  the  water  layer  greatly  Influences  water  surface  reflection  spall  velocity 
and  late-time  history  very  close  to  the  rock-water  Interface. 

Range  «  10.6  m:  Motion  histories  shown  In  Figures  37a-c  are  similar 
to  the  trends  noted  above.  In  this  case,  the  shock  peak  arrival  and  the 
spall  cutoff  peak  are  slightly  separated  (greater  distance  from  the  Interface). 
Less  obviously,  the  damping  of  the  late-time  spall  velocities  for  the  under¬ 
water  events  (Figures  37b-c)  Is  less  severe  relative  to  the  coral  36/00  event 
(Figure  37a)  at  this  location  further  below  the  Interface. 

Range  *  9.49  m:  Figure  36a  reveals  that  this  location  Is  beneath  the 
narrow  region  of  Intensely  reinforced  spall  launch  near  the  Interface;  the 
rock  Interface  reflection  peak  "R"  Is  lower  than  the  Initial  peak  for 
both  the  coral  36/00  event  and  the  underwater  cases.  Note  that  the  coral 
36/00  "free-surface"  calculation  again  enters  freefall  upon  arrival  of  the 
rock  reflection,  whereas  the  underwater  events  are  damped  until  the  water 
surface  reflection  arrives.  However,  the  relative  damping  becomes  still  less 
outside  of  the  "Intense  spall"  region,  and  the  calculations  with  and  without 
water  overburden  are  beginning  to  converge. 

Range  *  9.01  m:  Figures  35a  and  35b  show  certain  changes  In  phenomenology 
well  away  from  the  Interface.  The  "rock"  spall  arrival  "R"  remains  weaker  for 
the  water-layer  calculations,  but  a  brief  freefall  phase  occurs  after  this 
arrival  (before  damping  due  to  recompaction  of  the  overlying  rock  and  water 
begins).  Additional  damping  effects  of  the  water  layer  become  much  less 
Important  relative  to  event  "36/00".  However,  the  velocity  decrease  for  36/36 
(Figure  35b)  Is  still  quite  severe  until  the  very  late- time  gas  acceleration 
pulse  arrives  and  returns  local  velocities  almost  to  the  level  which  prevailed 
before  the  damping  occurred. 
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Range  ■  8.01  m:  Further  modi fl cations  In  relative  water  layer  effects 
are  visible  at  this  range.  The  coral  36/00  and  36/05  waveforms  (Figure  34a) 
are  quite  similar,  undergoing  freefall  spallation  after  the  "R"  reflection 
arrives.  No  water  layer  damping  effects  act  upon  coral  36/05  within  the  time- 
scale  of  the  problem.  Coral  36/10  also  reveals  a  freefall  phase,  which 
terminates  with  recompactlon  of  overlying  material  at  9  msec.  However,  the 
velocity  Is  decreased  by  only  about  10  m/sec.  Approximate  freefall  motion 
develops  at  later  times,  although  slight  noise  pulses  occur  due  to  Inter¬ 
actions  with  adjacent  zones.  Figure  34b  shows  a  more  significant  damping  for 
event  36/36  (velocity  decreased  from  about  50  to  35  m/sec).  In  spite  of  the 
damping  and  deceleration,  no  water  surface  reflection  (late-time  freefall 
phase)  Is  apparent  In  Figure  34b;  the  cavity  gas  acceleration  pulse  (outward- 
moving)  actually  precedes  the  water  surface  reflection  (Inward-moving),  and 
completely  overwhelms  any  further  dynamic  effects  of  the  second  or  water  surface 
reflection.  (NOTE:  Figures  34a-b  and  subsequent  figures  also  contain  waveforms 
for  coral  "30  ft"  calculations;  these  are  not  directly  coeiparable  with  coral 
"36  ft"  at  the  same  range,  due  to  the  difference  In  Interface  locations.) 

Range  ■  6.01  m:  At  this  range,  water  surface  reflection  effects  are 
becoming  negligible  In  all  cases,  but  recompactlon  still  Influences  material 
dynamics.  Note  that  the  waveforms  and  "rock  Interface"  reflection  spall  are 
closely  similar  for  coral  36/00,  36/05,  36/10,  and  36/36  (Figures  33a-b). 
However,  the  cavity  gas  acceleration  pulse  arrives  fractionally  sooner  for 
the  underwater  events,  due  to  the  slightly  lower  spall  velocities  of  material 
moving  away  from  the  cavity.  The  coral  36/36  calculation  (Figure  33b)  reveals 
a  significant  late-time  deceleration  due  to  recompactlon  of  overlying  layers 
after  the  gas  acceleration  pulse  arrives.  The  velocity  decreases  back  to 
approximately  Its  Initial  "freefall"  value.  Thus,  the  gas  acceleration  pulse 
has  only  a  temporary  Influence  on  the  waveforms  and  does  not  greatly  modify 
overall  mound  dynamics  for  these  calculations.  Similar  effects  have  been 
noted  at  Intermediate  or  close  ranges  In  the  other  media  as  well. 

Range  3  5.01  m:  Here  the  waveforms  are  closely  comparable  for  all 
calculations  (Figures  32a-b).  Again,  spall  velocities  are  a  little  lower 
and  cavity  gas  acceleration  begins  slightly  sooner  for  the  underwater  events. 
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Thc  gas -accel era t 1 on - 1 nduced  velocities  are  low  by  about  the  sane  amount 
as  the  spall  velocities;  this  Is  obviously  a  residual  effect  of  the 
Imperfect  reflection  at  the  rock-water  1nterface(  and  persists  throughout 
the  calculation.  The  velocity  difference  Involved  (at  R  *  5.01  m)  Is 
negligible.  Once  more,  the  coral  36/36  calculation,  which  was  carried  to 
late  time,  shows  a  recompaction  and  deceleration  phase  which  roughly  cancels 
the  Increase  due  to  the  gas  acceleration  pulse  (Figure  32b). 

Ranges  4.01  and  3.01  m:  Interactions  parallel  the  above  description, 
but  the  freefall  phase  between  the  rock  surface  reflection  "R"  and  the  gas 
acceleration  pulse  "C"  gets  progressively  shorter  at  ranges  close  to  the 
cavity  (Figures  31a-b  and  30b-c).  Finally,  near  the  cavity  boundary,  the  "R"  and 
"C"  events  coincide  and  there  Is  no  true  "spall"  phase  (Figure  30a,  Range  - 
2.01  m).  The  gas-acceleration  pulse  Is,  of  course.  Initiated  when  the  first 
spall  reflection  returns  to  the  cavity  Interface.  The  redaction  phase, 
which  returns  late-time  velocities  to  approximately  the  "R"  spall  value.  Is 
likewise  evident  (Figures  31b,  30c). 

Some  Interesting  aspects  of  the  late-time  deceleration  are  also  evident 
In  the  shallower  coral  "30  ft"  calculations.  In  this  case,  successive  Inter¬ 
actions  occur  earlier  due  to  the  shallower  rock  layer  depth,  and  at  shorter 
Intervals  due  to  the  higher  velocities  Involved.  The  deceleration  for  coral 
30/10  Is  visible  In  Figures  32a,  31a,  and  30b;  velocity  returns  to  the  "spall 
reflection"  value.  The  deceleration  for  coral  30/36  appears  In  Figures  33b, 

32b,  31b,  and  30c.  The  relative  velocity  decrease  Is  slightly  greater  than 
for  36/36  (higher  material  velocities  are  Involved  and  more  violent  Interactions 
can  occur).  The  velocity  decrease  for  30/36  occurs  relatively  soon  after 
the  gas  acceleration  pulse  arrival,  almost  Immediately  at  R  =  5.01-6.01  m 
(Figure  33b,  32b).  The  pulse  develops  a  rather  complex  double  structure 
Involving  acceleration  and  subsequent  deceleration  at  successive  ranges 
moving  away  from  the  cavity  (Figure  32b,  33b,  34b,  35e).  At  locations  near 
the  rock-water  Interface,  part  of  the  energy  transferred  by  gas  acceleration 
remains  with  the  material,  and  velocity  does  not  decrease  all  the  way  back 
to  the  "residual"  value  before  gas  acceleration  (Figures  34b,  35e;  ranges  of 
7.98  and  8.98  m).  Similar  but  less  well-defined  effects  are  detected  near 
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the  coral  36/36  Interface  (Figures  35b,  36b,  37c,  etc.;  ranges  >  9  an). 
Unfortunately,  the  Interactions  occur  later  for  coral  36/36  and  the  cal¬ 
culations  were  not  run  long  enough  to  fully  define  very  late  time  effects. 

The  overall  dynamic  Interactions  for  coral  Df  «  30  ft  calculations 
generally  parallel  those  for  Dr  -  36  ft,  but  are  shifted  closer  to  the 
cavity  In  the  Intermedlate-to-long  range  region  (Interface  at  9.14  m)  and 
are  compressed  In  time.  The  dynamic  regimes  for  basalt  (Figures  39-47)  and 
granite  (Figures  48-55)  are  likewise  fundamentally  similar,  the  chief 
differences  being  In  the  shock  wavefront  contours,  which  are  more  rounded 
for  the  basalt  (deep  and  Ill-defined  surface  spall  layer), but  sharp 
and  "noisy"  for  granite  (complex  stress  relief  and  spall  Interactions  near 
an  Interface).  The  granite  36/36  calculations  actually  show  two  gas 
acceleration  phases.  Initiated  by  the  first  tensile  wave  reflection  from  the 
rock-water  Interface  and  by  the  subsequent  reflection  of  the  first  gas 
acceleration  pulse  from  the  rock-water  Interface.  This  Interaction  Is 
apparent  as  two  successive  f reef all  phases  near  the  cavity,  each  followed 
by  an  acceleration  pulse  and  subsequent  deceleration  (Figures  49b,  50b, 

51b).  These  signals  proceed  toward  the  outer  regions  of  the  problem  as 
two  distinct  acceleration  pulses  separated  by  about  13  msec  (Figures  52b, 

53b).  The  multiple-pulse  structure  occurs  because  of  the  very  high  assumed 
sonic  velocity  In  granite  and  the  low  spall  velocities  In  these  calculations 
(material  does  not  become  greatly  rarefied  by  spallation;  recompaction 
and  transmission  of  later  pulses  are  very  efficient).  Reflection  effects 
would  become  less  severe  for  an  "opt1mum"depth  calculation.  Late-time 

velocities  near  the  Inner  regions  of  the  mound  tend  to  return  almost  to  the 
Initial  "R"  spall  level,  with  a  slight  Increase  (Figure  50b),  In  spite  of 
the  second  pulse.  Second  or  subsequent  cavity  pulses  In  a  one-dimensional 
calculation  are  difficult  to  Interpret  physically  because  of  convergence 
effects.  Generally,  multiple  reflections  are  of  little  Interest  for 
shallow  cratering  calculations;  the  spall  launch  rarefies  mound  material 
(decreasing  multiple  pulse  transmission),  and  velocity  fields  are  already 
established  before  two  or  three  source- to-surf ace  transit  times  have  elapsed. 
Host  of  the  mound  material  enters  ballistic  freefall  and  strong  "hydrodynamic" 
Interactions  cease  to  be  Important. 
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In  summary ,  cratering  dynamics  for  underwater  events  are  controlled 
by  the  initial  (rock-water  Interface)  spall  reflection,  recompactlon  due 
to  the  water  burden,  first  gas  acceleration  pulse  (Initiated  by  the  "Inter¬ 
face”  spall  wave),  and  by  water  surface  spall  reflection.  The  water  layer 
effects  of  recompactlon  and  particularly  the  water  surface  spall  reflection 
are  most  significant  near  the  rock-water  Interface,  becoming  less  Intense 
within  the  rock  mound.  Deceleration  of  the  rock  Interface  due  to  recompactlon 
Increases  with  Increasing  water  l*yer  depth,  up  to  a  limiting  value  which 
corresponds  with  late-time  residual  mound  velocity.  The  limit  Is  apparently 
approached  In  these  calculations  for  Dw/Dr  *  1.0  (l.e.,  30/36  and  36/36  events). 
Presence  of  a  water  l*yer  causes  certain  other  dynamic  modifications  as  well. 

The  gas  acceleration  pulse  transmits  more  efficiently  and  arrives  sooner  In 
the  upper  regions  of  the  mound,  due  to  recompactlon  of  the  mound  material. 

This  effect  may  partially  compensate  for  the  water  layer  dancing  In  some 
cases.  The  recompactlon  phase  near  a  rock-water  Interface  ends  and  true 
f reefal 1  begins  when  the  water  surface  tensile  reflection  arrives  (except  for 
the  deep  water  layer  events  Dw/Dr  =  1.0;  In  these  cases,  the  material  has 
already  been  decelerated  to  residual  velocity  before  the  reflection  arrives, 
and  the  strong  gas  acceleration  pulse  may  even  precede  the  surface  reflection). 
Water  surface  reflection  Is  thus  Important  near  the  Interface  (except  for 
"deep  layer"  events);  the  reflection  becomes  progressively  less  1n*>ortant 
In  the  deeper  portions  of  the  rock  mound,  and  does  not  cause  an  Independent 
gas  acceleration  phase(i.e.,  gas  acceleration  Is  associated  with  the  first 
Interface  reflection  rather  than  the  water  surface  reflection).  Likewise, 
multiple  reflections  In  the  water  layer  and  multiple  gas  acceleration  pulses 
exert  little  Influence  for  optimum  depth  cratering  events  with  high  spall 
velocities.  This  conclusion  may  not  apply  to  very  deeply  burled  low-velocity 
events  not  characterized  by  strong  spallation. 

Water-layer  Waveforms.  Selected  velocity  and  pressure  waveforms 
within  the  overlying  water  layer  have  been  plotted  In  Appendix  F  (coral 
configurations.  Figures  F1-F3;  basalt,  Figures  F4-F6;  granite,  Figures  F7-F8). 
Waveforms  within  the  rock  layer  adjacent  to  the  Interface  have  been  Included 
In  Figures  F1-F6,  to  demonstrate  velocity  and  pressure  continuity  across  the 
interface.  The  beginning  of  the  f reefal 1  phase  (water  surface  tensile 
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reflection,  propagating  downward  through  the  water  layer)  Is  Indicated  by 
an  arrow  In  the  velocity  waveforms.  Water  spall  velocities  are  compared 
with  detailed  "acoustic1*  calculations  and  with  experimental  results  In 
Chapters  VIII  and  IX.  The  gas  acceleration  pulse  arrival  near  the  Interface 
Is  also  visible  at  late  time  In  Figures  F3,  F6,  and  F8  (coral  36/36,  basalt 
36/36  and  granite  36/36),  but  does  not  have  time  to  recompact  and  propagate 
through  most  of  the  rarefied  (spalled)  water  layer. 

Water  waveforms  for  the  "coarse"  and  "fine"  water  zone  granite 
calculations  are  presented  In  Figures  F7  and  F8.  The  "fine"  calculations 
show  later  arrival  times,  shorter  rise  time,  higher  "peak"  values,  slight 
differences  In  noise  behavior  and  late-time  arrivals,  and  higher  spall 
velocities.  Water  layer  effects  for  granite  will  be  c'alyzed  In  Chapter  x. 

VII.  HOUND  DYNAMICS 

Peak  velocities  correspond  with  long-term  spall  velocities  only  In 
the  narrow  "cutoff"  region  near  a  free  surface,  where  tne  spall  arrival  Is 
higher  than  the  Initial  peak  and  no  subsequent  recompactlon  effects  occur. 

All  other  "peaks"  throughout  the  problem  history  are  transient  values  at 
the  shock  front  or  the  temporary  "undamped"  values  before  recompactlon. 

In  order  to  analyze  long-term  mound  dynamics,  we  should  Isolate  the  late¬ 
time  residual  velocity  and  spall  velocity  data  from  the  waveforms,  or, 
better  still,  determine  the  "averaged"  dynamic  parameters  within  appropriate 
problem  regions  (l.e.,  cavity  pressure,  overall  mound  momentum  and  kinetic 
energy,  etc.).  This  section  summarizes  dynamic  parameters  for  the  eighteen 
vertical  calculations.  The  overall  mound  dynamic  trends  may  also  be  correlated 
with  features  appearing  In  the  flner-resolutlon  Individual  waveforms 
(Figures  30-55,  Appendices  B-E). 

Residual  velocities.  Calculated  spall  velocities  and  late-time 
residual  velocities  within  the  rock  mound  are  shown  In  Figures  56  (coral, 

Df  =  36  ft),  57  (coral,  Dr  =  30  ft),  58  (basalt,  Dr  =  36  ft),  and  59  (granite, 
Dr  *  36  ft).  These  velocities  were  arbitrarily  selected  at  times  near 
20  msec.  Thus,  the  gas  acceleration  phase  and  the  very  late-time  recompactlon 


Interface,  36 


Residual  velocities,  Tugboat  coral 


Figure  58.  Residual  velocities,  Kaalualu  weak  basalt 
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after  gas  acceleration  are  excluded  for  the  outer  regions  of  some  of  the 
curves.  The  errors  Incurred  are  not  Urge,  since  velocity  changes 
after  20  msec  are  relatively  minor,  and,  as  we  have  seen,  the  gas  accelera¬ 
tion  velocity  pulses  are  oscillatory  and  tend  to  be  self-damping  (l.e.,  the 

gas  acceleration  pulse  Just  restores  velocity  lost  during  the  previous 
recompaction). 

The  granite  calculations  form  a  special  case,  due  to  the  "deeper- 
than -optimum"  burial,  low  spall  velocities,  and  high  sonic  velocity.  The 
gas  acceleration  pulse  Is  Intense,  and  propagates  through  about  8  m 
of  the  mound  (for  granite  36/00  and  36/10)  or  all  of  the  mound  (granite 
36/36)  prior  to  20  msec.  Likewise,  severe  recompaction  effects  After  the 
first  gas  acceleration  phase  and  other  late-time  Interactions  are  partially 
Included  for  granite  36/36  at  ranges  closer  than  6  m.  These  facts  explain 
the  somewhat  peculiar  behavior  shown  In  Figure  59.  Granite  36/36  residual 
velocities  are  lower  than  the  other  curves  at  close  ranges  because  of  the 
severe  deceleration  Included  Immediately  after  the  gas  acceleration  phase 
(cf.  Figures  51a  and  51b).  Granite  36/00  and  36/10  lie  close  together, 
with  36/10  slightly  lower  due  to  Imperfect  spall  reflection  at  the  Interface 
and  correspondingly  decreased  late-time  velocities.  At  ranges  beyond  about 
9  m,  the  situation  Is  reversed:  the  granite  36/00  velocities  are  quite  high 
(undamped  spallation  near  a  free  surface),  while  the  granite  36/10  veloci¬ 
ties  are  substantially  decreased  (water  layer  damping,  relieved  by  the  water 
surface  reflection  at  later  times).  The  gas  acceleration  pulse  Is  not 
Included  In  these  velocities.  The  Interactions  are  displayed  In  Figures 
52a,  53a,  54a,  and  55a.  The  granite  36/36  velocities,  on  the  other  hand, 
are  very  strongly  damped  by  the  water  layer  but  also  Include  a  strong  gas 
acceleration  pulse  (Figures  52b,  53b,  54b,  and  55b).  The  gas  acceleration 
pulse  reaches  the  outer  problem  regions  due  to  strong  recompaction  of  the  rock. 
The  initial  gas  acceleration  peak  compensates  for  the  strong  water  damping, 
and  granite  36/36  residual  velocities  (with  gas  acceleration)  are  equal 
to  the  granite  36/10  residual  velocities  (no  gas  acceleration)  beyond  9  m. 

This  comparison  may  be  a  little  misleading,  since  damping  of  the  granite 
36/36  gas  acceleration  velocities  In  the  outer  rock  .rones  will  continue  at 
very  late  times.  Thus,  we  see  that  very  late  t1m»  interactions  and  multiple 
pulse  reflections  may  become  significant  for  deeply  burled  events  with  very 


-157- 


low  spall  velocities.  These  considerations  may  be  Ignored  for  spall - 
dominated  cratering  events  In  saturated  media  (l.e.,  coral  and  basalt 
calculations  with  weak  late-time  Interactions). 

Discounting  the  minor  granite  difficulties*  residual  velocity  curves 
In  Figures  56  through  59  present  a  consistent  picture  of  underwater  mound 
dynamics.  The  presence  of  a  water  layer  decreases  the  residual  velocities 
dramatically  within  2  m  of  the  Interface.  Most  of  the  decrease  occurs 
for  even  the  shallowest  water  layers  (see  30/05*  36/05,  and  36/10  calcula¬ 
tions),  with  slight  additional  decreases  for  water  depths  up  to  36  ft 
The  dramatic  decrease  near  the  Interface  Is  caused  primarily  by  Imperfect 
spall  reflection  for  the  water  layer  calculations  (not  layer-depth  dependent), 
and  secondarily  by  water  layer  recompactlon  (depends  on  layer  depth  and 
elapsed  time).  As  expected,  the  velocity  decrease  Is  most  severe  for  coral 
(slight  Impedance  mismatch,  low  reflection  factor  at  the  Interface),  less 
severe  for  basalt,  and  least  severe  for  granite  (strong  Impedance  mismatch, 
almost  perfect  "free  surface"  reflection).  Water  layer  effects  decrease 
further  below  the  Interface,  and  the  velocities  for  both  free  surface  and 
underwater  calculations  tend  to  converge  on  high  values  (s  100-140  m/sec, 
coral  and  basalt).  The  exception  to  this  rule  occurs  for  the  coral  30/10 
and  30/36  events.  Figure  57.  Their  close-range  velocities  are  somewhat 
damped  relative  to  coral  30/00  and  30/05  because  "very  late  time"  deceler¬ 
ation  (after  the  cavity  gas  acceleration  pulse)has  already  occurred 
at  20  msec  (see  Figures  32a-b,  and  31 a-b) .  Similar  deceleration  and  recompac¬ 
tlon  of  the  mound  for  coral  30/00  and  30/05  will  proceed  at  slightly  later 
times.  The  close-range  differences  are  largely  apparent  rather  than  real. 

Note  that  the  late-time  recompactlon  begins  relatively  early  for  Dr  *=  30  ft 
events  because  of  the  shallow  scaled  burial,  very  high  spall  velocities, 
and  correspondingly  rapid  Interactions.  The  deep  mound  velocities  remain 
extremely  high  even  after  recompactlon  has  begun,  as  expected  for  shallow 
spall -dominated  cratering  events.  Late-time  mound  recompactlon  effects  close 
to  the  cavity  are  not  expected  to  modify  any  conclusions  presented  here  for 
"optimum  depth"  spall  cratering  events. 
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Cavlty  growth  and  close-range  Interactions.  The  time-dependent 
cavity  radii  "Rcay"  and  cavity  pressures  "Pcav"  are  plotted  In  Figures  60-62 
and  Figures  63-65  (a-b),  respectively.  The  corresponding  horizontal  calcula¬ 
tion  curves  are  Included  for  comparison  (free-fleld  cavity  growth  In  an 
Infinite  medium;  no  relief  pulse  reflected  from  an  Interface  or  free  surface). 
Initial  cavity  radius  Is  1.13  m  («  explosive  radius),  and  growth  begins 
about  0.3  msec.  The  cavity  pressures  are  plotted  beginning  at  0.1  msec, 
while  the  detonation  Is  still  In  progress.  Pressures  Increase  until  about 
0.26  msec,  when  the  peak  shock  enters  the  surrounding  medium,  then  drop  off 
as  the  cavity  expands.  All  pressures  are  average  values  weighted  by  zone 
volume  (zone  volumes  are  smallest  hear  the  center  of  the  spherical 
calculation;  volume  weighting  of  pressures  throughout  the  cavity  compensates 
for  this  effect).  Each  set  of  pressure  curves  has  been  plotted  at  compressed 
scale  (10  ^-10  3  W>ar,  Figures  53-65a),  and  also  on  an  expanded  scale  to  better 
reveal  late-time  effects  (10"3-l(f5  Mbar,  Figures  63-65b). 

Late-time  cavity  expansion  and  cavity  pressures  are  descriptive 
of  the  gas  acceleration  coupling  and  of  dynamic  effects  Immediately  adjacent 
to  the  cavity.  Cavity  growth  follows  "free-fleld"  expansion  (Infinite,  medium) 
until  the  spall  reflection  returns  from  the  nearest  free  surface  or  Inter¬ 
face.  Growth  Is  then  accelerated  toward  the  free  surface  and  the  gas  acceler¬ 
ation  phase  begins.  Deviation  from  free-fleld  behavior  begins  soonest 
for  a  shallow  event  (nearby  Interface);  the  cavity  expansion  Is  also  most 
rapid  and  pressure  dropoff  Is  fastest  for  a  shallow  burial  depth.  As  burial 
depth  Increases,  accelerated  cavity  expansion  begins  later  and  occurs  more 
slowly.  In  the  limit  (horizontal  or  free-fleld  event,  no  Interface),  cavity 
growth  proceeds  until  the  pressure  approximately  equals  the  lithostatlc  head 
plus  restraining  material  strength  of  the  adjacent  rock.  The  cavity  then 
begins  to  stabilize,  although  some  overshoot  and  "ringing"  effects  are 
expected  at  late  time.  Cavity  growth  rates  and  pressures  thus  reveal  the 
effects  of  Interface  reflections,  overburden  (including  water  layer  depths), 
and  material  strength  on  gas  acceleration  and  late-time  Interactions. 
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Figure  60.  Cavity  radius  as  a  function  of  time  (coral  configurations) 
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Figure  6l.  Cavity  radius  as  a  function  of  time  (basalt  configurations). 
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Figure  62.  Cavity  radius  as  a  function  of  time  (granite  configurations 
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Figure  63«.  Cavity  pressure  as  a  function  of  time  (coral  configurations). 
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Figure  65a. 


Cavity  pressure  as  a  function  of  time  (granite  configurations). 
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Cavity  growth  for  coral  Is  depicted  In  Figure  60,  while  the 
pressures  are  plotted  In  Figures  63a-b.  Cavity  expansion  for  the 
horizontal  calculation  ceases  after  about  12  msec  (cavity  pressure 

5  10  Mbar  or  0.1  kbar  •=  ultimate  strength  of  coral).  The  vertical  calcu¬ 
lations  begin  to  deviate  after  about  7  to  8  msec,  when  the  first  Inter¬ 
face  reflection  Impinges  on  the  cavity  for  the  "30  ft"  configurations. 

(The  "horizontal"  curves  differ  very  slightly  before  7  msec  because  of  the 
coarse  zoning  used.)  Cavity  growth  and  pressure  decline  are  most  rapid 
for  the  shallow  "Dr  *  30  ft"  calculations,  and  are  somewhat  less  rapid 
for  the  "Dr  *  36  ft"  events.  Water  layer  effects  are  less  Important 
than  rock  burial  depth  "Dr".  The  presence  of  a  water  layer  damps  cavity 
growth  and  Increases  cavity  pressure  somewhat.  However,  at  times  earlier 
than  14  msec,  the  curves  for  coral  30/05,  30/10  and  30/36  are  virtually 
Identical,  as  are  the  curves  for  36/05,  36/10,  and  36/36.  This  fact 
Indicates  that  the  presence  of  the  water  layer  Influences  early- time  cavity 
growth  and  gas  acceleration,  but  the  depth  of  the  water  layer  Is  unimportant. 


It  is  evident  that  imperfect  spall  reflection  from  the  rock-water  Interface 
causes  the  decreased  cavity  growth  rates  relative  to  the  "free-surface" 
events  (coral  30/00  or  36/00).  At  times  later  thati  14  to  17  msec,  water 
layer  depth  dependence  becomes  apparent  In  the  coral  30/10  and  30/36 
calculations.  The  thick  water  layers  damp  cavity  growth  and  Increase 
pressure.  This  effect  is  caused  by  the  close-range  damping  discussed  In 
the  previous  section.  Due  to  the  shallow  burial  depths  and  high  spall 
velocities  for  the  "30  ft"  calculations,  mound  recompaction  occurs  at 
relatively  early  times.  Thus,  the  gas  acceleration  phase  near  the  cavity 
Is  concluded  and  cavity  growth  Is  damped  before  20  msec  (termination  of 
plotted  curves).  This  effect  is  visible  In  Figures  30  a-b-cj  recompaction 
and  deceleration  begin  at  about  14  msec  (30/36)or  17  msec  (30/10),  the 
same  times  as  the  deviations  visible  In  Figures  60  and  63.  As  noted  pre¬ 
viously,  mound  recompaction  near  the  cavity  will  also  occur  (somewhat  later) 
for  the  corai  30/00  and  30/05  events,  thereby  partly  compensating  for  the 
apparent  differences  In  cavity  growth.  The  calculations  were  not  carried 
long  enough  to  analyze  very  late  time  cavity  damping  after  the  gas  acceleration 
phase  for  coral  30/00  and  30/05.  To  place  the  "water  layer"  effects  In 
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proper  perspective,  we  note  that  the  coral  30/36  cavity  radii  decline  and 
the  pressures  Increase  almost  to  the  levels  of  the  36/00  calculation  at 
very  late  times.  However,  the  coral  30/36  event  actually  has  37*  more 
overburden  weight  than  coral  36/00.  In  addition,  the  curves  should  separate 
again  at  later  times,  as  cavity  damping  begins  for  the  36/00  event.  Rock 
burial  depth  and  rock  overburden  are  more  important  than  water  layer  over¬ 
burden  In  determining  cavity  growth,  and  the  first  Interface  reflection  Is 
the  principal  factor  controlling  cavity  Interactions  In  these  calculations. 

Cavity  radii  and  pressures  for  basalt  are  presented  In  Figures  61 
and  64  a-b.  Horizontal  cavity  growth  decelerates  at  about  6  msec  (?c 
2  0-5  kbar  :  ultimate  strength).  The  vertical  calculations  36/00,  36/10, 
and  36/36  are  almost  Identical  up  to  20  msec.  In  accord  with  the  similar 
close-range  waveforms  (Figure  29a-b).  Some  cavity  damping  would  be 
expected  for  basalt  36/36  after  20  msec  (Figures  39b,  gas  acceleration 
phase  ends  and  recompaction  Is  apparent  near  25  msec),  and  at  later  times 
for  the  other  events.  The  basalt  36/10  and  36/36  growth  rates  are  marginally 
slower  and  pressures  slightly  higher  than  for  36/00,  due  to  Inperfect 
reflection  at  the  Interface. 

Granite  cavity  expansion  and  pressures  are  shown  In  Figures  62  and 
65a-b.  Horizontal  cavity  growth  ends  abruptly  at  5  msec,  and  the  pressure 
Is  almost  Immediately  stabilized  (Pcav  a  4  kbar  :  ultimate  strength).  This 
behavior  Is  expected  for  a  very  high-strength  rock.  Cavity  growth  Is  most 
rapid  for  granite  36/00,  somewhat  slower  for  36/10  and  36/36  (spall  reflection 
slightly  weaker),  lhe  cavity  growth  rate  Is  somewhat  damped  and  pressures 
are  Increased  for  granite  36/36  after  13  msec  (compare  with  36/10  curves). 

This  effect  is  again  attributed  to  cavity  damping  after  the  first  gas 
acceleration  phase;  Figure  48b  shows  rather  severe  recompaction  effects 
near  the  cavity  for  36/36  during  the  later  stages  of  the  calculation. 

To  conclude,  cavity  growth  and  gas  acceleration  effects  are  dominated 
uy  the  first  interface  reflection  arrival  at  the  cavity.  Cavity  growth 
is  most  rapid  for  a  strong  "free-surface"  spall  reflection,  and  Is  somewhat 
retarded  for  a  weaker  rock-water  Interface  reflection.  The  relative  degree 
of  retardation  Is  greatest  for  a  small  impedance  mismatch,  which  causes 


-170- 


Inefflclent  reflection  (l.e.,  coral-water  Interface).  Secondary  effects  of 
water  layer  depth  appear  only  at  late  tines,  after  the  first  gas  acceleration 
phase  Is  complete.  Recompaction  and  damping  of  cavity  growth  occur  earlier 
for  deep  water  layers.  This  deceleration  may  be  partly  compensated  at  later 
times,  when  cavity  damping  must  occur  even  for  the  "shallow  water  layer" 
or  "free  surface"  events.  In  any  case,  close- range  mound  recompaction  after 
the  gas  acceleration  phase  has  little  Influence  on  optimum  depth  cratering. 

It  Is  therefore  expected  that  the  presence  of  a  water  layer  will  exert 
primary  control  over  cratering  dynamics,  and  layer  thickness  will  be  a 
second-order  effect>  As  a  corollary  of  these  conclusions,  water  layer  depth 
for  moderate- to- thick  water  layers  has  far  less  effect  on  cavity  growth  and 
gas  acceleration  than  an  equivalent  mass  of  rock  overburden.  A  rather  abrupt 
change  In  dynamics  Is  expected  for  shallow  water  layers,  with  much  less  sensi¬ 
tivity  to  Increasing  overburden  at  greater  water  depths.  No  direct  parallel 
should  be  drawn  between  "rock  overburden"  and  water  layer  thickness. 

Hound  kinetic  energies.  Cratering  dynamic  effects  are  perhaps  best 
described  by  the  mound  kinetic  energy. ^  The  kinetic  energy  Is  a  velocity- 
sensitive  and,  therefore,  burial-depth  sensitive  parameter.  Kinetic  energy 
Is  Initially  transferred  from  the  explosive  to  surrounding  material  by  the 
detonation  shock  front  and  by  subsequent  spherical  expansion  of  the  cav- 

1ty*  After  the  Initial  shock  wave  Is  reflected  at  an  Interface  or  free 
surface,  spall  launch  of  the  mound  material  begins.  The  kinetic  energy  then 
Increases  as  spallation  spreads  throughout  the  mound.  This  Increase  occurs 
soonest  and  Is  greatest  for  shallowly-burled  events.  As  the  depth  of  burial 
Is  Increased,  spall  becomes  a  less  efficient  energy  transfer  mechanism  and 
the  fraction  of  the  total  problem  energy  converted  to  kinetic  energy  decreases. 
Finally,  for  horizontal  or  freerfleld  calculations,  no  spallation  occurs. 

Kinetic  energy  achieves  a  maximum  at  early  time,  then  declines  continuously 
as  material  motion  is  damped  In  the  Infinite  medium. 

The  kinetic  energy  content  of  spalled  mound  material  In  ballistic 
freefall  Is  approximately  constant  (Ignoring  gravitation).  Thus,  If  the 
entire  mound  entered  ballistic  freefall,  one  would  expect  the  kinetic  energy 
to  rise  to  a  peak  value  (during  the  spall  phase)  and  remain  constant  thereafter. 
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In  reality,  the  gas  acceleration  phase,  beginning  when  the  spall  wave  reaches 
the  cavity,  and  continuing  energy  transfer  from  the  cavity  at  late  time  cause 
the  medium  kinetic  energy  to  Increase  gradually.  The  gas  acceleration  effects 
spread  slowly  throughout  the  mound,  and  may  not  reach  the  near-surface  spall 
material.  In  addition,  material  Interactions  may  change  the  relative  energy 
contents  of  various  regions  of  the  mound.  For  example,  mound  recompaction  and 
reassembly  after  the  gas  acceleration  pulse  arrival  will  decelerate  lower 
parts  of  the  mound  and  transfer  energy  to  the  upper  sections  (see  the  previous 
discussions  of  recompaction  and  cavity  growth  damping).  Interactions  of  this 
type  are  of  particular  Interest  In  multilayer  cratering  problems. 

Other  energy-related  parameters  are  also  useful  In  describing  dynamic 
effects.  The  ratio  of  kinetic  energy  to  total  problem  energy  provides  a 
measure  of  the  efficiency  of  material  ejection.  Its  behavior  Is,  of  course, 
similar  to  kinetic  energy.  This  parameter  provides  a  convenient  means  of 
comparing  different  classes  of  problems.  The  ratio  of  kinetic  energy  to 
Internal  energy  (KE/IE)  likewise  gives  an  estimate  of  the  fraction  of  the 
locally  available  energy  (In  the  entire  mound  or  a  given  region)  which  goes 
Into  material  motion.  The  KE/IE  ratio  Is  high  for  shallowly  burled  events, 
and  extremely  high  In  the  spall  region  near  a  free  surface.  Material  pro¬ 
perties  also  influence  the  ratio:  other  factors  being  equal,  the  Internal  energy 
content  Is  highest  (and  KE/IE  ratio  Is  lowest)  for  stiff,  high-strength, 
shear-supporting  materials;  the  Internal  energy  content  decreases  for  weaker 
material  and  for  fluids.  The  KE/IE  ratio  becomes  small  for  deeply  burled 
events  In  all  media.  In  the  horizontal  problems,  material  motion  Is  strongly 
damped  at  very  late  times  and  the  Internal  energy  content  approaches  total 
problem  energy. 

Multilayer  dynamics  may  be  examined  by  dividing  the  problem  Into 
various  material  regions  and  examining  kinetic  energy  behavior  In  each  region. 
This  procedure  Is  particularly  appropriate  here  for  three  reasons:  (1)  Cratering 
material  ejection  for  underwater  events  depends  almost  entirely  on  velocities 
and  energies  In  the  underlying  rock.  The  dynamics  of  the  overlying  water  layer 
(after  Initial  shock  reflection,  damping,  and  recompaction  at  the  Interface 
are  complete)  do  not  Influence  the  dynamics  of  the  rock  material  and  the  rock 
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ejecta  trajectories.  Thus,  cratering  dynamics  after  the  Initial  strong 
interaction  phase  may  be  analyzed  solely  In  terms  of  the  underlying  rock 
layer.  Indeed,  the  entire  water  layer  for  most  of  the  shallow  events  In 
this  study  (Dwater  *  5  or  10  ft)  enters  freefall  when  the  water  surface 
reflection  spall  phase  Is  complete,  and  does  not  interact  again  with  the 
rock  layer  within  about  20  msec  (problem  termination  time).  (2)  The  problem 
sets  discussed  here  use  a  constant  rock  layer  depth,  In  combination  with 
variable  water  layers.  Hence,  the  rock  kinetic  energies  for  problems  In  a 
given  series  may  be  directly  compared  to  reveal  water  overburden  effects. 

Since  the  kinetic  energy  Is  velocity-sensitive,  damping  and  deceleration 
effects  which  influence  rock  layer  motion  will  be  readily  evident.  (3)  Gravita¬ 
tional  energies  are  negligible  (*  10‘3  x  KE)  over  the  time  and  distance 
scales  of  these  problems.  Even  for  the  rapidly-moving  upper  spalled  regions, 
gravitational  energy  accounts  for  only  one  percent  or  so  of  the  zone  energy 
at  late  times. 

Kinetic  energy  content  In  rock  as  a  function  of  time  for  all  problems 
Is  plotted  In  Figures  66-69  (solid  lines).  The  kinetic  energy  content  In  the 
overlying  water  layers  Is  also  given  for  comparison  (dashed  lines).  The 
various  problem  sets  are  displayed  separately:  coral  (Dr  =  30  ft)  in  Figure 
66,  coral  (Dr  =  36  ft)  In  Figure  67,  granite  (Df  =  36  ft)  In  Figure  68,  basalt 
(horizontal  and  36/00)  In  Figure  69a,  and  basalt  (36/10  and  36/36)  In  Figure  69b. 
The  basalt  diagrams  also  show  the  internal  energy  in  rock  and  water  and  the 
total  cavity  energy  for  all  problems.  Note  that  each  figure  gives  the  total 
problem  Input  energy  after  detonation  Is  complete  (horizontal  line  at 
6.65  x  10  Mbar-cm3)  and  the  arrows  at  the  bottom  of  the  figure  Indicate  the 
approximate  (acoustic)  shock  arrival  times  at  each  problem  Interface  or  boundary.* 
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All  of  the  kinetic  energy  figures  show  trends  consistent  with 
previous  discussions,  although  the  water  layer  damping  effects  are  some¬ 
what  exaggerated  due  to  the  extreme  velocity  sensitivity.  Primary  features 
of  the  kinetic  energy  behavior  are  summarized  as  follows: 

(1)  Kinetic  energy  curves  for  all  problems  In  a  given  medium  are 
Identical  until  reflection  occurs  at  the  rock  surface  or  rock-water  Inter¬ 
face.  The  kinetic  energy  for  rock  free-surface  problems  (30/00  and  36/00) 
then  rises  steeply  throughout  the  spall  phase  and  tends  to  level  off  at 
later  times.  The  rock  kinetic  energy  for  underwater  problems  continues  to 
decline  after  rock-water  Interface  reflection  due  to  water  layer  damping  at 
the  Interface.  However,  the  kinetic  energies  then  begin  to  Increase  as 
spallation  spreads  through  the  rock  mound  and  as  the  water  surface  reflection 
returns  to  relieve  damping  near  the  rock-water  Interface.  The  kinetic  energies 
finally  become  almost  constant,  with  slight  variations  at  late  times.  Horizontal 
problem  kinetic  energies  decline  steeply  at  late  times,  as  expected. 

(2)  The  kinetic  energy  In  water  rises  steeply  after  the  shock  wave  reaches 
the  rock^water  Interface.  This  Increase  continues  as  spall  reflection  occurs 
at  the  rock-water  Interface.  Ballistic  spallation  then  spreads  throughout 
the  water  layer  and  the  kinetic  energy  curves  level  off  (water  layer  In 
freefall).  Water  dynamics  are  most  clearly  displayed  In  the  ,Maep"  30/36  or 
36/36  problems.  Kinetic  energy  In  water  rises  steeply  as  the  shock  enters 
the  water  (time  30/00/1,  or  36/QOA),  Increases  more  slowly  as  the  shock  moves 
through  the  water  layer,  jumps  suddenly  as  the  shocx  is  reflected  at  the  free 
surface  (time  30/36A  or  36/36/1),  increases  rapidly  as  spallation  spreads 
through  the  water  layer,  then  levels  off  during  ballistic  freefall.  Further 
small  increases  may  occur  at  very  late  times  If  the  rock  layer  again  catches 
up  with  and  recompacts  the  bottom  part  of  the  water  spall  (gas  acceleration 
pulse  reaches  water  layer).  (3)  Kinetic  energy  content  in  rock  is  decreased 
significantly  by  the  presence  of  an  overlying  water  layer,  with  most  of  the 
decrease  occurring  even  for  the  shallowest  water  depths  (1^  =  5  or  10  ft). 
Additional  water  layer  depth  has  little  Influence  on  the  rock  kinetic  energy, 
and  likewise  does  not  Increase  the  water  layer  kinetic  energies  dramatically. 
(4)  The  rock  kinetic  energy  curves  are  slowly  varying  functions  beyond 
10-15  msec  (after  Initial  spall  phase  is  complete)  and  are  not  extremely 
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sensltlve  to  late-time  Interactions.  This  trend  clearly  shows  that  overall 
mound  dynamics  are  not  crucially  dependent  on  late-time  effects,  and  the  use 
of  one-dimensional  vertical  calculations  and  kinetic  energies  to  analyze 
mound  dynamics  Is  quite  appropriate.  (5)  The  late- time  effects  which  do 
appear  consist  primarily  of  a  slight  continuing  Increase  In  rock  energy  due 
to  gas  acceleration.  This  phase  Is  followed  by  a  very  late-time  decrease 
In  the  rock  kinetic  energy  for  coral  30/36  (Figure  66),  coral  36/36  (Figure  67), 
basalt  36/36  (Figure  69b),  and  granite  36/36  (Figure  68).  The  water  kinetic 
energy  then  Increases  for  these  same  problems  at  still  later  times.  This 
Interaction  Is  explained  as  recompaction  of  the  mound  following  the  gas 
acceleration  phase  (see  previous  discussion).  The  gas  acceleration  and  subse¬ 
quent  recompaction  effects  spread  slowly  upward  through  the  rock  mound, 
eventually  transferring  kinetic  energy  to  the  upper  regions  of  the  problem 
(water  layer).  Note  that  the  granite  kinetic  energy  curve  actually  shows 
two  gradual  Increases  followed  by  two  deceleration  phases  (Figure  68); 
these  correspond  with  the  two  separate  gas  acceleration  pulses  and  subsequent 
recompaction  phases  noted  earlier  (cf . ,  Figure  50b  and  others).  Only  the 
first  gas  acceleration  pulse  reaches  the  water  surface  layer  for  granite 
36/36  (single  late-time  increase  In  water  kinetic  energy). 

The  shallow  water  layer  problems  were  not  carried  late  enough  to 
show  these  mound  recompaction  effects  or  late-time  Interactions  with  the 
water  layer.  Recompacti on  would  occur  at  very  late  times  for  these  conflgur 
rations,  Lut  the  dispersion  of  the  rapidly-moving  spalled  rock  and  water 
layers  obviously  decreases  the  Intensity  and  speed  of  propagation  of  late-time 
effects. 

The  Internal  and  cavity  energies  are  also  shown  for  the  basalt  events 
(Figures  69a-b).  The  cavity  energy  content  at  late  times  Is  highest  for 
the  horizontal  calculation,  and  Is  lower  for  the  vertical  configurations 
(energy  Is  transferred  from  the  cavity  to  basalt  during  the  gas  acceleration 
phase).  The  cavity  energy  curves  are  similar  for  vertical  problems  36/00, 

36/10,  and  36/36  (as  were  the  cavity  radii  and  pressures  discussed  earlier). 
Again,  this  behavior  supports  the  conclusion  that  the  water  layers  have 
little  influence  on  gas  acceleration  and  late-time  close-range  effects  for 
basalt.  The  Internal  energy  curve  for  basalt  (horizontal  problem)  shows  a 


steady  Increase  at  late  times,  corresponding  with  the  decrease  in  kinetic 
energy.  Internal  energy  approaches  total  problem  energy  after  20  msec. 

The  basalt  vertical  problems  reveal  almost-constant  internal  energies  at 
late  times;  event  36/00  has  the  lowest  internal  energies,  while  the  under¬ 
water  problems  36/10  and  36/36  are  successively  a  little  higher.  However, 
the  difference  In  Internal  energy  between  these  three  problems  is  extremely 
small,  much  smaller  than  the  loss  of  kinetic  energy  for  the  underwater 
problems.  The  relationships  between  kinetic  and  Internal  energy  are  more 
clearly  shown  In  Figure  70,  which  plots  the  ratio  of  kinetic  energy  to 
internal  energy  (KE/IE)  for  both  basalt  and  water  layers.  The  ratio  in 
basalt  Is  highest  for  36/00,  increasing  dramatically  during  the  spall  phase. 

The  ratios  for  36/10  and  36/36  decline  due  to  water  damping  In  the  same  time 
period,  converging  on  lower  levels  at  late  time.  It  Is  evident  that  the 
presence  of  a  water  layer  causes  only  a  slight  increase  of  internal  energy 
but  a  large  decrease  of  kinetic  energy  and  KE/IE  ratio  in  the  underlying 
basalt.  Almost  all  of  the  kinetic  energy  lost  from  basalt  goes  into  the 
water  layer  (rather  than  into  internal  energy  of  basalt).  This  situation  is 
expected  to  occur  for  all  shallowly-burled  spall-dominated  problems.  The 
internal  energy  curves  for  the  water  layers  (Figure  69b)  reveal  interesting 
trends  not  visible  in  the  kinetic  energy  curves  alone.  Internal  energy 
increases  until  spall  reflection  at  the  water  surface  occurs,  then  decreases 
rapidly  as  spallation  spreads  throughout  the  water  layer,  finally  reaching 
a  constant  value  when  the  entire  water  layer  enters  freefall.  These  relation¬ 
ships  are  also  illustrated  in  Figure  70  (KE/IE  ratio  in  water  layer).  The 
ratio  of  kinetic  to  internal  energy  actually  remains  constant  or  declines 
slightly  as  the  shock  propagates  through  the  water  layer  (times  from  36/00A  to  36/1 0/v 
or  36/00a  to  36/36a).  However,  after  reflection  at  the  water  surface  (see 
times  after  36/10'  or  36/36:  ),  the  internal  energy  decreases  and  the  ratio  of  KE/IE 
increases  rapidly,  again  approaching  a  constant  value  at  late  time.  The 
late-time  internal  energy  content  in  the  rapidly-moving  water  spall  layer 
is  quite  small  and  the  ratio  of  KE/IE  is  very  large  as  expected.  The  water 
internal  energy  curve  for  basalt  36/36  shows  a  slight  increase  near  the  end 
of  the  problem,  corresponding  with  the  water  kinetic  energy  increase  at  the 
same  time  (gas  acceleration  pulse  and  recompaction  effects  reach  bottom  of  the 
water  layer). 


Ratio  of  kinetic  energy  to  internal  energy 


Some  of  the  representative  late- time  energy  data  have  been  abstracted 
from  Figures  66-69,  and  are  summarized  In  Tables  8  and  9.  The  tables 
also  list  corresponding  data  for  three  subsidiary  sample  problems  (lower 
three  entries  of  each  table).  These  results  will  be  discussed  In  Chapter  X. 

All  energies  are  selected  at  times  later  than  two  shock  transit  times  through¬ 
out  the  problem,  after  spallation  Is  complete  and  the  energy  curves  have"1eveled 
off."  In  most  cases,  a  two-value  range  Is  listed  for  each  quantity,  repre¬ 
senting  the  approximate  variability  at  late  times.  Table  8  lists  the  ratios 
of  kinetic  energy  (KE)  to  total  Input  energy  for  each  problem;  this  ratio 
Is  given  for  the  rock  layer  (KE  (rock)/total  energy),  water  layer  (KE  (water)/ 
total  energy),  and  both  layers  combined  (total  KE/total  energy).  Also, 
the  relative  partition  of  kinetic  energy  between  rock  and  water  Is  listed. 

Table  9  presents  the  ratios  of  kinetic  energy  to  Internal  energy;  the  ratio 
Is  compiled  for  the  rock  layer  [KE  (rock)/IE  (rock)]  and  for  the  water  layer 
[KE  (water)/IE  (water)].  Rock  values  list  an  approximate  range  of  varia¬ 
bility  at  late  times.  The  water  layer  ratios  tend  to  converge  on  a  constant 
value  at  late  times,  since  most  or  all  of  the  water  layer  Is  In  ballistic  freefall. 

The  ratios  of  kinetic  energy  In  rock  to  total  problem  energy 
(Table  8)  are  highest  for  coral  30/00  (about  0.45,  or  453!  of  total  energy 
Is  converted  into  kinetic  energy)  and  for  coral  36/00  (about  0.4  or  40%), 
becoming  progressively  lower  for  basalt  36/00  (38%)  and  granite  36/00  (15%). 

These  ratios  decline  somewhat  as  water  layers  of  Increasing  thickness  are 
added  to  the  problems,  and  the  ratio  (kinetic  energy  In  water)/(total  energy) 
Increases  correspondingly.  It  Is  quite  interesting  to  note  that  the  ratio 
(total  kinetic  energy)/(total  energy)  for  each  set  of  problems  declines  by 
only  about  0.02  to  0.06  (2  to  6%)  as  thicker  water  layers  are  added  (up  to 
Dw  =  36  ft).  Thus,  the  addition  of  water  does  not  greatly  decrease  the 
total  available  kinetic  energy,  in  spite  of  the  increase  In  overburden. 

The  addition  of  6  ft  of  coral  (30/00  to  36/00)  modified  the  total 
available  kinetic  energy  about  as  much  as  36  ft  of  water  overburden  for 
similar  problems.  Since  total  kinetic  energy  does  not  change  dramatically 
with  water  depth,  we  may  conveniently  deduce  water  layer  effects  by  comparing 
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Table  8  Ratio  of  Kinetic  Energy  to  Total  Problem  Input  Energy 
(late  times,  after  rock  layer  spallation  Is  complete)?* 


Basalt  35/00 


Basalt  36/10 


Basalt  36/36 


Granite  36/00 
Granite  36/10 
Granite  36/36 

Coral  41.31* ** *** 

j- - 

Coral  55.11* 

i - 

Basalt  39.83* 


*  Equivalent  overburden  comparison  calculations  (no 

**  Two  values  are  listed  In  each  case,  corresponding 
variability  at  late  times. 

***  Total  KE  =  Rock  KE  +  Water  KE 


Lower  Laye~r) 

.185 

.199 


rr  Layei 

.167 

.167 


.149 

1.0 

0 

.134 

.138 

.866 

.870 

.134 

.130 

.133 

.11 

.842 

.764 

.158 

.236 

.352 

.366 

.526 

.544 

.474 

.456 

.280 

.295 

.296 

.342 

.704 

.658 

.349 

.363 

.705 

.719 

.295 

.281 

water  layer)  -  see  text.  Chap  X. 
to  the  approximate  range  of 


TABlE  9.  Kinetic  Energy/ Internal  Energy  Ratios  at  Late  Times 
(  t  =  2  transit  times) . 


CONFIGURATION 


Coral  30/00 
30/05 
30/10 
30/36 
36/00 
36/05 
36/10 
36/36 


Basalt  36/00 
36/10 
36/36 


Granite  36/00 
36/10 
36/36 


Coral  41.31  c 
Coral  55.11* 


Basalt  39.83 


KE/IE  RATIO**  ROCK 

0.86 

-  0.95 

.49 

-  .58 

.46 

-  .48 

.46 

-  .3 

0.73 

-  0.78 

.46  - 

43  -  .48 

.45  - 

38  -  .41 

.45 

-  .33 

0.66 

-  0.74 

.44  -  . 

43  -  .48 

.44 

-  .42 

0.17 

-  0.20 

.13 

-  .16 

.13  - 

•  15  -0.1 

RATIO, 

LOWER  LAYER 

.331  -  .37$ 

.151  -  .168 

.396 

-  .480 

LATE-TIME  KE/IE 
RATIO  IN  WATER 


RATIO,  UPPER  LAYER 


* 

Equivalent  overburden  comparison  calculations  (no  water  layer) 
see  text.  Chapter  X. 

** 

A  range  of  values  is  listed,  corresponding  to  the  approximate 
variation  at  late  times,  after  two  transit  times. 
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the  kinetic  energy  partition  for  each  set  of  problems.  Using  average 
values,  the  fraction  of  kinetic  energy  partitioned  Into  the  rock  declines 
from  1.0  to  0.44  for  the  series  coral  30/00  to  coral  30/36  (the  kinetic 
energy  partition  In  water  correspondingly  Increases  from  0  to  0.56). 

The  partition  In  rock  decreases  from  1.0  to  0.45  for  coral  36/00  to  36/36, 
from  1.0  to  0.62  for  basalt  36/00  to  36/36,  and  from  1.0  to  0.80  for  granite 
36/00  to  36/36.  In  all  cases,  the  greatest  decrease  of  rock  kinetic  energy 
occurs  for  the  shallow  water  layer  (Dw  =  5  or  10  ft),  with  small  additional 
decreases  at  greater  water  depths.  The  water  layer  causes  the  largest 
kinetic  energy  decrease  In  rock  for  the  coral  events  (slight  Impedance 
mismatch  at  the  rock-water  Interface);  the  decrease  Is  much  less  severe 
for  basalt  and  for  granite  (gross  Impedance  mismatch,  efficient  tensile 
reflection  at  the  Interface).  This  situation  has  a  very  significant  effect 
on  the  KE  (rock)/total  energy  ratio:  The  "rock"  kinetic  energy  for  coral 
36/00  Is  greater  than  kinetic  energy  for  basalt  36/00,  but  the  basalt  kinetic 
energy  In  rock  becomes  greater  for  events  36/10  and  36/36.  The  Implications 
for  underwater  cratering  In  high  and  low  Impedance  rock  are  obvious.  Kinetic 
energy  partitioning  Into  the  rock  Is  even  more  efficient  for  the  granite 
36/10  and  36/36  events  (z  0.86  to  0.80).  Note,  however,  that  the  kinetic 
energy  (rock)/total  energy  ratios  for  all  granite  problems  are  extremely 
small  (=  0.15  to  0.1).  The  rock  kinetic  energy  for  granite  36/00  Is  less 
than  the  rock  kinetic  energies  for  coral  36/36  and  basalt  36/36,  with  deep 
water  layers.  This  result  once  more  verifies  that  the  high-strength  granite 
configurations  are  at  deeper-than-optlmum  burial  and  do  not  model  realistic 
cratering  events.  Nonetheless,  the  fundamental  result  obtained  here,  high 
kinetic  energy  partition  into  the  rock  layer  for  high  Impedance  rock,  is 
expected  to  hold  at  all  burial  depths. 

The  KE/IE  ratios  (Table  9)  show  trends  similar  to  the  above.  The 
ratios  are  highest  for  coral  30/00  and  36/00,  declining  for  the  other  media. 
The  KE/IE  ratios  for  high-strength  granite  are  quite  low,  again  revealing 
inefficient  material  ejection  and  excessively  deep  burial.  Most  of  the 
explosive  energy  Is  effectively  lost  as  Internal  energy  of  the  medium, 
and  does  not  contribute  to  material  ejection. 
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General  conclusions  of  the  kinetic  energy  analysis  may  be  stated 
as  follows:  (1)  Presence  of  a  water  layer  Is  more  Important  than  water 
depth  In  determining  overall  mound  dynamics  (even  for  relatively  deep 
water  layers,  Dw/Dr  =  1.0);  (2)  The  kinetic  energy  partitioned  into  rock 
depends  primarily  on  the  acoustic  mismatch  at  the  rock-water  interface, 
and  only  secondarily  on  water  depth;  (3)  Increasing  water  layer  depth 
does  not  rapidly  modify  the  total  kinetic  energy  available  to  eject 
material;  by  contrast,  the  available  kinetic  energy  decreases  precipitously 
with  Increasing  burial  depth  In  a  single-layer  rock  medium;  and  (4)  Following 
conclusions  (1)  -  (3),  the  water  layer  should  never  be  treated  as  an  "equivalent 
weight"  of  rock  overburden  In  cratering  problems. 

VIII.  SPALL  VELOCITIES-ACOUSTIC  ANALYSIS 

This  chapter  compares  SOC-calculated  shock  propagation  and  spall 
velocities  In  the  water  layer  with  the  results  of  a  simple  acoustic  trans¬ 
mission  analysis.  It  Is  assumed  that  the  shock  front  Is  transmitted  and 
reflected  as  a  low-amplitude  (acoustic)  wave  at  the  rock-water  Interface. 

Angle  of  refraction  and  radius  of  curvature  are  described  by  geometric  laws. 
Further,  the  front  Is  assumed  to  propagate  at  sonic  velocity  (linear  Elastic 
medium)  and  with  a  constant  peak  attenuation  rate  In  the  water  layer.  These 
assumptions  are  reasonable  because  the  peak  amplitudes  near  the  Interface 
are  moderate  (P  =  1  kbar)  and  peak  waveforms  arc-  relatively  well-resolved. 

For  acoustic  waves,  the  angle  of  reflection  at  an  interface  Is  equal 
to  the  angle  of  Influence.  The  angle  of  the  refraction  Into  the  overlying 
layer  Is  given  by,  c 

cos  ew  =  ^  COS  6r 

where  "e"  Is  measured  between  the  direction  of  shock  propagation  and  the 
Interface;  e  =  9C°  for  normal  Incidence. 

0w  =  angle  of  refraction  Into  upper  (water)  layer 
8r  *  angle  of  incidence  In  lower  (rock)  layer 
Cyf  =  front  velocity  of  propagation  In  upper  layer 
Cr  *  front  velocity  of  propagation  in  lower  layer. 
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Continuity  of  pressure  and  particle  velocity  across  the  interface  may 

applied  to  derive  transmission  and  reflection  factors.  The  reflected 

pressure  at  the  interface  Is  determined  to  be, 

pc  sine  -pc  sine 

P  (reflected)  -  P, 


be 

peak 


Pf  =  Incident  peak  pressure  In  lower  (rock)  layer 
pCr  *  acoustic  impedance  of  lower  layer 
pCty  *  acoustic  Impedance  of  upper  layer 

Similarly,  the  peak  particle  velocity  after  reflection  at  the  Interface  is, 

V  (reflected)  =  V  pCwsin9r~pCrsin9w 

pcwslner+pCrs1n6w 

Vr  =  Incident  peak  particle  velocity  in  the  lower  (rock)  layer. 

The  quantity  "sin  0"  Is  of  course  equal  to  1.0  for  normal  incidence.  Note 

that  the  velocity  or  pressure  reflection  factor  is  negative  at  a  hlgh-to- 

low  Impedance  Interface  (pCr>pCw,  Cf>Cw,  0w>0f).  Thus,  an  Incident  coirpresslon 

wave  Is  reflected  as  a  downward-travelling  tensile  wave  (upward  velocity  component). 

The  spall  velocity  at  the  Interface  is  obtained  by  summing  the  vertical  velocity 

components: 

Vspall  (rock)  =  Cvr+v  (reflected)]  s1n0p  <  2Vf 

Note  that  this  rock  interface  spall  velocity  is  a  peak  value;  velocity  near 
the  interface  may  be  damped  by  the  water  layer  at  later  times,  as  the  trans¬ 
mitted  shock  propagates  away  from  the  Interface  and  velocities  In  the  water 
layer  decrease. 

The  peak  pressure  transmitted  Into  the  upper  layer  at  the  Interface  is, 
2pCws1n6r 

Pw  “  Pr  ^ptwsiner+pCrs1n0w^ 

Pw  =  transmitted  peak  pressure  in  the  upper  (water)  layer 

The  term  In  brackets  will  be  denoted  as  the  pressure  transmission 
factor  Into  the  water  layer.  The  transmission  factor  decreases  for  increasing 
Impedance  mismatch  between  rock  and  water  and  for  Increasing  Incident  angle  "0r"; 
Figure  71  presents  a  plot  of  the  rock-water  pressure  transmission  factor  as  a 
function  of  angle  for  the  coral,  basalt,  and  granite  media  used  in  this  study. 


ic  pressure  transmission  factor  as  a  function  of  angle  of  incidence 
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The  Impedances  of  rock  and  water  differ,  so  the  velocity  transmission 
factor  Into  the  water  layer  Is  not  the  same  as  the  pressure  transmission 
factor.  Peak  particle  velocity  and  peak  pressure  at  the  shock  front  In  the 
water  layer  are  related  by  the  acoustic  equation, 

P  =  pC  V 
w  M  w  w 

Vw  *  peak  particle  velocity  at  the  front. 

Thus,  the  transmission  equation  may  be  expressed  In  terms  of  velocity  as, 


'w 


pCr  2pCws1n0r 

Vr  pCwsiner+pCrslnew 


Note  also  that  either  peak  velocity  or  peak  pressure  at  the  shock  front  may 
be  determined  If  the  other  quantity  Is  known.  For  the  properties  of  water 
assumed  In  these  problems, 

pCw  5  66.3  Eil^p  =  217.5 
and  the  peak  pressure  Is  given  by, 


Pw  *  pCwVw  =  217-5VW  (Vw  In  m/sec,  Pw  In  psl) 

This  relationship  applies  at  the  shock  front  (peak  values),  but  not  to 
spall  velocities  In  the  surface  cutoff  layer  adjacent  to  an  Interface,  where 
the  tensile  reflection  truncates  the  peak  pressure  of  the  Incident  wave. 

For  water.  It  has  been  shown^®-^  that  peak  attenuation  of  weak  pulses 
and  even  of  relatively  strong  shock  fronts  occurs  at  the  rate  of  where  R 

Is  the  range  from  center  for  a  spherically  divergent  front,  or.more  generally, 
the  radius  of  curvature  of  the  front.  Hence,  attenuation  of  the  transmitted 
shock  In  the  water  layer  may  be  adequately  modeled  by  a  power  law  dependence. 

Spall  launch  velocities  of  the  water  surface  may  likewise  be  calculated, 
since  perfect  tensile  reflection  occurs  at  a  free  surface  and  the  vertical 
component  of  particle  velocity  Is  doubled: 

Vspal  1  (water)  =  2Vat  surface )s1new  =2.P(wate*%  at  surface)  s1n0^ 

w 

The  peak  pressures  and  velocities  at  or  above  the  rock -water  Interface 
and  the  peak  surface  spall  velocities  may  thus  be  acoustically  simulated  If 
the  following  parameters  are  available:  Peak  pressure  or  particle  velocity  of 


\ 
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the  shock  front  Incident  at  the  Interface,  angle  of  incidence,  impedances 

and  sonic  velocities  of  the  two  layers,  and  attenuation  rate  in  the  upper 

layer.  All  of  these  quantities  are  known  except  the  peak  pressure  (particle 

velocity)  at  the  Interface,  which  must  be  derived  from  SOC  calculations  for  a 

T55l 

given  medium.  Preston  and  Wheel erLJ  have  demonstrated  that  peak  shock 
front  parameters  in  rock  may  be  accurately  fitted  over  a  broad  range  interval 
by  power  law  equations.  However,  separate  equations  are  needed  for  the  "high 
pressure"  regime  (range  <  a  fixed  transition  range)  and  the  "low  pressure" 
or  acoustic  regime  (range  >  transition  range).  Approximate  relations  have 
been  derived  for  the  coral,  basalt,  and  granite  calculations.  It  was  convenient 
In  this  case  to  fit  peak  particle  velocities  at  the  shock  front.  The  horizontal 
calculations  were  first  used,  but  results  proved  slightly  Inaccurate  (z  10%  low) 
due  to  the  coarse  zoning.  Small  corrections  were  applied  to  better  match  the 
finely-zoned  vertical  calculations.  Fitted  equations  take  the  form, 
bl 

Vj  =  ajR  (for  Range  R  <  transition  range  Rtr) 
v2  -  a2Rb2  (for  R  >  Rtr) 

V,  =  V2  at  Rtr 
a  and  b  are  constants 

Peak  velocity  relations  for  the  three  media  ar?  presented  in  Table  10.  Note 
that  attenuation  rates  "b^1  and  "b2"  In  each  midlum  are  always  more  rapid  for 
the  close-in  region  ("bi")  than  at  long  ranges  ("b2").  Comparing  materials, 
the  attenuation  rates  in  both  regions  are  most  rapid  for  granite.  Power  law 
fits  proved  slightly  Inaccurate  at  very  close  ranges  (R  <  10  ft)  for  these  SOC 
calculations;  all  of  the  events  discussed  here  were  burled  more  deeply  than 
10  ft  so  the  close-range  inaccuracies  have  no  effect  at  the  Interface. 

Given  a  suitable  model  of  the  rock-layer  peak  particle  velocities, 
the  velocity  or  stress  at  any  Interface  location  may  be  determined  and  acoustic 
calculations  for  the  ray  transmitted  at  that  location  may  be  performed.  Note 
that  these  calculations  are  applicable  for  any  Intercept  angle  (as  well  as 
for  the  vertical  direction),  thus  giving  a  two-dimensional  picture  of  shock 
transmission  into  the  upper  layer.  Figure  72  schematically  depicts  the 


rock 
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transmission  of  a  vertical  ray  and  an  oblique  ray  in  a  two-layer  configu¬ 
ration.  (The  figure  shows  a  third  overlying  layer  of  air  as  well;  we  will 
treat  the  water  surface  as  a  free  surface  for  present  purposes.) 

Acoustic  calculations  have  been  performed  for  events  at  a  variety 
of  rock  and  water  layer  depths,  in  coral,  basalt,  and  granite.  Velocity 
fits  from  Table  10  provided  the  input  rock  velocities  "Vr".  Results  of 
these  calculations  may  be  conveniently  compared  with  the  "SOC"  interface 
and  surface  spall  velocities  (vertical  direction  only).  However,  It  should 
be  remembered  that  the  SOC  one-dimensional  calculations  necessarily  impose 
the  condition  of  a  spherically-divergent  geometry.  No  corrections  are 
Included  for  refraction  effects  and  change  In  wavefront  curvature  at  an 
Interface.  The  error  thus  Introduced  will  be  negligible  very  close  to  the 
interface  (thin  water  layer),  and  will  become  progressively  larger  as  the 
transmitted  shock  propagates  away  from  the  Interface.  In  order  to  compare 
"SOC"  and  "acoustic"  velocities.  It  was  necessary  to  make  two  sets  of  acoustic 
calculations.  The  first  set  assumes  spherical  divergence  of  the  transmitted 
shock  in  the  water  layer,*  and  will  be  denoted  as  the  "linear"  or  "uncorrected" 
set.  The  second  set  of  calculations  Includes  the  correct  treatment  of  refraction 
and  change  In  divergence  at  the  interface  for  all  raysJ27,26^  These  calculations 
will  be  known  as  the  "corrected"  set.  The  two  sets  are  Identical  at  the  rock 
interface,  but  the  water  surface  spall  velocities  deviate  progressively  for 
increasing  water  layer  depths. 

Spall  velocities  above  shot  point.  Table  11  compares  the  spall  velo¬ 
cities  at  the  rock-water  interface  and  at  the  water  surface  for  the  SOC  calcula¬ 
tions  and  acoustic  calculations.  All  of  these  velocities  are  in  the  vertical 
direction,  directly  above  shot  point.  In  addition  to  peak  values,  the  SOC 
tabulations  list  a  range  of  spall  velocities  for  the  near-surface  high-velocity 


*  For  example,  the  vertical  (normal)  ray  transmitted  into  the  water  layer 

attenuates  as  [D/(Dr+Dw)]l *13  beyond  the  rock-water  interface.  No  change  in 
the  radius  of  curvature  of  the  wavefront  occurs  at  the  interface.  Corres¬ 
ponding  equations  may  be  written  for  other  angles  of  incidence  in  the  case 
of  "two-dimensional"  calculations  (oblique  incidence). 
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spall  region  (see  Appendix  A).  The  first  column  of  SOC  results  gives  the 
rock  Interface  spall  velocities  for  each  configuration.  The  next  two 
columns  l*st  the  water  surface  spall  velocities  (In  the  high-velocity  spall 
region,  and  over  an  arbitrary  1.5  m  depth  Interval  below  the  surface). 

Two  sets  of  velocities  are  listed  In  the  case  of  granite,  for  the  coarse  and  fine 
water  zone  SOC  calculations.  The  acoustic  section  of  the  table  gives  the 
rock  interface  spall  velocity  (Identical  for  both  uncorrected  and  corrected 
calculations),  the  water  surface  spall  velocity  for  linear  or  uncorrected 
calculations,  and  the  water  surface  spall  velocity  for  corrected  calculations. 

The  rock  interface  velocities  agree  quite  closely  between  SOC  and  acoustic 
results,  as  would  be  expected  (SOC  calculations  were  used  to  model  peak 
shock  front  parameters  In  the  rock  layer).  However,  the  SOC  Interface  velocity 
for  granite  (36/10  and  36/36)  Is  somewhat  higher  than  the  acoustically-calculated 
value  for  reflection  at  a  rock-water  Interface.  Comparing  the  water  surface 
peak  spall  velocities  for  coral  and  basalt  events,  we  find  that  the  acoustic 
linear  (uncorrected)  values  are  In  excellent  agreement  with  the  corresponding 
SOC  results.  The  corrected  acoustic  water  surface  velocities  are  much  higher 
than  the  uncorrected  and  SOC  results  (refraction  at  a  high-  to  low-velocity 
interface  decreases  divergence  of  the  shock  front,  thus  Increasing  water 
layer  pressures  and  spall  velocities).  Largest  errors  occur  for  the  deep 
water  layer  configurations.  The  errors  for  0.  =  36  ft  are  approximately  21-26% 
(coral)  or  36%  (basalt,  strong  Impedance  mismatch).  The  acoustically- 
calculated  water  surface  velocities  for  granite  configurations  are  all 
significantly  higher  than  the  SOC  values;  agreement  Improves  somewhat  for 
the  "fine"  water  zone  calculations.  In  summary,  the  uncorrected  acoustic 
calculations  match  the  SOC  results  for  basalt  and  coral  configurations,  but 
not  for  granite.  The  corrected  acoustic  velocities  are  higher,  showing  that  the 
assumption  of  spherical  divergence  causes  errors  in  the  SOC  water  layer  calculations. 
Errors  are  greatest  for  deep  water  layers  (far  from  the  interface)  and  for 
a  large  acoustic  mismatch  at  the  Interface. 

The  disagreement  between  acoustic  and  SOC  calculations  for  granite 
is  closely  related  to  the  attenuation  effect  mentioned  previously  (Chapter  VI, 
peak  velocities  and  pressures).  It  was  noted  that  peak  velocities  beyond 
the  granite-water  Interface  suffered  a  sudden  drop,  then  continued  to  attenuate 
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rapldly  throughout  the  water  layer.  The  use  of  "fine11  (half-size)  water  zones 
reduced  the  sudden  drop  but  rapid  attenuation  remained.  It  Is  not  surprising  that 
the  constant  "R"1-13"  attenuation  rate  of  the  acoustic  calculations  generates 
higher  water  spall  velocities  than  the  SOC  calculations.  The  water  layer 
attenuation  rates  for  "SOC"  results  are  examined  In  Table  12.  Table  12a  pre¬ 
sents  the  approximate  attenuation  rate  of  peak  particle  velocity  at  the  shock 
front  (excluding  the  narrow  surface  spall  layer  where  peaks  are  reinforced 
by  the  tensile  reflection).  The  water  attenuation  above  coral  and  basalt  varies 
approximately  as  R  to  R  1-2,  rates  which  are  slightly  greater  than  the  "acoustic" 
assumption.  The  rate  above  granite  Is  R’1-55.  The  attenuation  rates  of  water 
surface  spall  velocity  between  the  various  calculations  (l.e. ,  coral  36/05- 
36/10-36/36.  etc.),  show  a  similar  trend  (Table  12b).*  The  granite-water 
calculations  once  more  show  an  attenuation  rate  on  the  order  of  R-^*5.  This 
question  will  be  further  examined  In  Chapter  X. 

More  detailed  results  of  the  acoustic  calculations  for  the  vertical 
direction  are  presented  In  Figures  73-79: 

Figure  73  shows  the  rock  Interface  peak  spall  velocities  as  a  function 
of  burial  depth  In  rock.  "Dr".  Acoustically-calculated  velocities  are  shown 
for  the  case  of  no  water  layer  (free  surface)  and  for  the  case  of  an  overlying 
water  layer  (partial  reflection  at  interface,  spall  velocities  are  lower,  decrease  does 
not  depend  on  water  layer  depth).  SOC-calc.ilated  velocities  (with  and  without 
water  layers)  are  Included  for  comparison.  Some  additional  SOC  velocities 
(from  Chapter  X)  have  also  been  added  to  demonstrate  agreement  at  other  burial  depths. 

Figures  74-77  depict  the  corrected  water  surface  spall  velocities 
as  a  function  of  water  layer  depth.  The  water  surface  velocities  for 
all  three  media  are  given.  Figure  74  Is  for  a  rock  burial  depth  Dr  =  24  ft. 

Figure  75  for  a  depth  Dr  =  30  ft,  Figure  76  for  Dp  =  36  ft,  and  Figure  77  for 
Dr  =  42  ft.  This  encompasses  the  range  of  burial  depths  comnonly  u: 2d  with  10-ton 
excavation  charges.  Note  that  the  water  spall  velocities  are  generally  highest  for 
coral,  lower  for  basalt,  and  lower  still  for  granite.  However,  coral  calcula¬ 
tion  velocities  decrease  most  rapidly  with  Increasing  water  layer  depth.  This 
★ 

Free-surface  spall  velocity  Is  approximately  twice  the  local  peak  particle 
velocity  at  the  shock  front;  therefore,  the  attenuation  rate  of  surface  spall 
velocity  between  the  various  water-layer  depth  calculations  should  match  the 
attenuation  rate  of  peak  particle  velocity  In  the  water  layer. 
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Figure  7b.  Water  surface  spall  velocities  as  a  function  of  water  depth  for 


Figure  75.  Water  surface  spall  velocities  as  a  function  of  water  depth,  for 


Water  surface  spall  velocities  as  a  function  of  water  dep 


trend  again  Illustrates  the  Importance  of  Interface  refraction  effects: 

Strong  refraction  at  the  basalt  or  granite  Interface  causes  the  shock 
front  to  diverge  more  slowly.  Spall  velocities  far  from  the  interface 
correspondingly  decrease  less  rapidly.  This  difference  between  coral  and  basalt 
behavior  would  not  be  apparent  for  uncorrected  calculations  (spherical 
divergence). 

Figures  78  and  79  show  both  the  corrected  and  linear  (uncorrected) 
water  surface  spall  velocities  for  coral,  basalt,  and  granite.  The  SOC- 
calculated  peak  spall  velocity  data  points  are  plotted  for  comparison. 

The  coral  (C)  and  basalt  (B)  points  agree  with  the  linear  acoustic  calcu¬ 
lations,  as  discussed  above.  The  granite  points  (6^  =  fine  zone,  G£  = 
coarse  zone)  are  low.  In  order  to  examine  attenuation  effects,  the  linear 
(uncorrected)  acoustic  calculations  for  granite  were  redone  with  a  water- 
layer  attenuation  rate  of  R-1.55.  The  resultant  curve  is  shown  In  Figure  79. 
Agreement  with  the  "fine-zone"  granite  calculations  Is  greatly  Improved, 
with  SOC  data  points  slightly  below  the  revised  acoustic  calculations. 

Empirical  data.  Water  surface  velocities  have  been  measured  for 
the  Tugboat  Phase  I  underwater  events  In  coral. Data  are  listed  In 
Table  13.  Four  of  these  events  were  of  one-ton  yield;  the  rock  burial  depths 
"Dr"  and  water  layer  depths  "Dw"  have  been  cube-root  scaled  to  10  tons.  In  order 
to  render  the  results  comparable  with  calculations.  All  velocity  measurements 
were  made  using  late-time  low-resolution  photographic  techniques  (time  >  90  msec, 
mound  height  -  10  to  100  ft).  For  this  reason,  the  data  cannot  be  directly 
compared  with  early-time  peak  surface  spall  velocities.  The  water  medium 
has  very  low  tensile  strength  and  spalls  upward  In  thin  layers  when  a  shock 
wave  Is  reflected  at  the  surface.  These  thin  layers,  unlike  massive  rock 
spall,  suffer  enormous  air  drag  deceleration  at  early  times.  Indeed,  the 
Initial  surface  deceleration  may  be  on  the  order  of  hundreds  of  g's,  even 
for  greater  ranges  and  lower  velocities  than  those  encountered  here.^-56^ 

The  rapid  deceleration  ceases  and  a  stable  freefall  configuration  Is  achieved 
only  after  an  appreciable  thickness  of  water  Is  accreted  to  the  bottom  of 
the  spall.  Thus,  velocities  measured  over  distance  scales  on  the  order  of 
many  feet  and  time  scales  much  greater  than  the  initial  spall  launch  time 


Measured  Ground  Zero  Water  Surface  Vertical  Velocities  At 
Mound  Height  >10  to  20  ft). 
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will  bo  representative  of  the  "average"  accreted  velocity  throughout 
several  feet  of  the  water  mound  (rather  than  the  Instantaneous  peak  very 
near  the  free  surface). 

The  SOC-calculated  rarvge  of  water  peak  spall  velocities  In  the  high- 
velocity  near-surface  region  and  In  the  top  1.5  m  of  the  water  layer  are 
plotted  against  water  layer  depth  In  Figures  80  (coral,  Dp  *  30  ft),  81 
(coral  Dr  -  36  ft),  82  (basalt),  and  83  (granite).  The  corrected  acoustic  peak 
spall  velocity  curves  are  also  shown.  Data  points  for  the  Tugboat  coral 
events  at  the  closest  comparable  rock  depths  "Dr"  are  presented  in  Figures  80 
(Dr  s  30  ft)  and  81  (Df  *  36  ft).  It  Is  evident  that  the  measurements  fall 
below  the  "SOC"  and  "acoustic"  peak  surface  velocities,  and  are  more  nearly 
comparable  with  the  velocities  about  1  to  1.5  m  below  the  surface.  Note  that 
event  "ID"  In  Figure  81  may  not  be  legitimately  compared  with  the  other  data 
and  calculations,  since  it  was  at  a  very  great  scaled  depth  of  Dp  *  44.4  ft. 

The  differences  between  Tugboat  data  and  the  corrected  acoustic  calcula¬ 
tions  are  more  accurately  compared  in  Figure  84.  Acoustic  calculations  were 
performed  for  the  precise  (scaled)  rock  and  water  depths  of  each  Tugboat  con¬ 
figuration.  Figure  84  plots  the  apparent  error  "V(acoust1c)/V(measured)." 

The  acoustic  peak  velocities  are  about  1.9  times  as  high  as  the  measurements 
for  the  one-ton  events  (IA-ID),  about  1.37  times  as  high  for  the  10-ton 
event  (IE).  This  fact  In  itself  indicates  a  photographic  time-scale  and 
resolution  problem,  since  accretion  occurs  over  the  smallest  distance  and 
shortest  time  Interval  for  the  "small  scale"  one-ton  events.  However, 

Figure  84  again  indicates  that  measured  velocities  are  representative  of  the 
top  layers  of  water  spall  rather  than  the  peak  surface  velocity. 

A  final  comparison  between  calculation  and  experiment  is  presented 
in  Figures  85-87.  These  Figures  show  the  SOC-calculated  spall  velocity  (or 
residual  velocity)  profiles  throughout  the  water  mound  after  spallation  is 
complete.  The  corrected  acoustic  surface  velocities  are  also  shown  as 
crosses  (somewhat  higher  than  the  SOC  surface  velocities).  The  SOC  calcu¬ 
lations  reveal  the  characteristic  features  of  slowly  changing  velocities  deep 
within  the  mound,  followed  by  a  sharp  increase  in  spall  velocity  within 
about  2  m  of  the  free  surface.  The  measured  velocity  for  Tugboat  IB  is 


calculated  water  surface  layer  spall  velocities  for  basalt 


Figure  84.  Water  surface  spall  velocity,  apparent  errors:  V  (acoustic)/V  (measured) 
vs  rock  depth,  D  (scaled  to  10  tons). 
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Indlcated  In  Figure  85  (compare  with  coral  30/10).  Measurements  for  Tugboat 
IC  and  IE  are  marked  In  Figure  86  (compare  with  coral  36/10).  The  late¬ 
time  measured  velocities  Indeed  appear  consistent  with  the  Interpretation 
of  "averaged"  velocities  throughout  the  accreted  water  spall  mound.  Water 
velocity  profiles  for  basalt  and  granite  calculations  are  shown  In  Figure  87 
(no  measurements  available).  These  data  are  slightly  less  accurate  than 
coral  due  to  the  larger  refraction  errors  (high -impedance  media)  and  the 
resolution  difficulties  for  the  granite-water  calculations. 

"Two-dimensional"  calculations.  As  noted  previously,  the  acoustic 
transmission  equations  may  be  used  to  simulate  shock  transmission  for 
Intercept  angles  other  than  the  normal  (vertical)  direction.  Thus,  trans¬ 
mitted  pressure  and  velocity  may  be  determined  anywhere  within  the  water  layer; 
peak  vertical  spall  velocities  may  likewise  be  calculated  for  any  location 
along  the  rock  Interface  and  water  surface,  giving  an  approximate  "two- 
dimensional"  picture  of  mound  growth.  Corrected  acoustic  calculations  have 
been  performed  for  a  variety  of  configurations;  sample  results  are  shown  In 
Figures  88  (coral,  Dp  -  30  ft),  89  (coral,  Dr  =  36  ft),  90  (basalt,  =  36  ft), 
and  91  (granite,  Dr  =  36  ft).  Each  diagram  plots  the  peak  vertical  spall 
velocity  profiles  of  the  rock  free  surface  (no  water),  rock  Interface  (over- 
lying  water  layer.  Imperfect  reflection  at  the  Interface),  water  surface 
(for  Dw  =  10  ft),  and  water  surface  (for  Dw  *  36  ft).  Note  that  the 
water  surface  spall  velocities  frequently  exceed  the  underlying  rock  spall 
velocities  at  the  same  surface  range,  even  for  the  36  ft  water  layers 
(Figures  89-91).  This  situation  Is  a  result  of  the  low  Impedance  of  water 
and  of  the  Interface  refraction  effects  which  decrease  divergence  In  the 
water  layer.  These  acoustic  calculations  give  only  peak  particle  velocities  at 
the  shock  front  and  spall  velocities  at  an  Interface;  thus,  they  are  not  equiva¬ 
lent  to  true  two-dimensional  dynamic  calculations  used  to  simulate  the  entire 
history  of  the  mound  and  predict  crater  dimensions.  Currently,  there  are  no 
available  mound  growth  data  or  two-dimensional  dynamic  calculations  with  which 
Figures  88-91  may  be  compared. 


Figure  88.  Peak  surface  spall  velocity  profiles  for  coral,  =  30  ft. 
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Figure  90-  Peak  surface  spall  velocity  profiles  for  basalt,  D  =36  ft. 
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IX.  AIRBLAST 

In  the  foregoing  discussion*  the  water  surface  above  a  submerged 
cratering  event  has  been  treated  as  a  free  surface.  This  approximation  Is 
a  reasonable  one  to  make  when  studying  the  Instantaneous  peak  spall  velocities 
or  overall  dynamics  of  the  water  layer  at  early  times,  since  water  has  a  much 

greater  density  and  Impedance  than  air.  However,  the  alrblast  shock  trans¬ 
mitted  Into  the  overlying  air  Is  also  of  Interest  for  safety  predictions. 
Window  breakage  and  undesirable  noise  effects  are  often  limiting  factors 
In  designing  underwater  excavation  shots.  This  section  discusses  alrblast 
prediction  for  underwater  events  and  provides  a  comparison  with  empirical 
data. 


Explosive  cratering  events  usually  generate  two  alrblast  pulses. 

The  first  of  these  Is  the  spall -Induced  shock,  coupled  Into  the  air  by 
the  piston  effect  of  the  spalled  rock  (or  water)  surface.  The  Intensity 
of  the  spall-induced  alrblast  Is  closely  related  to  the  peak  surface  spall 
velocity,  and  ray  be  accurately  predicted  If  the  surface  velocity  profile 
Is  known.1  •  J  The  second  or  gas- vent-induced  pulse  originates  when  the 
material  mound  begins  to  disintegrate  and  cavity  gases  vent  to  the  surround¬ 
ing  atmosphere.  The  Intensity  of  gas  vent  alrblast  Is  related  to  the  pres¬ 
sure  and  volume  of  cavity  gases  at  vent  time;  this  pulse  may  be  weak  or 
absent  If  cavity  pressure  Is  near  or  below  atmospheric  pressure  when  venting 
occurs  (1. e. ,  late  vent  time  or  deeply  burled  event).  It  has  been  observed 
that  vent  times  are  relatively  late  and  gas  vent  alrblast  Is  negligible  for 
events  near  or  below  optimum  burial  depth  In  saturated  and  saturated  sub¬ 
merged  media.  This  trend  Is  presumably  attributable  to  the  coherency  of  the 
saturated  material  mound,  which  efficiently  contains  the  cavity  gas  bubble 
until  late  time.  Since  gas-vent  alrblast  Is  negligible  for  typical  cratering 

events  In  saturated  media,  this  discussion  will  be  concerned  with  the  predominant 
spall-induced  pulse. 

Due  to  the  very  long  duration  or  wavelength  of  the  alrblast  signal 
from  a  large  excavation  event,  window  breakage  and  other  distant  damage  effects 
are  related  primarily  to  the  peak  amplitude  (rather  than  Impulse  or  duration). 
Thus,  safety  criteria  may  be  established  on  the  basis  of  the  peak  overpressure 
AP  (greatest  increase  above  the  local  ambient  air  pressure).  Surface  spall 
velocities  may  be  used  to  predict  the  peak  alrblast  overpressure,  both  In  the 
Immediate  vicinity  of  the  spall  mound  ("near  field'),  and  at  ranges  beyond  this 
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vicinity  ("far  field"). l-57^  Specifically,  the  peak  spall-induced  over¬ 
pressure  "AP0"  directly  above  the  center  of  the  mound  (surface  ground  zero) 
is  determined  from  the  acoustic  equation. 

4Po  5  pCa  vspa1 1  <at  SGZ)  El"- 
pCa  =  acoustic  impedance  of  air 

h  0.0184  ^s-i^ec  at  sea  level 

^spall^7)  s  >,ea*<  sPa^  velocity  of  the  mound  at  surface  ground  zero. 

The  far  field  peak  airblast  overpressure  "AP"  as  a  function  of  horizontal 
range  "R"  is  given  by, 

DOB 

AP  *  APQa  -jp  where  DOB  =  total  depth  of  burial  of  detonation 

or,  for  a  two-layer  configuration, 

D  +D 

AP  =  APQa  Eqn.  II 

AP0  ■  peak  local  airblast  overpressure  at  surface  ground  zero 
(near  field). 

a  s  dimensionless  "source  strength"  parameter  which  is  related  to 
the  overall  or  integrated  mound  velocity  profile. 

R  =  horizontal  range  from  surface  ground  zero  (measured  along  the 
rock  or  water  surface). 

Note  the  different  use  of  the  symbol  "R"  in  this  connection.  The  above 
equation  applies  only  in  the  horizontal  direction  (i.e.,  close  to  surface 
level),  and  only  at  far  field  ranges  beyond  the  rapidly-rising  spall  mound 
[R»(Dr+Dw)]. 

As  might  be  expected,  the  source  strength  parameter  "a"  is  dependent 
on  the  width  of  the  surface  velocity  field  and  the  integral  of  the  spall  velocity 
profile  throughout  the  field.  For  the  cases  under  consideration  here,  "a" 
may  be  satisfactorily  approximated  by  a  cross-sectional  integral  of  the 
(symmetric)  surface  spall  velocity  profile: 
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Vspal1^R^  =  Peak  surface  spall  velocity  as  a  function  of  range  "R" 
from  SGZ. 

This  relationship  has  been  found  to  be  accurate  for  explosive  cratering  events 
in  d  wide  variety  of  media.  In  practice,  the  profile  Integral  Is  carried 
c  ly  over  the  "high  velocity"  Inner  portion  of  the  velocity  field,  which  Is 
entirely  responsible  for  creating  the  airblast  pulse.  The  surface  velocity 
declines  rapidly  away  from  SGZ,  and  the  outer  "wings"  of  the  velocity  profile 
do  not  contribute  significantly  to  the  Integral. 

Since  the  impedance  of  air  Is  much  lower  than  that  of  water,  tensile 
reflection  of  an  Incident  compressional  shock  at  the  water  surface  (water- 
air  Interface)  Is  almost  perfect.  Therefore,  the  vertical  conponent  of 
particle  velocity  is  very  nearly  doubled  at  the  Interface,  and  water  surface 
velocities  calculated  on  the  basis  of  the  "free  surface"  assumption  may  be 
used  to  predict  airblast.  Water  surface  spall  velocities  for  the  underwater 
cratering  media  Included  In  this  study  are  obtainable  from  either  "SOC"  or 
acoustic  calculations.  The  acoustic  calculations  are  corrected  for  refraction 
effects,  and  can  easily  be  performed  for  a  variety  of  different  water  depths 
(only  ten  "SOC"  underwater  calculations  are  available).  Also,  acoustic 
calculations  may  be  used  to  generate  the  entire  surface  spall  velocity  profile 
(as  well  as  the  velocity  In  the  vertical  direction,  at  SGZ).  Accordingly, 
the  corrected  acoustic  results  have  been  adopted  for  airblast  predictions. 

Water  surface  spall  velocities  as  a  function  of  water  layer  depth 
were  presented  In  Figures  74  through  77.  These  velocities  have  been  used  to 
calculate  "AP0"  (Eqn.  I).  Likewise,  the  acoustic  water  surface  spall  velocity 
profiles  were  Integrated  to  obtain  "a"  (Eqn.  III).  (Typical  examples  of  the 
velocity  profiles  were  shown  in  Figures  88-91.)  The  results  of  these  calcu¬ 
lations  are  displayed  in  Figures  92-94: 

Figure  92  shows  "4P0"  (psl)  and  "a"  (dimensionless)  as  functions  of  overlying 
water  layer  depth  for  events  buried  in  coral  (Dr  ■  30,36,  and  42  ft).  The  "aP  > 
scale  is  calibrated  along  the  left  vertical  axis,  the  V  scale  along  the  right  axis 

Figure  93  presents  "4P0"  and  V  as  functions  of  water  depth  for  events 
In  basalt  (Dp  =  30,36,  and  42  ft). 


Corrected  acoustic  calculation  -  airblast  parameters,  A?  and  a  for  basalt 
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Figure  94  gives  "AP0"  and  "a"  as  functions  of  water  depth  for 
events  In  granite  (Df  «  30  and  36  ft). 

Note  that  the  water  depths  "Dw"  and  rock  depths  "Dr"  In  Figures  92-94 
are  for  a  10- ton  AANS  charge  yield  (15.9  tons  energy).  These  depths  may  be 
cube-root  scaled  to  other  yields  of  Interest.  The  water  depth  range  and 
sample  rock  depths  used  here  were  chosen  to  represent  typical  cratering  events 
near  optimum  burial  depth.  Acoustic  calculations  and  alrblast  predictions 
may  also  be  performed  for  other  burial  depth  combinations  of  Interest. 


Tugboat  alrblast  data.  Alrblast  measurements  were  obtained  for  the 
five  Tugboat  Phase  I  single-charge  experiments  In  submerged  coral. Data 
show  certain  deviations  from  acoustic  behavior  (Eqn.  II):  The  peak  overpressure 
attenuation  rates  far  from  surface  ground  zero  are  typically  on  the  order  of 
R  «  Instead  of  R  Rapid  attenuation  Is  observed  when  air  temperature 
decreases  with  increasing  altitude  (sound  rays  are  refracted  away  from  surface 
level,  progressively  decreasing  the  surface  overpressures  at  long  ranges). 

Such  conditions  often  exist  at  seashore  locations.  Reverse  refraction  may 
also  occur  at  higher  altitudes  or  at  an  Inversion  layer,  bending  the  rays 
back  toward  ground  level  and  creating  a  "focus”  (area  of  Increased  overpressure) 
several  thousand  feet  from  ground  zero.  Tugboat  measurements  show  Indications 


of  both  effects:  rapid  attenuation  at  Intermediate  ranges,  with  one  "high"  data 
point  at  long  range  (In  three  cases).  These  deviations  must  be  taken  Into 
account  when  comparing  measurements  with  acoustic  predictions. 


Corrected  acoustic  calculations  have  been  performed  and  alrblast 
predictions  have  been  obtained  for  each  of  the  five  Tugboat  Phase  I  configu¬ 
rations  (coral  +  water  layer).  The  acoustic  results  and  measured  data  are 
summarized  In  Table  14.  This  tabulation  lists  the  approximate  charge  weight 
and  true  depth  of  burial  for  each  event.  Since  four  of  the  events  were  one- 
ton  charges,  the  length  dimensions  (D^.D^,  and  R)  have  been  cube-root  scaled  to  a 
charge  weight  of  10  tons.  Next,  the  acoustic  calculations  for  the  given  rock 
and  water  depth  are  presented:  Water  spall  velocity  "Vspall"  at  SGZ,  "AP0"  In  psl, 
and  a  .  Finally,  the  measured  alrblast  data  have  been  least-squares  fitted 
by  an  equation  of  the  form, 


D  +D 

..  /  r  W\n 
a  (-R— )■ 


Eqn.  IV. 
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where  "AP0"  Is  obtained  from  the  acoustic  predictions;  "a"  and  "n"  are  the 
least-squares  constants.  The  complete  fitted  equations  and  the  "a"  constant 
are  listed  In  the  last  two  columns  of  Table  14.  The  longest-range  "high" 
data  point (which  shows  evidence  of  a  focus )has  been  omitted  from  the  fit  In 
three  cases. 

The  acoustic  alrblast  predictions  listed  In  Table  14  may  be  expressed  as. 


AP 

*1 


a 


D+D 
r  w 

~ir~ 


Eqn.  V. 


aPq  and  a  =  acoustic  prediction  parameters. 

This  form  Is  directly  comparable  with  the  least  squares  fit  given  above 
(Eqn.  IV).  Substituting  the  acoustically-predicted  "AP0"  and  "a"  values 
and  the  layer  depths  ("Dr"  and  "Dw")  for  each  event  Into  Eqn.  V,  we  obtain 
a  predicted  alrblast  function  "AP/AP  "  vs  "(Dr+Dw)/R"  (straight  line  of  slope 
"R-l").  The  acoustic  prediction  R‘T  lines  and  the  least  squares  fits  (dashed 
lines)  for  each  of  the  five  Tugboat  events  are  presented  In  Figures  95-99. 

The  measured  data  points  are  also  shown.  In  this  case,  It  has  been  convenient 
to  plot  "AP/AP0"  against  "R/(Dr+Dw)"  [rather  than  ”(Dr+Dw)/R"];  this  convention 
Is  followed  because  alrblast  overpressures  are  normally  plotted  as  decreasing 
overpressure  vs.  Increasing  range  "R".  It  Is  evident  In  all  figures  that  the 
data  and  the  least  squares  fits  are  lower  and  show  a  faster  attenuation  rate 
than  the  acoustic  "R-l"  predictions.  The  average  attenuation  rate  for  the 
least  squares  lines  Is  about  "R’^*^®"  (constant  n  *  +1.175  In  fitted  equation). 
Rapid  attenuation  Is  attributed  to  refraction  of  the  alrblast  pulse  away  from 
ground  level. 


Due  to  the  differing  attenuation  rates,  Tugboat  measurements  and 
acoustic  predictions  cannot  be  directly  compared  at  long  ranges.  No  data 
points  are  available  at  very  close  ranges,  near  the  high-velocity  spall  mound. 
However,  an  approximate  comparison  may  be  based  on  close-range  predictions 
and  the  expected  behavior  of  the  spall-induced  alrblast  signal.  The  alrblast 
pulse  from  a  buried  cratering  event  originates  within  the  range  Interval 
R/(Dr+Dw)  j  1  to  3  [i.e.,  (Dr+Dw)/R  {  1  to  0.33].  Beyond  this  interval, 


R/(D  +D  ) 
'  r  w 


Figure  95.  Airbl.at  overpressures,  TUGBOAT  IA  event  (in  saturated  coral  overlain  by 


surface  spall  velocities  are  usually  much  less  than  1/2  to  1/10  of  the  SGZ 
velocity,  and  most  of  the  airblast  pulse  energy  may  be  expected  to  have 
decoupled  from  the  ground  (or  water)  surface.  Thus,  it  Is  legitimate  to  expect 
that  refraction  effects  and  faster  attenuation  will  become  apparent  at 
ranges  R/(Dr+Dw)  M  to  3.  Accordingly,  straight  lines  of  slope  "R-1*2" 
have  been  extrapolated  from  the  acoustic  "R-^"  predictions  at  R/(D  +D  )  =  1 
and  at  R/(Dr+Dw)  =  3.  These  lines  encompass  most  of  the  data  points  and 
the  least-squares  fits  within  the  region  of  the  data  for  four  of  the  five 
events  (Figures  95-99).  Event  ID  (Figure  98)  still  lies  below  the  prediction. 
Note  that  ID  also  shows  the  fastest  apparent  attenuation  rate  (R“^  *2549) . 

The  measurements  and  acoustic  predictions  may  also  be  compared  by 
extrapolating  the  least  squares  fits  inward,  toward  the  high-velocity  spall 
region.  The  Intercept  of  the  fitted  line  at  R/(Dr+Dw)  =  1  is  equal  to  the 
constant  "a"  In  Eqn.  IV.. This  Intercept  may  be  compared  with  the  acoustic  Pre¬ 
diction  of  "a"  for  each  event.  The  predicted  a-values  are  listed  In  column  9 
of  Table  14,  while  the  least-squares  intercepts  are  given  In  the  last  column. 
Four  of  the  predicted  values  appear  slightly  higher  than  the  extrapolated  fits 
(Fig.  96-99),  although  the  event  IA  prediction  Is  somewhat  lower  (Fig.  95). 
These  results  may  indicate  that  either  the  acoustically  predicted  "aP0's"  (and 
spall  velocities)  or  the  a-values  are  slightly  too  high. 

As  noted  In  the  previous  section,  photographically-measured  late¬ 
time  velocities  for  the  Tugboat  events  are  somewhat  lower  than  the  corrected 
acoustic  predictions  (peak  spall  velocity).  The  measured  velocities  and  "aP  " 
values  calculated  from  these  velocities  are  listed  in  Table  15.  The  Tugboat 
airblast  data  have  been  reanalyzed,  using  the  new  APQ  values  to  determine 
AP/APg  .  Least-squares  fits  to  the  revised  data  give  the  a-intercept  values 
listed  In  the  last  column  of  Table  15.  The  theoretical  maximum  value  of 
acoustic  "a"  (for  spall  launch  by  a  spherical  non-dissipative  shock  wave) 

Is  equal  to  0.5.  Four  of  the  five  a's  in  Table  15  exceed  0.5  (measured 
surface  velocity  for  the  fifth  event,  Tugboat  IE,  was  In  much  closer  agreement 
with  the  corrected  acoustic  predictions).  The  systematic  deviation  indicates 
that  measured  velocities  and  corresponding  APQ- values  are  too  low.  The 
corrected  acoustic  predictions  are  apparently  more  consistent  with  measured 


Tugboat  Phase  I  Events,  Calculated  Alrblast  Parameters  For  Measured  (Photographic) 
Surface  Velocities 
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airblast  data.  This  result  supports  the  earlier  conclusion  that  measured 
velocities  do  not  represent  the  true  early-time  peak  spall  velocity  at  the 
water  surface. 

It  Is  concluded  that  the  acoustic  technique  can  be  used  to  predict  alrblast 
from  submerged  cratering  events;  however,  allowance  must  be  made  for  the  variations 
In  long-range  alrblast  pulse  attenuation  rate  which  occur  under  some  circumstances. 

X.  SUBSIDIARY  CALCULATIONS 

In  an  effort  to  resolve  some  of  the  uncertainties  and  pursue  Interesting 
results  revealed  by  the  first  suite  of  SOC  problems,  several  additional  SOC 
calculations  were  run.  These  Include:  (a)  Modified  coral  and  granite  under¬ 
water  configurations,  Intended  to  examine  methods  for  better  modeling  shock 
propagation  In  the  water  layer;  and  (b)  coral  and  basalt  "equivalent  overburden" 
calculations  (no  water  layer),  undertaken  to  compare  the  dynamic  effects  of  a 
water  layer  wltn  those  of  an  equivalent  mass  of  rock.  This  section  reviews 
the  subsidiary  calculations  and  compares  results  with  the  previous  set  of 
problems. 

Shock  resolution  In  the  water  layer.  Significant  difficulties  Involving 
excessive  shock  damping,  rapid  attenuation,  and  disagreement  with  the  spherical 
acoustic  calculations  In  the  water  layer  were  encountered  only  for  the  granite 
underwater  events.  It  therefore  seems  probable  that  these  difficulties  are 
related  to  unusual  granite  medium  properties  (stiff,  high  strength),  errors  In  the 
acoustic  analysis,  or  errors  In  the  water  layer  SOC  modeling.  Shock  resolution  In 
the  water  layer  presents  an  Immediate  candidate  as  a  causative  mechanism.  Wave¬ 
forms  have  been  examined  for  the  36/36  ft  (deep  water  layer)  vertical  calculations 
in  coral,  basalt,  and  granite.  The  velocity  waveforms  at  the  rock  Interface 
(R  ■  10.93  m),  in  water  just  beyond  the  rock  interface  (R  *  11,3  id),  and  at 
long  range  In  the  water  layer  (R  =  19.9  m)  are  presented  in  Figures  100  (coral 
36/36),  101  (basalt  36/36),  102  (granite  36/36  -  coarse  water  zones),  and  103 
(granite  36/36  -  fine  water  zones).  The  vertical  (velocity)  scales  In  these 
figures  are  arbitrary;  all  peaks  have  been  plotted  about  the  same  height  to 
permit  comparison  of  rise  times.  The  horizontal  scale  consists  of  an  arbitrary 
10  msec  time  Interval,  and  the  waveforms  for  each  medium  have  been  moved  close 


10-msec  time  interval 


Figure  102.  Velocity  waveforms:  granite  36/36 ,  coarse  water  zone 
(arbitrary  scales). 
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In  granite, 
R  =  10.93  m 


Figure  103 • 


Velocity  waveforms:  granite  36/36,  fine  water  zones 
(arbitrary  scales). 
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together.  It  Is  evident  that  the  moderately  sharp  velocity  peaks  for  coral 
(Figure  100)  and  the  somewhat  rounded  waveforms  of  basalt  (Figure  101) 
undergo  no  dramatic  modification  upon  transmission  at  the  interface, or  in 
travelling  about  9  m  through  the  water  layer.  A  quite  different 
situation  occurs  for  the  sharp,  short-duration  and  rather  noisy  pulses  In 
granite  (Figures  102,  103).  An  abrupt  increase  In  rise  time  and  an  even  more 
noticeable  broadening  of  the  peak  is  apparent  just  beyond  the  water  Inter¬ 
face.  More  gradual  changes  continue  throughout  the  water  layer.  Waveforms  near 
the  surface  (R  =  19.93  m)  resemble  the  broad,  smooth  peaks  for  coral  and 
basalt  calculations  (note  the  slight  differences  in  time  scales  between  figures). 
Smearing  of  the  front  and  rise  time  Increase  are  more  noticeable  for  the  "coarse 
water  zone"  granite  calculation,  but  occur  even  in  the  "fine  zone"  case.  The 
"fine  zone"  waveforms  (half-size  water  zones)  show  evidence  of  noise  trans¬ 
mission  and  "ringing"  on  t.ie  declining  edge  of  the  water-layer  waveforms 
(Figure  103).  This  effect  is  also  visible  In  the  waveforms  shown  In  Appendix  F 
(Figures  F7-F8). 

A  more  quantitative  comparison  of  rock  and  water  layer  waveforms  for 
the  36/36  ft  calculations  Is  given  In  Table  16;  this  tabulation  presents 
the  approximate  velocity  peak  rise  times  in  the  top  part  of  the  rock  layer 
and  in  the  water  layer,  and  the  pressure  or  velocity  peak  half-widths*  near  the 
rock  Interface  and  in  the  water  layer.  The  quantities  are  listed  for  coral 
36/36,  basalt  36/36,  and  granite  36/36  (coarse  and  fine  zone)  calculations, 
and  also  for  additional  calculations  "coral  36/36"  witn  very  low  damping  in 
the  water  layer, and  "modified  granite  36/36"  with  equal  size  rock  and  water 
zones  of  width  *  0.099  m  (discussed  in  greater  detail  below).  The  half-widths 
of  the  rapidly-declining  pressure  peaks  were  measured  for  coral  and  basalt; 
however,  the  "granite"  pressure  peaks  were  so  rapidly  damped  and  noisy  that 
velocity  peak  half-widths  were  used  instead.  It  is  most  useful  to  compare 


*  Pulse  width  in  msec  at  one-half  of  peak  anplitude. 
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the  relative  trends  In  rise  time  and  peak  half-width  for  the  three  media. 

The  rise  times  undergo  no  significant  modification  very  close  to  the  rock- 
water  Interface  for  coral  (1.0  msec)  and  basalt  (1.2  msec),  and  only  a 
gradual  change  while  traversing  the  thick  water  layer.  Indeed,  the  wave¬ 
form  parameters  remain  almost  constant  for  the  sample  "coral  36/36"  low- 
damping  calculation.  The  pressure  pulse  half-widths  likewise  show  only  moderate 
changes  near  the  Interface  and  In  the  water  layer.  It  is  Interesting  to 
note  that  the  factor  by  which  the  pressure  pulse  half-width  increases  In  the 

water  layer  Is  roughly  comparable  for  basalt  and  coral,  in  spite  of 
the  enormous  water  zone  size  used  In  the  "basalt"  calculations.  This  fact 
Is  attributed  to  the  relatively  smooth  and  well -resolved  waveforms  entering 
the  water  layer.  It  Indicates  that  excellent  resolution  was  achieved  in  the 
water  layer  for  these  events  and  explains  the  close  agreement  with  acoustic 
calculations. 

The  granite  "coarse"  and  "fine"  zone  calculations  (Table  16)  reveal 
a  dramatic  Increase  In  rise  time  close  to  the  rock-water  Interface.  This 
change  occurs  over  a  very  few  calculatlonal  zones  (the  exact  distance  cannot 
be  defined  due  to  Inadequate  resolution  In  the  data  edits)  and  is  most 
severe  for  the  coarse-zcned  calculation.  Its  immediate  effect  is  to  smear 
the  shock  front  and  severely  decrease  the  peak.  The  consequent  sudden  decline 
of  particle  velocities  near  the  Interface  was  discussed  in  Chapter  VI 
(Figures  A38,  A39,  A46,  A47),  and  obviously  results  from  an  abrupt  loss  of 
resolution.  In  order  to  reduce  this  effect,  a  modified  granite  36/36  calcu- 
latlon  was  performed  with  somewhat  coarser  rock  zones  and  with  equal-size 
water  zones.  This  calculation  eliminated  the  sudden  Increase  In  rise  time 
near  the  Interface,  and  the  peak  velocities  likewise  ceased  to  decrease 
abruptly  at  the  Interface.  However,  the  general  trend  of  increasing  rise 
time  and  rapid  attenuation  throughout  the  entire  water  layer  persisted. 

Rise  time  within  the  water  layer  increased  by  a  factor  of  2.5,  In  spite  of 
the  elimination  of  the  discontinuous  behavior  near  the  Interface.  The 
attenuation  rate  of  peak  velocity  also  remained  close  to  R"^'5  in  the  water 
layer,  comparable  with  the  R’^’^  factor  for  the  other  granite  events  (beyond 
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the  immediate  vicinity  of  the  Interface).  In  summary,  rise  times  for  all 
three  of  the  granite  calculations  tend  to  Increase  from  a  few  tenths  of  a 
millisecond  to  about  1  or  2  msec  within  the  water  layer,  and  correspondingly 
rapid  attenuation  of  the  peaks  Is  observed.  This  trend  Is  not  noticed  for  the 
coral  and  basalt  calculations,  which  fall  In  close  agreement  with  predictions 
for  compresslonal  acoustic  wave  transmission.  Possible  reasons  for  the 
differing  behavior  are  examined  below. 

Modified  coral  36/36  calculation,  low  damping  in  the  water  layer. 

Good  shock  resolution  and  near-acoustic  behavior  in  the  water  layer  were 
achieved  for  the  coral  and  basalt  underwater  events.  In  order  to  more  pre¬ 
cisely  define  the  factors  which  Influence  the  best  multilayer  calculations, 
further  SOC  runs  were  undertaken  for  the  coral  36/36  configuration.  The 
problem  format  and  the  Impedance-matched  rock  and  water  zoning  were  Identical, 
but  very  low  water  viscosity  factors  v  z  0.05  to  0.07  were  used.  Viscosity 
factors  lower  than  s  0.1  normally  cause  peak  overshoot  and  severe  ringing 
behind  the  shock  front.  In  this  case,  it  was  noted  that  noise  from  the  coral 
layer  was  transmitted  undamped  Into  the  water,  but  remarkably  little  evidence 
of  overshoot  or  ringing  was  apparent  (relative  to  the  "normal"  coral  36/36 
calculation).  Noticeable  ringing  began  to  develop  several  meters  beyond 
the  rock-water  interface  when  the  viscosity  factor  was  reduced  to  0.05.  The 
success  of  these  calculations  is  attributed  to  the  gradual  rise  times  and 
smooth  shock  contours  near  the  coral -water  interface.  At  v  =  0.05,  at 
water  layer  velocity  attenuation  rate  of  about  R"1,03  was  achieved  (compared 
with  spherical  nondlssipatlve  attenuation  rate  of  R”1*0,  and  with  the  measured 
rate  for  water  and  earlier  SOC-calculated  rates  ~R”1,13  to  R-1,2).  For 
problems  of  this  type,  shock  attenuation  in  the  water  layer  may  be  adjusted 
and  nondissipati ve  acoustic  behavior  may  be  simulated  with  no  serious  adverse 
effect  on  the  solution.  Modifying  the  water  layer  viscosity  did  not  signifi¬ 
cantly  Influence  the  Initial  reflection  at  the  rock-water  Interface  nor  did 
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1t  change  the  calculated  dynamic  behavior  of  the  rock  layer.  Spall 
velocities  near  the  rock  Interface  were  60-70  m/sec,  similar  to  the  previous 
SOC  calculation  and  to  the  acoustically-calculated  interface  peak  spall  velocity  of 
71.25  m/sec.  However,  shock  attenuation  In  the  water  layer  was  decreased, 
and  the  water  surface  spall  velocity  Increased  correspondingly.  Water 
surface  layer  spall  velocities  were  38-70  m/sec,  compared  with  34-62  m/sec 
for  the  earlier  coral  36/36  calculation  and  with  a  peak  water  surface  spall 
velocity  of  65.1  m/sec  for  the  spherical  (uncorrected)  R“^3  acoustic  calcu¬ 
lation  (Table  11).  A  revised  acoustic  calculation  (uncorrected)  was  per¬ 
formed  with  a  "water"  attenuation  rate  of  R"1*03,  to  match  the  lower  rate 
attained  In  this  problem.  The  acoustic  peak  water  surface  spall  velocity 

was  found  to  be  69.8  m/sec.  In  excellent  agreement  with  the  70  m/sec  SOC 
result. 

Additional  calculations  were  attempted  with  still  lower  water  viscosity 
factors,  but  ringing  and  overshoot  were  encountered  at  ranges  close  to  the 
Interface.  The  use  of  lower  viscosity  factors  Is  not  justified,  since 
essentially  acoustic  attenuation  can  be  achieved  In  this  case.  It  Is 
concluded  that,  for  smoothly-rising  "resolvable"  shock  waveforms  typical 
of  coral  and  basalt,  the  viscosity  factor  controls  attenuation  In  the  water 
layer.  Loss  of  resolution  due  to  impedance  matching  and  water  zone  size  Is 
not  a  significant  effect.  Ho  serious  discontinuity  in  shock  behavior  occurs 
at  the  Interface.  These  conclusions  are  consistent  with  the  earlier  dis¬ 
cussion  of  rise  times  and  peak  widths.  The  viscosity  factor  may  be  used 
to  control  the  water  layer  calculation  and  achieve  any  desired  attenuation 
rate  for  underwater  calculations  in  coral-  and  basalt- type  media. 

"Standard"  and  modified  granite  underwater  calculations  -  zone  size 
effects.  Granite  underwater  events  present  the  opposite  situation  of  very 
sharply-rising,  noisy  waveforms  at  the  interface  and  resolution-limited 
propagation  within  the  water.  Attenuation  rates  greatly  exceed  the 
acoustic  level  throughout  the  water  layer.  Low  water  viscosity  factors 
cannot  be  used  to  reduce  attenuation  because  of  the  sharply  rising  and 
noisy  shock  input  from  the  rock  layer.  Likewise,  the  use  of  very  fine  water  zones 


creates  a  severe  acoustic  mismatch  at  the  Interface,  which  also  causes 
ringing  for  sharp-fronted  shock  waves.  The  succeeding  paragraphs  discuss 
in  more  detail  some  of  the  results  of  granite  underwater  calculations. 

The  initial  granite  36/10  and  36/36  calculations  were  conducted  with 
the  standard  rock  zone  size  of  0.04  m  and  impedance-matched  water  zones 
of  about  0.39  m  (Chap. Ill,  Table  2).  These  coarse  water  zone  configur¬ 
ations  produced  a  sudden  drop  In  particle  velocity  and  a  lack  of  pressure 
restoral  beyond  the  rock-water  Interface  (just  above  the  "low-pressure" 
rock  zones  in  the  spall  cutoff  region  below  the  Interface).  A  corresponding 
loss  of  shock  resolution  and  sudden  peak-smearing  were  noted  near  the  Inter¬ 
face.  At  slightly  longer  ranges,  the  particle  velocities  approached  a 
constant  but  rapid  attenuation  rate  of  about  R-1*55,  and  gradual  smearing 
of  the  wavefront  continued  throughout  the  10  ft  and  36  ft  water  layers. 

The  two  granite  calculations  were  repeated  using  halved  water  zone  size 
(*  0.196  m),  with  the  result  that  the  sudden  velocity  dropoff  and  the  degree 
of  smearing  of  the  front  near  the  Interface  were  reduced  significantly. 

However,  these  fine-zone  calculations  continued  to  show  gradual  loss  of 
resolution  and  a  rapid  attenuation  rate  In  the  water  layer.  The  fine-zoned 
calculations  also  began  to  reveal  impedance-mismatch  effects  In  the  water 
layer.  Impedance  mismatch  may  be  easily  distinguished  from  artificial 
viscosity  effects,  since  mismatch  tends  to  cause  erratic  ringing  on  the 
declining  edge  of  the  waveform  close  to  the  Interface,  while  Inadequate 
damping  gives  rise  to  gradually-developing  overshoot  and  oscillations  at 
longer  ranges  (for  a  low-amplitude  wave).  Water  waveforms  near  the  Inter¬ 
face  (R  =  11.4  m)  for  granite  36/36  are  shown  in  Figures  104  (coarse  zones) 
and  105  (fine  zones).  Artificial  viscosity  factor  for  the  water  layer  was 
v  =  0.1  in  both  cases.  The  fine-zoned  calculation  gives  a  higher  pressure 
and  better  peak  resolution,  but  some  erratic  ringing  is  apparent  In  the  late¬ 
time  waveform  (impedance  mismatch  =  factor  of  2). 

A  closer  examination  of  the  coarse  and  fine  zone  calculations  demonstrated 
that  the  major  differences  were  primarily  confined  to  the  water  layer  and  had 
little  influence  on  the  underlying  rock.  Waveforms  In  the  rock  layer  were  almost 
identical  except  for  slight  changes  In  tensile  reflection  arrival  time  and  duration 


Granite  36/36,  pressure  waveform  in  water  layer  at 
R  =  11.4  m;  v  =  0.1,  fine  zones. 
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of  the  f reef all  phase,  and  In  the  late-time  residual  velocities.  The  beginning  of 
the  "water  overburden"  deceleration  (recompactlon)  phase  well  below  the  rock-water 
Interface  and  the  subsequent  arrival  of  the  upward-moving  gas  acceleration  pulse 
were  slightly  earlier  for  the  "fine"  zone  calculations.  The  overall  kinetic 
and  Internal  energies  In  the  rock  layer  were  modified  by  less  than  3%, 
and  did  not  diverge  at  late  times.  However,  the  water  layer  velocities 
and  energies  began  to  differ  noticeably  as  the  water  surface  spall  reflection 
moved  downward  through  the  water  layer.  Kinetic  energies  Increased  by  about 
4%  and  Internal  energies  decreased  by  up  to  10%  at  late  times  for  the  fine- 
zoned  calculations,  due  to  better  shock  resolution  and  higher  spall  velocities 
In  the  water. 

The  granite  36/36  calculation  was  repeated  in  modified  form  with 
equal  rock  and  water  zone  sizes  of  0.099  m.  Shock  rise  time  in  the  rock 
layer  was  slightly  increased  due  to  the  coarser  rock  zones,  but  the  peak 
values  and  shock  width  near  R  s  10  m  (close  to  the  rock-water  interface) 
were  only  negligibly  changed.  This  fact  Indicates  that  somewhat  coarser 
rock  zoning  could  have  been  used  with  no  deleterious  effect  in  the  vertical 
calculations. 

The  equal  zone  sizes  and  finer  water  zones  exerted  a  considerable 
Influence  on  the  water  layer  calculations,  as  discussed  previously.  The 
abrupt  loss  of  peak  resolution  at  the  rock-water  Interface  and  the  corres¬ 
ponding  sudden  decline  of  velocity  in  the  water  layer  just  beyond  the  interface 
both  disappeared.  However,  the  slowly-increasing  rise  time  (to  about  1  msec) 
and  rapid  overall  attenuation  In  the  water  layer  persisted.  The  use  of  equal 
water  zones  also  caused  almost  undamped  transmission  of  "granite"  noise 
throughout  the  thick  water  layer,  with  the  addition  of  severe  ringing  near 
the  interface  due  to  the  gross  Impedance  mismatch. 

The  elimination  of  discontinuous  behavior  near  the  interface  and 
improvement  of  water-layer  resolution  slightly  affected  other  aspects  of  the  cal¬ 
culations.  The  previous  granite  underwater  events  (Table  11)  gave  a  rock 
interface  peak  spall  velocity  of  about  32  m/sec,  approximately  equal  to  the 
peak  spall  velocity  with  no  water  overburden  (granite  36/00).  The  peak 
velocity  was  thus  not  reduced  as  expected  for  imperfect  reflection  at  the 
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gran  <te-wa  tor  interface.  The  early-time  effects  of  the  water  layer  were  not 
properly  calculated  due  to  the  large  water  zone  size  (zone  width  greater  than 
the  effective  wavelength  of  the  incoming  shock).  The  modified  granite  36/36 
calculation  with  equal  zoning  produced  a  velocity  range  of  18  to  30  m/sec  in  the 
rock  surface  spall  layer.  The  peak  velocity  of  30  m/sec  lies  close  to  the  acoustic 
value  of  29-29.6  m/sec  for  reflection  at  a  granite-water  Interface  (Table  11). 
Hence,  the  modified  SOC  calculation  Is  seen  to  be  more  consistent  with  expected 
spall  velocity.  The  difference  In  spall  velocities  between  the  SOC  calculations 
declines  at  locations  further  below  the  Interface;  even  at  the  Interface,  differ¬ 
ences  In  peak  velocity  disappear  soon  after  shock  arrival,  as  the  peak  Is  damped 
by  water-layer  deceleration  of  the  rock  material. 

The  modified  water  layer  calculation  also  Influences  water  spall 
velocities.  The  surface  spall  layer  is  somewhat  thinner  and  more  clearly 
defined  due  to  Improved  resolution  of  the  shock  front.  Surface  layer  spall 
velocities  for  the  modified  granite  36/36  configuration  are  11.5  to  20  m/sec, 
significantly  higher  than  the  previous  coarse-zoned  (9.7-16.1  m/sec)  and 
f1ne-zoned(10.6-18  m/sec)  calculations.  However,  the  peak  spall  velocity 
remains  lower  than  the  spherical  (uncorrected)  acoustic  value  of  26.5  m/sec 
(Table  11,  Figure  79;  assuming  R  ‘  0  attenuation).  Substituting  a  water- 
layer  attenuation  rate  of  R  into  the  spherical  acoustic  analysis 
gives  a  spall  velocity  of  20.0  m/sec  (Figure  79),  In  precise  agreement  with 
the  modified  SOC  result. 

We  see  that  the  apparent  disagreements  between  "SOC"  and  "acoustic" 
calculations  for  granite  are  attributable  to  discontinuous  behavior  near 
the  Interface  (eliminated  In  the  modified  calculation)  and  rapid  attenuation 
through  the  water  layer  (not  eliminated).  Rapid  attenuation  Is  not  solely 
attributable  tc  water  zone  sizes,  since  near-acoustic  behavior  was  attained 
In  coral  calculations  (water  zones  =0.1  m)  and  even  in  the  basalt  calculations 
with  coarse  water  zoning  (=0.254  m).  Granite  events  consistently  show  rapid 
attenuation  well  beyond  the  rock-water  Interface,  even  when  equal-size  rock 
and  ^a ter  zones  are  used.  This  may  represent  a  physical  inconsistency  In  the 
results  or  a  real  trend  dependent  on  material  properties.  The  underlying 


granite  is  a  stiff  high-velocity  rock  which  transmits  sharp  waveforms.  It 
possesses  a  high  shear  strength  and  a  low  Poisson's  ratio.  Thus,  the  material 
supports  large  shear  stresses  which  build  up  even  at  close  ranges,  and  shows  a 
high  shear-wave  velocity.  The  radial  and  tangential  stress  components  differ 
considerably  at  ranges  £10n,  near  the  Interface  range  (see,  for  example,  peak 
stresses  plotted  In  Figure  A15,  Appendix  A).  This  situation  differs  from  the 
coral  and  basalt  media,  for  which  the  low  ultimate  strength  severely  limits  shear 
stress  and  behavior  Is  more  nearly  fluid.  The  acoustic  analysis  used  here  takes 
account  only  of  peak  compresslonal  wave  transmission  and  refracvlon  at  an 
Interface.  Shear  wave  effects  and  shear  wave  reflection  at  the  Interface  are 
not  considered.  The  errors  thus  Incurred  may  Influence  peak  attenuation  beyond 
the  Interface  for  the  stiff,  shear-supporting  granite  material.  The  situation 
Is  further  complicated  In  this  case  because  of  the  failure  behavior  and 
spallation  taking  place  near  the  Interface.  The  granite  rapidly  undergoes 
shear  failure  at  a  pressure  and  stress  level  below  Its  ultimate  strength, 
falls  In  tension,  spalls,  and  then  Is  reloaded  In  shear  by  water  layer 
recompaction.  The  behavior  at  and  Immediately  after  shock  arrival  may  thus 
be  quite  complex.  Rapid  attenuation  encountered  In  other  multilayer  one¬ 
dimensional  calculations^8^  has  been  attributed  to  unloading  and  energy 
feedback  Into  the  low-pressure  reflected  rarefaction  region  beneath  the 
Interface.  This  effect  could  become  Important  for  granite,  a  strong  material 
with  a  very  narrow  high-velocity  spall  zone  at  the  Interface  and  rapid 
velocity  damping  beneath  the  Interface.  Two-dimensional  calculations  may  be 
required  to  solve  the  general  layer  transmission  problem  for  granite  and 
similar  stiff,  high-strength  materials.  Part  of  the  apparent  inconsistency 
In  attenuation  Is  undoubtedly  due  to  the  shear  reflection  effects,  which  can 
be  studied  by  a  more  comprehensive  acoustic  treatment.  Part  may  be  due  to 
failure  and  other  nonlinear  material  behavior  which  can  be  properly  analyzed 
only  by  complete  calculations. 

To  sunmarlze,  a  judicious  choice  of  problem  format  and  zone  size 
producer  valid  two-layer  calculations  for  the  coral  and  basalt  submerged 
cratering  cor  guratlons  examined  here.  Peak  pressures  and  velocities  for  the 
water  layer  end  other  safety-related  parameters  such  as  alrblast  can  be 
predicted  by  a  simple  acoustic  analysis.  The  granite  medium  imposes  rather 
more  severe  calculatlonal  restrictions  due  to  the  large  density  and  sonic 
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velocity  mismatch  at  the  Interface,  and  due  to  the  shear  supporting  character¬ 
istics  and  other  properties  of  the  material.  Fortunately,  the  precise  behavior 
In  the  water  layer  exerts  very  little  Influence  on  the  underlying  rock. 
Lagranglan  codes  with  coarse  water  zoning  may  be  confidently  used  to  predict 
shock  transmission,  long-term  dynamics,  and  cratering  for  the  rock  layer. 

Equivalent  overburden  calculations.  Attests  have  been  made  to 
treat  the  multilayer  cratering  problem  with  simple  "overburden"  models. 

Possible  approaches  Include  the  use  of  a  single  layer  of  the  same  total 
depth  and  average  density  (weighted  by  layer  thickness)  as  the  series  of 
layers  In  the  actual  configuration.  The  two-layer  underwater  case  has 
been  analyzed  by  assuming  that  the  water  layer  acts  as  a  depth  of  rock 
equivalent  to  the  overburden  weight  (llthostatlc  head)  of  the  water  layer  [12]; 
[Equivalent  "rock"  depth  of  water  layer  <=  Dw  (p*/ pr);  total  equivalent  burial 
depth  of  shot  Dr  +  Dw  (pw/pr)]*  Dimensional  arguments  may  be  cited  In  sup¬ 
port  of  this  approach.  The  complexity  of  multilayer  shock  Interactions 
provides  certain  counter-arguments.  Consider  the  cases  of  an  undervater 
cratering  event  and  of  Its  "equivalent  overburden"  analog  In  rock  alone.* 

The  underwater  event  experiences  an  Initial  rock  layer  spall  phase  caused 
by  partial  tensile  reflection  at  the  rock-water  Interface.  This  Is  followed 
by  a  delayed  spall  pulse  from  the  water  surface,  transmitted  through  the 
water  layer.  The  gas  acceleration  phase  begins  at  relatively  early  time, 
when  the  Interface  spall  reflection  reaches  the  cavity.  In  the11  equivalent 
overburden" case,  the  Initial  shock  wave  Is  perfectly  reflected  as  a  single 
tensile  wave  at  the  rock  free  surface.  Due  to  the  added  rock  overburden, 
this  pulse  arrives  at  a  given  "rock"  location  relatively  later  than  the 
corresponding  rock-water  Interface  reflection,  but  much  earlier  than  the 
corresponding  water  surface  reflection  (which  must  travel  through  the  thicker, 
low-velocity  water  layer).  The  gas  acceleration  phase  also  begins  somewhat 
later  because  of  the  Increased  rock  depth.  Comparison  is  further  complicated 
by  the  fact  that  the  underwater  event  experiences  early-time  spall  rarefaction 


"THe^equIvalent  overburden"  event  has  additional  rock  overburden  thinner 
thickness**^  (J^/p  )Ut  °f  CqUa1  overburden  w^9ht;  added  rock  overburden 
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of  the  rock  layer,  followed  successively  by  recompaction  due  to  water 
overburden  and  relief  from  the  water  surface  reflection.  These  depth 
dependent  interactions  begin  Immediately  after  Interface  spall  reflection, 
and  affect  the  later  dynamic  history  and  transmission  of  all  subsequent 
pulses  In  the  rock  layer.  The  Interactions  are  not  Identical  to  the  delayed 
single  spall  pulse  and  the  more  continuous  damping  which  occur  for  the  equivalent 
overburden  event  in  rock  alone.  Finally,  the  water  layer  and  rock  over¬ 
burden  problems  differ  with  regard  to  strength  effects,  shock  attenuation, 
and  spall  velocities  In  the  upper  (water  or  rock)  region.  Even  Ignoring 
these  Inherently  different  aspects  of  solid  medium  behavior,  such  as 
ability  to  support  shear  and  the  material  strength.  It  seems  Improbable 
that  the  complex  double  pulse  spall  and  recompaction  phenomena  for  the  under¬ 
water  event  would  produce  exactly  the  same  dynamic  results  as  the  single¬ 
reflection  overburden  case.  One  might  expect  the  greatest  discrepancies 
to  occur  for  thick  water  layers  (l.e.,  largest  difference  between  water 
layer  and  rock  overburden  thicknesses,  greatest  difference  between  reflec¬ 
tion  arrival  times,  severe  recompaction  effects  near  the  rock-water  Inter¬ 
face  for  the  underwater  event),  and  for  large  Impedance  mismatch  at  the' rock- 
water  Interface  (largest  differences  between  depths  and  sonic  velocities  of 
the  water  and  rock  overburden  layers,  greatest  difference  between  pulse 
arrival  times,  longest  relative  delay  between  first  and  second  spall 
pulses  for  the  underwater  event). 

A  set  of  three  "equivalent  overburden"  SOC  calculations  has  been 
Implemented  to  examine  the  dynamic  effects  of  rock  and  water  overburden. 

These  include  coral  41,31  (rock  depth  Df  =  41.31  ft,  simulates  coral  36/10), 

-ora1  55,11  (°r  =  55*11  ft*  simulates  coral  36/36),  and  basalt  39.83 

(Dr  -  39.83  ft,  simulates  basalt  36/10).  In  all  cases,  the  added  rock  depth 

was  chosen  to  match  the  overburden  depth  of  the  corresponding  water  layer. 

The  "equivalent  overburden"  rock  calculations  used  constant  rock  zoning  through¬ 
out  (homogeneous  rock  material,  no  Impedance  mismatch);  however,  the  rock  was 
arbitrarily  divided  into  two  sections  for  purposes  of  analysis.  These  sections 
encompass  the  "lower"  36  ft  and  the  "upper"  overburden  region,  corresponding 
with  the  "rock"  and  "water"  layers  of  the  underwater  calculations. 


In  comparing  spall  and  late-time  residual  velocities  of  the  various 
calculations,  the  following  trends  are  noted: 

For  the  "shallow"  (36/10  and  equivalent)  configurations,  the  water 
layer  spall  velocities  are  generally  higher  than  the  "upper"  region  rock 
spall  velocities  of  the  equivalent  overburden  events.  This  effect  Is 
expected  from  the  lower  Impedance  of  the  water  medium.  More  Important  to 
cratering  dynamics,  the  residual  velocities  near  R=36  ft  (near  the  Interface) 
are  relatively  lower  for  the  underwater  events  (due  to  water- 
layer  damping  of  the  Initial  peak  spall  velocity  near  the  Interface,  and 
due  to  the  nearness  of  the  rock  free  surface  for  the  equivalent  overburden 
events);  the  residual  velocities  then  become  relatively  higher  for  the  under¬ 
water  events  further  below  the  Interface.  Velocity  waveforms  are  Illustrated 
In  Figure  106a  (basalt  36/10,  and  equivalent  overburden  calculation  basalt 
39.83).  The  residual  velocity  after  water  recompaction  near  the  Interface 
(R  *  10.6  m)  for  basalt  36/10  Is  dramatically  lower  than  the  spall  velocity 
at  the  same  location  for  basalt  39.83.  The  water  layer  damping  effects  are 
quite  evident.  At  R  =  9.06  m,  well  below  the  rock-water  Interface,  this 
trend  is  reversed;  the  basalt  36/10  late-time  residual  velocity  (after  slight 
water  damping)  Is  higher  than  the  "equi valent"overburden"  spall  velocity 
for  basalt  39.83.  This  relationship  persists  throughout  the  lower  regions 
of  the  rock  layer. 

Somewhat  different  behavior  is  noted  in  the  "deep"  calculations  (coral 
36/36  and  coral  55.11).  In  this  case,  the  lower  portion  of  the  rock  layer 
for  coral  55.11  is  far  removed  from  the  free  surface.  Initial  peak  velocities 
within  the  lower  36  ft  of  the  rock  are  strongly  damped  before  the  much- 
delayed  free-surface  spall  reflection  arrives.  The  spall  velocities  are 
accordingly  very  low.  The  earlier-arriving  interface  spall  reflection  for 
coral  36/36  proves  more  efficient  in  ejecting  material,  despite  water  layer 
damping  effects.  Velocity  waveforms  near  the  interface  (R  =  10.6  m  and  9.01  m) 
are  presented  in  Figure  106b.  In  both  cases,  the  late-time  residual  velocities 
for  coral  36/36  (t  =  10  to  20  msec)  are  seen  to  be  significantly  higher  than 
the  corresponding  single-layer  rock  spall  velocities.  Evidently,  the  underwater 
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Basalt  3o/lO  and  equivalent  overburden  calculation 
velocity  waveforms  at  R  =  10.6  m  and  R  =  9-09  m. 
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Figure  106b.  Coral  36/36  and  equivalent  overburden  calculation 
velocity  waveforms  at  R  =  10.6  m  and  R  =  9-01  a. 
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36/36  event  produces  more  efficient  material  ejection  within  the  critical 
rock  region  (R  <  36  ft  or  10.97  m)  than  the  deeply-burled  55.11  ft 
"equivalent  overburden"  event.  The  difference  with  regard  to  cratering 
effectiveness  will  be  still  greater,  since  there  Is  no  additional  rock 
overburden  above  36  ft  to  be  ejected  for  coral  36/36.  After  20  msec, 
another  modification  becomes  apparent:  The  early-arriving  gas  acceleration 
pulse  for  coral  36/36  reaches  the  upper  regions  of  the  rock  layer,  increas¬ 
ing  the  velocities  still  further  above  "coral  55.11".  Due  to  the  short  distance 
between  the  cavity  and  the  Interface,  gas  acceleration  occurs  sooner  and  Is 
expected  to  be  a  more  effective  ejection  mechanism  for  coral  36/36.  To 
summarize,  the  underwater  configuration  has  less  rock  material  to  eject 
and  the  velocities  within  the  critical  region  are  higher  than  for  the 
"equivalent  overburden"  event.  The  discrepancy  Is  greatest  for  "deep" 
overburden  layers  ( 1  .e . ,  36/36  ft  comparison).  The  coral  36/36  ft  configur¬ 
ation  thus  simulates  a  more  shallowly  burled  "rock"  event,  and  Is  effectively 
closer  to  optimum  depth  of  burial  than  Its  55.11  ft  analog.  The  "equivalent 
overburden"  analogy  Is  conservative,  tending  to  underestimate  cratering 
efficiency  and  crater  dimensions  for  submerged  events  relative  to  their  dry¬ 
land  counterparts. 

The  overall  mound  dynamic  parameters  have  also  been  examined  for  the 
equivalent  overburden  calculations.  Residual  velocities  at  late  time  are 
plotted  In  igures  107  (coral)  and  108  (basalt).  Other  events  in  the 
respective  media  are  shown  for  comparison.  The  coral  41.31  velocities  are 
much  less  than  coral  36/10  in  the  lower  parts  of  the  mound,  becoming  greater 
within  2  m  of  the  36  ft  (10.97  m)  location  of  the  coral  36/10  inter¬ 
face.  The  coral  55.11  velocities  are  lower  than  coral  36/36  throughout  the 
"36  foot"  region,  and  are  significantly  less  than  50  m/sec  for  most  of  the 
mound  (even  within  the  deeper  "high-velocity"  regions).  These  velocities 
were  determined  after  the  Initial  spallation  phase  was  complete.  Coral  55.11 
Is  obviously  below  optimum  cratering  depth.  The  basalt  39.83  spall  velocities 
(Figure  108)  are  somewhat  less  than  basalt  36/10  in  the  lower  mound,  but  rise 


Residual  velocities  for  coral  36/OO,  coral  36/10,  coral  36/36,  and  equivalent  over¬ 
burden  calculations  coral  41.31  (simulates  36/10)  and  coral  55.11  (simulates  36/36) 


burden  calculation  basalt  39*83  (simulates  36/10). 
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rapldly  within  1.5  m  of  the  36  ft  location.  This  behavior  Is  attributable 
to  the  free  surface  at  39.83  ft  (located  only  slightly  beyond  the  36  ft 
range).  The  overall  velocity  profiles  within  the  36  ft  layer  for 
basalt  36/10  and  39.83  are  much  more  similar  than  the  corresponding  profiles 
for  the  coral  events.  Thus,  the  basalt  "equivalent  overburden"  calculation 
(with  large  Impedance  mismatch)  did  a  better  job  of  modeling  the  underwater 
event  than  the  coral  calculation  (smaller  Impedance  mismatch).  The  conclu¬ 
sion  applies  only  for  these  two  specific  media.  A  much  larger  discrepancy 
would  be  expected  between  basalt  36/36  (deep)  event  and  Its  equivalent  over¬ 
burden  calculation.  A  "deep"  overburden  calculation  for  basalt  would  provide 
useful  Information,  but  has  not  yet  been  performed  due  to  the  long  calculation 
time  required. 

Cavity  radii  and  cavity  pressures  for  all  of  the  equivalent  overburden 
calculations  are  shown  In  Figures  109  and  110,  respectively.  These  para¬ 
meters  give  a  good  Idea  of  the  relative  gas  acceleration  effects  for  the 
"overburden"  and"underwater"  events.  In  all  cases,  the  underwater  calculations 
reveal  progressively  greater  cavity  radii  and  lower  cavity  pressures  than  the 
equivalent  overburden  events  at  late  time.  Energy  transfer  by  gas  acceleration 
Is  considerably  more  effective  for  the  underwater  calculations,  and  the  closer 
location  of  the  rock-water  interface  obviously  has  a  significant  Influence 
on  this  phase  of  cratering  dynamics.  Note  that  the  radii  and  pressures  are 
virtually  identical  for  coral  36/10  and  36/36;  the  Interface  reflection 
completely  controls  cavity  growth  for  these  events.  The  discrepancy  between 
underwater  and  overburden  calculations  Is  rather  large  for  coral  36/36-coral 
55.11,  somewhat  less  for  coral  36/10-coral  41.31,  and  least  for  basalt  36/10- 
basalt  39.83  (time  >  10-15  msec).  Once  more,  the  basalt  overburden  calcu¬ 
lation  fell  closest  to  the  underwater  event,  and  the  basalt  curves  do  not 
diverge  dramatically  at  late  times. 

The  calculated  energies  as  a  function  of  time  for  equivalent  overburden 
and  underwater  calculations  are  depicted  in  Figures  llla-b  (coral)  and  112 
(basalt).  To  permit  meaningful  comparison,  the  problems  have  been  divided 
Into  a  "lower"  region  (below  36  ft),  and  an  "upper"  water  or  overburden  region 


Cavity  radius 


Cavity 


Kinetic  energy  (Mbar 


-40H- 


Figure  111a.  Kinetic  energy  in  the  lower  layer,  coral  underwater 
calculations  (coral  36/IO  and  3°/36 )>  anc*  comparable 
coral  calculations  with  the  same  amount  of  overburden 
(no  water) . 
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Figure  112.  Kinetic,  internal,  and  cavity  energy;  basalt  underwater 

(basalt  36/10),  and  comparable  basalt  calculation  with  the 
same  amount  of  overburden  (no  water). 
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( underwater  or  equivalent  overburden  events).  Coral  kinetic  energies  in 
the  lower  region  or  rock  layer  are  shown  In  Figure  111a.  This  kinetic  energy 
plot  clearly  reveals  the  overall  differences  In  relative  mound  velocities 
for  the  various  events.  Kinetic  energies  for  the  underwater  events  lnwe- 
dlately  deviate  above  the  equivalent  overburden  calculations  when  reflection 
occurs  at  the  rock-water  Interface,  and  remain  appreciably  higher  at  late 
times.  The  discrepancy  Is  greatest  for  the  "deep"  layer  events.  The  coral 
36/10  curve  lies  somewhat  above  and  almost  parallel  to  coral  41.31,  while 
coral  36/36  lies  far  above  coral  55.11.  The  relative  loss  of  kinetic  energy 
to  the  lower  36  ft  region  Is  very  much  less  between  coral  36/10  and  coral 
36/36  than  between  coral  41.31  and  coral  55.11.  The  addition  of  a  thick 
water  overburden  layer  therefore  has  less  effect  on  the  lower  problem  region 
than  the  addition  of  an  equivalent  depth  of  rock  overburden.  It  Is  signifi¬ 
cant  to  note  that  the  coral  36/36  rock  mound  kinetic  energy  Is  actually  much 
closer  to  the  near-optimum  coral  41.31  and  36/10  events  than  to  the  "deep" 
equivalent  overburden  analog  (coral  55.11).  The  great  Importance  of  the 
Interface  spall  reflection  and  lower  layer  spall  velocities  Is  quite  evident 
here,  particularly  for  the  "deep"  water  layer  event. 

Kinetic  energies  In  the  upper  region  or  water  layer  for  the  "coral" 
calculations  are  displayed  In  Figure  111b.  The  upper  section  may  be  thought 
of  as  a  high-velocity  spallation  layer,  most  of  which  enters  freefall  upon 
return  of  the  free  surface  reflection.  In  this  case,  the  discrepancies  between 
underwater  and  equivalent  overburden  events  are  less  evident  than  in  Figure  111a. 
The  "water  layer"  kinetic  energies  rise  more  gradually  and  achieve  the  "free- 
fall"  (constant)  value  later  than  the  equivalent  overburden  events  because 
of  the  greater  thickness  and  lower  transmission  velocity  of  the  water  layer. 

The  final  kinetic  energies  are  slightly  higher  for  the  underwater  events,  a 
fact  which  may  be  partly  attributed  to  the  zero  strength,  low  Internal  energy 
content,  and  efficient  shock  transmission  of  the  fluid  material. 

Energy  curves  for  basalt  36/10  and  basalt  39.83  are  plotted  In  Figure  112. 
In  this  case,  the  kinetic,  internal,  and  cavity  energies  are  all  shown.  Kinetic 
energy  In  the  lower  layer  is  slightly  greater  for  the  underwater  event  "36/10" 
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than  for  ’’39. 83",  although  the  discrepancy  is  much  less  than  for  coral. 

The  Internal  energy  In  the  lower  layer  is  slightly  smaller  for  basalt  36/10, 
again  indicating  more  efficient  conversion  to  kinetic  energy  and  better 
material  ejection  in  the  lower  region  for  the  underwater  case.  The  basalt 
36/10  water  layer  achieves  freefall  spallation  later  than  the  upper  layer 
of  basalt  39.83;  Its  final  kinetic  energy  Is  higher  and  the  internal  energy  is 
lower,  as  expected  for  a  fluid.  The  cavity  energy  content  for  basalt  36/10 
drops  below  basalt  39.83  after  the  Interface  reflection  Impinges  on  the 
cavity,  and  remains  lower  at  late  times.  Thus,  transfer  of  cavity  energy 
to  rock  Is  more  efficient  and  the  gas  acceleration  phase  Is  stronger  for 
basalt  36/10,  as  noted  previously. 


Late-time  energy  data  abstracted  from  the  equivalent  overburden  cal¬ 
culations  are  sunmarlzed  in  Tables  8  and  9  (see  Chapter  VU  with  accoi*>any1ng 
discussion).  Comparing  the  appropriate  pairs  of  calculations,  we  see  that 
coral  41.31  places  less  kinetic  energy  In  both  "lower"  and  "upper"  regions 
than  coral  36/10  does  In  the  corresponding  rock  and  water  layers  (Table  8). 
Hence,  the  ratio  of  total  KE/total  energy  Is  also  less  for  coral  41.31.  The 
kinetic  energy  partition  fractions  between  lower  and  upper  regions  show  that 
relatively  more  of  the  available  kinetic  energy  goes  Into  the  upper  region 
and  less  Into  the  lower  region  for  coral  41.31.  Comparing  coral  55.11  with 
co“al  36/36,  we  find  that  the  discrepancy  between  "lower  layer"  kinetic 
energies  Is  still  greater;  the  coral  55.11  "lower  layer"  energy  Is  less  than 
half  of  the  coral  36/36  energy  In  the  same  region.  However,  the  upper  region 
or  water  layer  kinetic  energies  In  the  two  cases  are  about  the  same.  The 
ratio  of  total  KE/total  energy  Is  obviously  smaller  for  coral  55.11  than  for 
"36/36",  and  the  partition  fraction  shows  that  much  less  of  the  available 
energy  goes  Into  the  lower  region  and  more  Into  the  upper  region  for  coral 
55.11.  These  observations  verify  and  extend  previously  noted  trends:  The 
"equivalent  overburden"  coral  calculations  give  less  kinetic  energy  In  the 
lower  region  (equivalent  to  the  rock  layer)  and  less  total  mound  kinetic 
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energy  than  the  corresponding  underwater  events.  Relatively  more  of  the 

available  energy  Is  partitioned  Into  the  upper  region.  Discrepancies  are 

worst  for  "deep"  layer  events  (l.e.,  a  thick  water  or  overburden  layer 

comprising  a  significant  portion  of  the  entire  problem).  The  "equivalent 

overburden"  events  thus  simulate  detonations  at  a  relatively  deeper  scaled  burial 

depth;  the  Inefficiency  of  kinetic  material  ejection  from  these  events  will  be 

further  accentuated  by  additional  rock  fallback  material  from  the  "upper"  rock 
region. 

Kinetic  energies  for  the  basalt  36/10  and  39.83  calculations  (Table  8) 
show  the  same  similarity  of  behavior  mentioned  earlier.  The  late-time  kinetic 
energies  In  the  rock  layer  or  lower  layer  are  approximately  equal,  while  the 
water  (upper)  layer  energies  and  total  energies  are  slightly  greater  for  the 
underwater  event.  This  similarity  would  probably  break  down  for  deeper  over¬ 
burden  layers,  and  the  underwater  event  would  again  be  expected  to  attain 
relatively  greater  kinetic  energies  In  the  lower  region. 

The  ratios  of  kinetic  energy  to  Internal  energy  ( KE/ IE)  In  the  lower 
and  upper  regions  are  listed  In  Table  9.  The  KE/IE  ratio  of  the  lower  or 
rock  layer  Is  less  for  coral  41.31  than  for  coral  36/10,  and  much  less  for 
coral  55.11  than  for  coral  36/36,  again  Indicating  Inefficient  material  ejection 
In  the  "equivalent  overburden"  cases.  Note  that  the  ratio  for  coral  55.11  Is 
less  than  half  that  for  coral  36/36  or  any  of  the  other  coral  events,  and  Is 
actually  comparable  with  the  "rock  layer"  ratio  for  the  granite  36/36  event 
(low  spall  velocities,  below  optimum  burial  depth).  This  ratio  Is  quite  use¬ 
ful  for  defining  the  overall  dynamic  behavior  of  the  mound,  and  particularly 
In  determining  the  ransltlon  from  efficient  kinetic  ejection  (high  KE/IE 
ratio),  to  predominant  crushing  and  mounding  of  material  (low  KE/IE  ratio). 

The  lower  layer  KE/IE  ratios  for  basalt  39.83  and  basalt  36/10  are  closely 
comparable.  In  the  upper  or  water  layer,  KE/IE  ratios  are  always  smaller 
for  the  equivalent  overburden  events,  although  the  discrepancy  Is  least  for 
the  basalt  calculations.  The  ratio  Is  less  meaningful  In  the  freely-launched 
spall  region,  since  all  material  has  high  ejection  velocities  and  very  low 
Internal  energy  content. 
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In  stannary,  significant  dynamic  differences  do  occur  between 
two-layer  underwater  events  and  the  equivalent  overburden  calculations. 

The  differences  are  greatest  for  thick  overburden  layers.  In  most  Instances, 
the  underwater  event  may  be  expected  to  be  relatively  more  efficient  for 
cratering.  The  discrepancies  between  two-layer  and  equivalent  overburden 
configurations  are  not  directly  correlated  with  the  degree  of  Impedance 
mismatch  at  the  Interface,  and  material  properties  play  an  Important  role  In 
defining  multilayer  Interactions.  Calculations  are  required  to  accurately 
assess  overburden  effects  for  each  specific  combination  of  media  and  layer 
depths.  The  easily-accomplished  "SOC"  vertical  calculations  are  very  useful 
for  determining  relative  layer-depth  effects  and  particularly  for  establishing 
"optimal"  or  "non-optlmal"  configurations  In  a  medium  of  known  properties. 

XI.  CONCLUSIONS 

Hydrodynamic  code  calculations  provide  a  viable  and  informative  approach 
to  the  problems  of  shock  propagation  and  material  dynamics  In  underwater 
cratering.  The  effects  of  water  overburden  and  layer  depth  on  the  overall 
mound  motion  may  be  defined  by  one-dimensional  vertical  calculations.  Indeed, 
such  calculations  very  accurately  depict  the  entire  time-dependent  dynamics 
of  the  underlying  rock  lqyer  throughout  the  period  of  Interest  for  crater 
formation.  High-velocity  spall  launch  Is  the  predominant  material  ejection 
mechanism  near  optimum  burial  depth  In  saturated  media.  The  tensile  waves 
reflected  from  the  Interface  and  free  surface  control  material  dynamics. 

Multiple  reflection  effects  and  detailed  dynamics  of  the  water  layer  have  little 
Influence  on  the  rock.  Thus,  one-dimensional  calculations  are  particularly 

appropriate  for  cratering  events  in  an  underwater  environment.  Likewise, 
peak  shock  parameters  and  safety-related  effects  In  the  water  layer  may  be 
predicted  on  the  basis  of  selected  one-dimensional  results  In  combination  with 
an  acoustic  analysis.  Two-dimensional  code  calculations  are  required  to 
determine  final  crater  dimensions,  detailed  particle  dynamics  within  the 
water  layer,  and  very  late-time  Interactions  for  deeply-buried  multilayer 
configurations. 
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In  view  of  the  progress  achieved  during  this  study,  it  is  suggested 

that  two-dimensional  comparison  calculations  for  one  of  the  media  and 
burial  depth  combinations  used  here  might  prove  worthwhile.  Understand¬ 
ing  of  the  multilayer  cratering  process  has  been  advanced,  and  the  ability 
to  perform  underwater  cratering  design  calculations  ^or  cohesive  rock  media 
is  within  reach. 
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Appendix  C 


Waveforms  for  the  Vertical  and  Vertical 
Underwater  Calculations  ir  Coral 
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FIGURE  C-1.  Velocity- waveform  as  a  function  of  tine  for  coral  calculation  COR  36/00 
Range  R  ■  3.01  ■. 


T  (msec) 


OR  36 


"V*;  XQR  34/05,  *  ■  6.01  *,  with  rezoning  after  Shock  crosses 
Interface  (linear  plot). 


T  (msec) 


FIGURE  C-26.  "V*;  COR  36/05/ R-6.Q1  «,  without  dezoning, constant  zwne  size; 

(linear  plot). 
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Appta^ix  B 


Wave forma  for  the  Vertical  and  Vertical 


Underwater  Calculation*  In  Granite 
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Appendix  F 


Hater  Layer  Waveform 
for  Vertical  Underwater  Calculations 

(Selected  Examples) 
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